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SUMMARY 


A computer program is developed to solve the compressible, laminar boundary-layer 
equations for two-dimensional flow, axisymmetric flow, and quasi-three-dimensional flows 
including the flow along the plane of symmetry, flow along the leading-edge attachment line, 
and swept wing flows with a conical flow approximation. The finite-difference numerical 
procedure used to solve the governing equations is second-order accurate. The flow over a 
wide range of speed, from subsonic to hypersonic speed with perfect gas assumption, can be 
calculated. Various wall boundary conditions, such as wall suction or blowing and hot or 
cold walls, can be applied. 

The results indicate that this boundary-layer code gives velocity and temperature profiles 
which are accurate, smooth, and continuous through the first and second normal derivatives. 
The code presented herein can be coupled with a stability analysis code and used to predict 
the onset of the boundary-layer transition which enables the assessment of the laminar flow 
control techniques. A user’s manual is also included. 
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NOMENCLATURE 


A, B 
a, b 

c 

c* 

c 

c fx 


Cfy 


stagnation point velocity gradients 

major and minor semiaxis lengths ol the ellipsoid of revolution 

pp/p e p c 

B/A 

Pel P or chord 

skin friction coefficient in the .r-direction based on the edge condition, 
Eq. (83) or Eq. (84) or Eq. (S5a) 

skin friction coefficient in the ^-direction based on the edge condition, 
Eq. (85b) 


C p Pressure coefficient 

c p specific heat 

E HI H e 

F ujUe 

k F 

G V jVref Oit Vy f\ rej 


9( G 

H total enthalpy 

h\,h 2 metric coefficients in the x and y coordinates, respectively, 

z, jf, k indices in the ,r, y , and r direction, respectively 
K coefficient of thermal conductivity(= c p fi/ Pr) 

K \ , K 2 geodesic curvature of the curves y = const . and x = const., respectively, 

Eq. (5) 

K\ 2 , K 2 i parameters defined in Eq. (6) 

/ length scale (to be used by the stability analysis codes) 

M Mach number 

mi,.., mis coefficients, defined in Eqs. (55), (61), (63), (68), (70) 
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N 

n u n 2 ,n 3 

P 

Pr 

r 

q w 

s 

T 

U, V , IU 

Vy 

Vye 

V 

x,y,z 

x,y,z 

/ I 

x ,y 
X 


logarithmic disturbance amplification ratio 
coefficients defined in Eqs. (55), (61), (63), (68), (70) 
pressure 
Prandtl number 

radius measured from the center of rotation for axisymmetric body (Fig. 1 ) 

heat transfer at the wall, Eqs, ( 89 )- (9 1 ) 

arc length measured along x direction (from x = 0 ). 

temperature 

velocity components in the x, y, and z directions 

dv/dy 

dvjdy 

total velocity 

bound ary - 1 ayer coord i n ates 

Cartesian coordinates for wing definition, Figs. 3, 4 , 5 . 
reference rectangular coordinates 

distance measured along the chord line {X/c — 1 at the trailing edge) 


a angle of attack 

7 ratio of specific heat 

Ax, A?/, A^grid spacing in the x,y,C directions, respectively. 

6 boundary-layer thickness; ( 2 ) 17 ^= 0.995 

S* displacement thickness in x-direction, defined in Eq. (87) or ( 8 Sa) 

6* displacement thickness in ^-direction, defined in Eq. ( 88 b) 

( transformed normal coordinate, Eq. (40) or Eq. (43) or Eq, ( 45 ) 

6 angle between x and y coordinates 

k A({k + 1)/A ((A’), or curvature 

A sweep angle 


viii 


I 


leading-edge sweep angle 


A, 

A 2 trailing-edge sweep angle 

p viscosity 

v kinematic viscosity, p/p 

p density 

<f> azimuthal angle, 0 and 7r on the windward and leeward plane of symmetry, 

respectively, see Fig. 2 

or stretching variable for an airfoil, cos -1 (l — X/c), Fig. 6 

subscript 

aw adiabatic wall 

e edge of the boundary-layer 

le leading-edge 

s stagnation point or stagnation line (leading-edge attachment line) 

t total 

w wall 

x partial differentiation with respect to x 

y partial differentiation with respect to y 

( partial differentiation with respect to ( 

oo free stream 

superscript 

' ordinary differentiation with respect to ( 
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PART I. NUMERICAL METHOD 


1. INTRODUCTION 

Over the past decade, the major emphasis in computational fluid mechanics has focused 
on numerically solving the Euler and Navier-Stokes ecpiations for increasingly complex aero- 
dynamic shapes. Navier-Stokes solutions, however, are generally much more expensive to 
obtain in terms of computer resources than the boundary-layer solutions, and while capable 
of simulating the physics of complex flows, they are often of low resolution due to grid point 
restrictions. Therefore, Navier-Stokes solution may not be adequate for the stability analy- 
sis which requires accurate solution inside the boundary layer. Furthermore, Navier-Stokes 
solutions are not essential for many design and analysis procedures. 

Renewed emphasis on drag reduction [1], laminar flow control [2], and transition pre- 
diction [1] has indicated the need to develop boundary-layer software that can be routinely 
applied to aerospace vehicles at a fraction of the cost associated with solutions obtained from 
the thin-layer Navier-Stokes equations. Research at the NASA Langley Research Center has 
resulted in the development and verification of three-dimensional boundary-layer procedures 
[3] [4] for application to aerospace configurations. 

In preliminary design, especially for the design of the laminar flow control wing, empe- 
nage, nacelle, or fuselage, it is convenient to use simplified two-dimensional or axisymmetric 
or quasi-three-dimensional flow calculations to limit computation expense. A boundary-layer 
code for the application to swept-back and tapered wings developed by Kaups-Cebeci [5] has 
been used for the design of a laminar flow control wing. This boundary-layer code, however, 
can only solve subsonic and swept-back wing flow with a spanwise conical flow approxima- 
tion. Also, this code yields discontinuous first and second normal derivatives of the velocity 
along with an oscillatory skin friction coefficient for some flow cases, especially when suction 
is applied. 



In this report, an efficient, second-order accurate, finite-difference method is used to 
solve the two-dimensional or axisymmetric or quasi-three-dimensional, compressible, lami- 
nar boundary-layer equations. The quasi-three-dimensional flows considered in the present 
report include flow along the plane of symmetry, flow along the leading-edge attachment 
line, and swept wing flows with conical flow approximations. The swept wing case includes 
swept-back tapered wings and swept-forward tapered wing with a. spanwise conical flow as- 
sumption, and swept-back delta (low aspect ratio) wings with a streamwise conical flow 
assumption which is more appropriate for the flow over a low aspect-ratio wing. 

Results are presented for several test cases. This method was used for the boundary-layer 
stability analysis of axisymmetric laminar-flow-control nacelles [6] and of the flow along the 
leading-edge attachment line [7]. The method is valid for perfect gas flows from subsonic to 
hypersonic Mach numbers. Interaction between the inviscid and viscous flow is not included. 
A user’s manual with a detailed description of the computer program for the boundary-layer 
analysis is presented in PART. II. 
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2. BOUNDARY-LAYER EQUATIONS 


Three-dimensional, compressible, laminar boundary-layer equations (in dimensional form) 
for the general nonorthogonal body-oriented coordinate system, from which the governing 
equations for two-dimensional, axisymmetric, and quasi-three-dimensional flows are derived, 
are as follows [3], [8], [9]: 


continuity equation 


d c) c) 

~{pu h 2 sin 0) + —(pvhi sinO) + — (pwhili 2 sin 0) = 0 
Ox ay dz 


^-momentum equation 


pu du pv du 

1 r\ ”f" j 


+ pw — pu 2 I\ \ cot 0 + pv 2 K 2 esc 0 + puvK\ 


h x dx h 2 dy f dz r 1 1 A 

esc 2 0 dp cot 9 esc 0 dp d . du , 

= 7 — 7T + 7 -f + 

hi dx h 2 dy dz dz 


^/-momentum equation 


pu dv pv dv dv 2 2 

~r~ o — I" + p^o— pv A 2 cot 0 + pu Ai esc 6 + puvh 2 \ 

h\ dx h-2 dy dz 

cot 9 esc 0 dp esc 2 8 dp d dv 

= hi fa ~ h 2 dy + dz^dP 


energy equation 


pu dtl pv OH dH dip dH 


1777 + r 2 ~a^ + 77 " 


dz Pr dz 


+ M 1 “ 771 


1 . d ,V 7 


Pr dz 2 


The metric coefficients h\ and h 2 are functions of x and y. Theta, 6 , is the angle between the 
x and y coordinates. The parameters f\\ and I\ 2 are the geodesic curvatures of the curves 


y = const and x ~ const respectively, where 


d n m dh 1 

— (h 2 cosO) - — — 


h\ h 2 sin 9 1 dx 


dy r 


I<2 = 


h\h 2 sin 0 1 dy 


0 ,, dhi 

— (ft] cos 0) - -7— - 


A', 9 = 


A 1 /i 2 sin 2 0 


(1 + cos 2 9)- 


2 cos 9- 


h 1 /) 2 sin 2 


7-2-pl + cos 0)— -2cos<?— 
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V is the total velocity and is given by 


V = ( ft + v 2 + 2uv cos ft) 


1/2 


( 7 ) 


The boundary conditions are 


= 6, u = u e (x,y), v = v e (x,y), II = II c 


dH 


2 = 0, u = v = 0, tv — w w , II = II w or ( — ) W = Q 


(8-a) 

(S.b) 


The pressure gradients are related to the inviscid velocities by the following equations at the 
edge of the boundary-layer: 


du, v e du 


€ ( yu e t'e yj 2 


+ 


— ttp A i cot 0 + v: I \2 esc ft -f- u e v t K\ 2 


h 2 Oy 


[hi dx h 2 dy 
esc ft dp cot ft esc ft dp 
hi dx ^ 

[ u t dv e v e dv e -2 r- /I 2r- n . I' 

{ h tip A 2 cot ft + Up A ! CSC ft + U e V e h 21 

( /t, a.r ft 2 dy 

cot ft CSC ft dp esc 2 ft dp 


(9. a) 


(9 b) 


ft ! dx hi dy 

The perfect gas equation of state and Sutherland’s viscosity are used to close the equation 
set. 

2.1 Stagnation Point 

To obtain the boundary-layer solutions at the three-dimensional stagnation point, the 
governing equations for three-dimensional laminar compressible flows in Cartesian coordi- 
nates are required and can be obtained by setting h\ = 1, h 2 = 1, and ft = 7r/2 in Eqs. 

(l)-(4). 


continuity equation 


x-momentum equation 


J-M + + £(rw) = o 


( 10 ) 


du du du dp d du 

pU ^ + pl % + PU Tz = -Vx + ¥z {fl d-z ) 


( 11 ) 


i 


y-momentum equation 


dv dv dv dp d dv 

pu ^ +p % +fn, Tz = -ar J + rJ l % ) 


energy equation 


dH dH dH 

pU -^ + pV ^ + pW ^ = 


d ( p dH 
dz\Pr dz 


+ y(l — 


Prr >2 2 ' 


where 


V = {u + u) 1/2 


( 12 ) 


(13) 


(14) 


The boundary conditions are 


z = 6, u = u e {x,y), v = v e (x,y), H = H e (15. a) 

d H 

z — 0, u = v = 0, w = w w , H = H w or (-r— )„, = 0 (15. b) 

oz 


2.2 Two-dimensional and Axisymmetric Flow 


The boundary-layer equations for two-dimensional and axisymmetric flow can be obtained 
from Eqs. (l)-(4) by substituting 0 = x/2, v = 0, h.\ = 1, h 2 = r 3 (j = 0 for 2-D flow; j — 1 
for axisymmetric flow). See Figure 1 lor the coordinate system for 2-D and axisymmetric 
flow. 

continuity equation 

~—{pur J ) + ~^-(pwi J ) — 0 (16) 

dx dz 

^-momentum equation 


du dv dp d du 

pu— + pw— - -tt- + zz-{Pzz-) 
dx dz dx dz dz 


(17) 


energy equation 


dH 


dH 


pu ^ + pw -dl 


ft du 
p 7- dz 


+ /*(i 


Pr'rh 2 ' 


(18) 


where 


V = u 
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( 19 ) 



The boundary conditions are 


(20. a) 

(20. b) 



u = u e (;r), II = H e 


u - 0, ir = w u 


H = H„ 


or 


,dH s 

{ dz )w ~° 


The pressure gradients are related to the inviscid velocity by the following equation at the 
edge of the boundary-layer: 


d u e 

PeUc-r- 
a x 


dp 

dx 


( 21 ) 


2.3 Flow along Plane of Symmetry and Along Leading-Edge Attachment Line 


Along the windward and leeward plane of symmetry and along the leading-edge attach- 
ment line, v and Ah esc 0 are zero and 0 is generally 7r/2.( 0 is retained for the general system) 
Consequently, each term in the {/-momentum equation(Eq.(3)) vanishes. However, partial 
differentiation of Eq. (3) with respect to y yields an equation for dv/dy. After differentiation 
and using the appropriate symmetry conditions (du/dy = dw/dy = d 2 v/dy 2 — dH/dy = 
dhx/dy = dh^jdy = 0) along with Eq. (9), the governing equations for the flow along these 
lines become (See Figure 2 for the coordinate system.): 
continuity equation 


d d 

-r— (/)«/* 2 sin 0) + pvyh i sin 0 + — {pwh\h 2 sin 6 ) = 0 
dx UZ 


( 22 ) 


x-momentum equation 


pudu du , , . „ u e du e 9 r , d , du, , 

T, Tx + pw F- ~ pu A ' cot 0 = pA T, SI - “* A 1 cot 0) + Tz^ (23) 


{/-momentum equation 


pu Ov 


h\ Ox 


5? + P w % 7 + {~ 2 v l + P ul, y + pu 2 


djKi esc 0) 

dy 


He dv V e l J ti 


/?i dx 


-»e '-'vy e ye j’ \ i 2^(^1 CSC 0) d dl)y 

Pei 7 — + T~ + UeV y eI\2l) + PtK 7TT + 


dy 


(24) 


energy equation 


pudll dH d ( p dll , Id ,V 2 ' 

T, lh + pw & = Fz FJI + " (1 “ 7v>fo ( T> 


(25) 


I 


where v y = dv/dy , v ye = dv e /dy and V = u along these lines. 
The boundary conditions for the plane of symmetry are 


= <5, 


z = 0, 


U = K e (.r), V y = V ye , // = H e 

u = v = = 0, u) = ie u ,, II = //,„ or 




(26. a) 
(26. b) 


2.4 Swept Wing Flow 


The governing equation for a three-dimensional compressible laminar flow in a polar 
coordinate system with conical flow assumption ( dp/dx = 0) can be obtained by substituting 
h\ = 1, h 2 = x , 6 = 7r/2 from the equations (l)-(4). Figures 3, 4, and 5 show the coordinate 
systems used for the swept-back wing with the spanwise conical flow assumption, the swept- 
forward wing with the spanwise conical flow assumption, and the swept-back delta (low 
aspect ratio) wing with the streamwise conical flow assumption, respectively. Fig. 6 shows 
the airfoil definition used for the swept wing cases, 
continuity equation 

o o o 

^ (/>«*) + ^(/w) + -qz(P wx ) = 0 ( 27 ) 

x-momentum equation 


da pv da da pv 2 d da 
P u ir_ + — -57. + P w ~nZ - — = 7 r(Pir> 


dx x dy 


dz 


X 


dz Kr dz 


(28) 


y - momentum equation 


dv pv dv dv puv 

pu-r~ + — tt- + P«J— + 

ox x ay dz x 


l dp d . dv 
x dy dz dz 


energy equation 


dH pv dH 


OH 


-£{ 


dip dH . 1.0 ,V 2 ' 

dz\Prdz + ^ Pr dz 2 


where 

V = (u + t>)' /2 


(29) 


(30) 


(31) 
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The boundary conditions are 


z = 6, u = u e (y), v = i > e (y), H = H, 

z = 0, 


dll 


u — v = 0, to = w w , H = H w or (-rr-)u» = 0 

02 


(32. a) 
(32. b) 


The pressure gradients are related to the inviscid velocities at the edge of the boundary-layer 
by the following equations: 

du. 


d y 

dy 


2 dp 
Pedy 


(33. a) 
(33. b) 


Along the stagnation line, t> = 0. Partial differentiation of Eq. (29) with respect to y 
yields an equation for dv/dy. After differentiation and substituting v = 0, the governing 
equations for the stagnation line become: 
continuity equation 


x-momentum equation 


c) d 

— (pu.r) + pv y + Q^(pwx) = 0 


du du d du 
pu— + pw— = —(/< — 
dx d: dz dz 


y-momentum equation 

di 


P .y dl'y OUVy 1 d 2 p d dv y 

pu-TT + ~ v y + P w ^r + — = — tt + V 

Ox x y dz x x dy 2 dz dz 


energy equation 


pu— 


du du d ( p dH 

d7 + pw d= ~ dz l Pr dz 





where 


V = 


u 


The boundary conditions are 


z = 6, 
z = 0. 


U - U e , Vy = Vyg , 


II = H e 


dll 


u = v y = 0, w = w w , II = H w or ( )u< = o 

dz 


(34) 

(35) 

(36) 

(37) 

(38) 

(39. a) 
(39. b) 
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3. TRANSFORMED EQUATIONS 


The boundary-layer equations are transformed to a coordinate system which allows a 
coupled solution of the continuity and momentum equations. In this report, three different 
transformations are used. 

(A) The following transformation is used for stagnation point, 2-D flow, axisymmetric. 
flow, and flow along the plane of symmetry. 


where 


x = x , y = y, 




h i dx 


together with the following definitions 


(40) 


(41) 


F = fc = ’ll /lie, G = g < = v/Vref, E = H/H e 


(42) 


For the stagnation point, V re j — v e ; for 2-D and axisymmetric flow, G = — 0; for the 

plane of symmetry, V re j — Ko, and G = = Vy/Voo. This transformation removes the 

singularity at x = 0. 

(B) The following transformation is used for the flow along the leading-edge attachment 
line. 


x = 


;r, 


V = Ih 



(43) 


and the definitions 


F - f ( = ti/tie, G = (M = v y /V^ E - H/H e (44) 


This transformation is used to avoid an infinite transformed normal coordinate (( e ) near 
x = 0; such an occurrence is possible when using the transformation given in Eq.(40). 

(C) The following transformation is used for the swept wing flow with conical flow as- 
sumption. 


x = ;r, y = y, 



(45) 
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and the definitions 


F — fa — u/ii € , G = Pc — v l v t) & — H/He (off the stagnation line) (46. a) 
F = = u/u e , G = </( = Vyfvye, E — 11/ II t (on the stagnation line)(46.b) 

This transformation is used to handle cases when u e is positive or negative. 


3.1 Three-Dimensional Stagnation Point 

The governing equations for three-dimensional laminar compressible flows in rectangular 
coordinates, Eqs. (10)-(13), are transformed using Eq. (40) and Howarth’s [10] inviscid ve- 
locity components near the stagnation point. For s approaching zero, the following ordinary 
differential equations are obtained. 

^-momentum equation 

( c.f")' + If" - u'f + 4 f"g + — = 0 (47) 

A P 

y-momentum equation 

{(■</)' + !</ - -Ag’f + -J'g + ~ = 0 (48) 

/I A A p 

energy equation 

(Ee’Y + V + ?-g)E' = 0 (49) 

Fr A 

The equations above are based on the assumptions that the outer flow is irrotational and 
that the inviscid velocity components near the stagnation point can be approximated by 

tf ( = Ax, v e = By (50) 


The primes denote ordinary differentiation with respect to (\ i.e., 


(If u , dy v 


d( M e ’ 


(J <K V e ’ 
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The boundary conditions are 


(52. a) 
(52. b) 


C = 0 : / = /' = g = g' = o, E f = 0 or E = E w 

C = Ce : /'=!, = £ = 1 

Two-dimensional and axisymmetric stagnation point equations can be obtained by sub- 
stituting C* = B/A — 0 (for a 2-D stagnation point) or C* = B/A = 1 (for an axisymmetric 
stagnation point). 


3.2 Two-Dimensional and Axisymmetric Flow 

Using the transformation given in Eq. (40) and the relations given by Eq. (21), the 
governing Eqs. (16)-( IS) are transformed to the following form: 
x-momentum equation 

F = k (53. a) 

(CF c ) ( + mj/Fc - m 2 F 2 + m n c - 77?. 13 F C = m w {FF x - F c f x ) (53. b) 

energy equation 

(??iF f ) c + » 2 F < + (n 3 ) c - 77i 13 F c = m 10 {FE x - E c f x ) (54) 


The coefficients ???i to 77 ? 13 and n\ to n 3 are as follows: 


= « O + 




+ 


5 


/) 1 H e 5 .r J ll 1 /? 2 \fPcPe 9 x 


{ll 2 \/PePe} 


77? 2 = 


77? 10 


3 

/?l« e <9.r 


/?! 

77? H = 7772 


77? 13 = 


(/>'«Q U , I PcUe* 
p ( u t V /'< 


(55. a) 

(55. b) 
(55. c) 
(55. d) 
(55. e) 
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(55. f) 
(55. g) 
(55. h) 


»i = 


C 

~p7- 


«2 
n 3 


mi / 

1 )FF 

U7 (l " 


Note that h\ = 1 and h 2 = r J , where 7=0 ior 2-D flow or 7 = 1 for axisymmetric flow. 
The boundary conditions are 


C = 0 : f = F = 0, E' = 0 or E = E w (56.a) 

C = Ce • ^ = 1, £ = 1 (56. b) 


3.3 Flow Along Plane of Symmetry 

The governing equations lor the plane of symmetry are transformed using the transfor- 
mation given in Eq. (40) and by defining 

F = /< = u/u e , G = gc - Vy/Voo, E = H/H e (57) 

^-momentum equation 

F = /< (58. a) 

(CF c ) c + ?)?i/F c - m 2 F 2 + m 6 F ( <7 + m n c - 777 13 F C = m 10 {FF x - F^/*) (58. b) 

y-momentum equation 

G = yc. (59. a) 

((?(?()( + ? 7 ?i JG^ — 7773G 2 — 777 ,j FG m^G^g — m$F 2 -f 77712c — niisG^ 

= m 10 (FG x — G q f x ) (59. b) 

energy equation 

(" 1 F( )( + 77 2 F^ + (»7 3 )c — 1»13 F<; = 777-io(F E x — E(f x ) (60) 

12 
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The coefficients »??i-n?i 3 , 111-113 are defined as follows: 


= 7: i 1 + 


S 0 U e 


+ 


5 


2 \ h\u e dx j k 1 h 2 sin Oy/p^dx 


{ h 2 sin Os/pefle] 


m-i = 


m 3 = 


777 4 


s du f 


h\u e dx 
sVc o 


— sA'l cot 0 


h 2 Ut 


A 21 + 


* aVo, 
/? iKco a* 


m & = 777.3 

S7i e a(A'i CSC 0) 


7770 = 


777-10 


777 n 


v: 


A 1 


dy 


,s #17 


mi 2 = 


/? 1 U f #.T 

U7 e #( A , esc 0 ) 


#«/ 


mi3 


77 1 = 


v; 


( pw) u , p'U e S 


s cot 0K\ 

+ 77”( 


1 dv ye l j/e . 

,, V, o 1 } rA 2 ii7 ye J 

V, x , h\ Ox U e li 2 


Pt •«< 


c 


Pl- 

772 = 777] / + 777(5 7/ 


,,3= «7 (1 


1 

P^ 


)PA’ C 


The boundary conditions are 

£ = 0 : f=F = g — G = 0, I'j' — 0 or E — E u 

C = Ce : 

where TV = dvjdy. 


E = L G = , £ = 1 


(61. a) 

(61 .b) 
(61. c) 

(61 .d) 
(61.e) 

(61 -f ) 
(61 -g) 
(61. h) 
(61 -i) 
(61 -j) 

(61 .k) 
(61.1) 
(61. m) 


(62. a) 
(62. b) 


3.4 Flow Along A Leading-Edge Attachment Line 

The governing equations tor the flow along the leading-edge attachment line are tians- 

formed using the transformation given in Eq. (43) to the same equations as for the flow 
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along the plane of symmetry (Eqs.(58)-(60)). However, because of the difference in the 
transformation, the coefficients m t to m 13 and to n 3 are defined now as follows: 


1 du t 1 d 

m i = 7T, T) •" T7 — • n {h 2 sm 0^p c Pe) 

2h\ U e OX III 1)2 Sill UsJprJO ox 
1 du * j - . n 

77 ? 2 = Ai cot 0 

h\ U t Ox 

Voo 

= 1 

h 2 iic 

m 4 = K 21 




S(/v, esc 0) 


m 10 = 7“ 

«-i 

1 

?»11 = - cot 0 A 1 

<xr 

« f d(I\'i cscO) I 1,1 dry _ t>J e 

m 12 — . ■ h 77“ ( ; 7 I ] h A 2 lVj,e) 

I'co a?/ Vqo /? 1 ar u c /?2 

(P«')«r I P(Ut 

»»-13 = — \ 

/V«e V /‘e 


n 2 = mi/ + (6 

" 3 = T 77 (l - ~k-^ FFi < 63 

The boundary conditions are the same as for the flow along the plane of symmetry and 
repeated here for completeness: 


C = 0 : f=F = g = G = 0, E' = 0 


E = E„ 


C — Cc : A = l, G’ = -f, £= 1 

* DC) 


if 


3.5 Swept Wing Flow 


The governing equations for the swept wing flow (off the stagnation line) are transformed 


using the transformation given in Eq. ( 45 ) to the following equations, 
x-momentum equation 

F = fc ( 65 . a) 

(CF ( )t + ?7 ?i/F c - m 2 F 2 + m 6 F c <y + m n c - m 13 F< - m 7 {GF y - F ( g y ) ( 65 . b) 

y- momentum equation 

G = g c (66. a) 

{CGq)q + it) i f G^ - m 3 G 3 - m 4 FG + m 6 G^g - )n 9 F 2 + m r2 c - m 13 6'( 

= m 7 (GG y — Gfg y ) (66.b) 

energy equation 

{n\E^)^ + U‘2 E^ + ( ?? 3 )< ~~ — m^{GE y — E^g y ) ( 67 ) 

The coefficients ??q to /?? 13 and ??i to 77.3 are defined as lollows: 




??>,, = 0 


?7?i2 


m 13 = 


1 d v e 
I u e | d. y 
(pw)w 

Pe | Me | \ 


+ 


I Pe | «e I X 


Pe 


c 


n i = 


Pi- 

n -2 = nnf + m 6 <7 


»3 


Cuj 

He 


(1 


1 

P 


-){fF c + GG ( (^) 2 )} 

r l ti e J 


The boundary conditions are 


C = 0 : f = F = g = G = 0, E' = 0 or E = E w 

C = Ce : F = l, G = 1 , E = 1 


(68-j) 

(6S.k) 

( 68 . 1 ) 

(68. m) 
(68. n) 
(68. o) 


(69. a) 
(69. b) 


The governing equations on the stagnation line (leading-edge attachment line) for the 
swept wing flow are transformed to the following equations: 
i-momentum equation 

F = f( (70. a.) 

(CF ( ) ( + mj/F<; - m 2 F 2 + m G F^g + m n c - m 13 F c = 0. (70. b) 


y-momentuni equation 


G — (j(_ (71. a) 
(CG\)( + m 1 / G< - m 3 G 2 - m A FG + m 6 G\g - m 9 F 2 + m 12 c - m 13 G'( = 0 (71. b) 


energy equation 

( E ( + n 2 £< T ( )c “ 777 13 = 0 (72) 

The coefficients ???!, ??? 2 , rn 9 , m n , m 12 , m J3 , n x , ?r 2 are the same as in Eq. (68). 

The remaining coefficients are defined as follows: 

1 d v e 


nii< - 


Me | dy 
16 


(73. a) 


i 


in 7 = 0 

(73,b) 

>n s - 0 

(73. c) 

Cu*. 1 , rr 

(73. d) 
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4. NUMERICAL METHOD 


All the equations listed in the previous section can be expressed in the following form: 


x-momentum equation 

F = / c (74. a) 

(CE c )c + m x f - m 2 F 2 - m 3 FG + m 6 F c g - m s G 2 + m u c - m 13 F c 

= n?io(/ F x — l\jx) + u?? (GFy — F^(jy ) (74. b) 

^-momentum equation 

^ = 3c (75. a) 

(GG<;)<; -f mifG{ — m 3 G 2 — m±FG -f m^G^g — m 9 F 2 -f m^c — m^G^ 

= n?io(EGj- — G^jx) + ???t( GGy — G^r/ y ) (75. b) 

(75x) 


energy equation 


("i£ C )c ?? 2 ^c T ( 7 ? 3 ) c ” ??? i3^c — m io(EVG — F^ fx) -f m7(G£ , y — E^g y ) (76. a) 

(76. b) 


Note that either m? or m 10 is zero for two-dimensional, axisymmetric, and quasi-three- 
dimensional flows. 

The above equations are linearized using Newton-Rhapson’s linearization technique [11]. 
The (-derivative terms are discretized using a central finite-difference scheme. An implicit 
second order backward finite-difference is used for the x- and y-derivative terms.. An implicit 
marching procedure, which is similar method used in Ref. [12], is used to solve the governing 
equations. 


For abbreviation, finite-difference operators are defined as 


AC(A:) ' /l ’ +1 ^'^ACjk-l )G k E’-l, 
AG + A(/,_i 


= 2, 3, .., kmax — 


(77.a) 


IS 


A c (C fc A c ft) = — 


' Fk+i - F k r F k -F k - , 

t-' A:+ 1 / 2 y** ^ ^ ^ 


-1 


AC* 


AC*—i 


f ^ F - Fi_i . f . 0 

FF = r it i=2 

axi 


(77,b) 
(77 -c) 


c r , { Aa*^ — (Aaq -fi ) 2 } t i + ( A x t + A;r t -i ) Fj- i Aa,-F t -2 ^ 3/77 r n 

vr a 1 3 4 \ / * . a \ A /A 1 A — \9 ' 


(A.r,-) 2 (A xi + A®i_i) - Ax^Axi + A.r.-i) 2 


F — F 1 

SyFj = 1 7 — if j=2 


Aj/j 


(77.e) 


Vi 


{Ay] - (Aft + A W -,) 2 } F, + (A gj + Ayj_i) 2 Fj_i ~ A y)Fj-i , f . f) 

(A2/j) 2 (A?/j + Aj/ 7 _i) — At/j ( Aj/j + Ay^i) 2 


where Ck ±\/2 “ |(C* + C7,*±i)i Ac^ — (,a-+i — Go A.t,- — X{ £*- 1 , Aj/j — j/j an ^ G 

j, and A* represent the ;y, and £ directions, respectively. The overlined quantities are the 

converged solution at the previous step (/ — 1, ? — 2, j — 1 , j — 2). 

The fin ite-di (Terence equations lor the Eqs. (74)-(/6) are linearized as tollows: 
x-momentum equation 


fk - /*- 1 - + F,_.) = 0 (78.a) 

A c {C k A ( F k ) + mxif^Fk + 8 ( T k f k ~ 7*V*) " rn 2 (2FF k - F 2 k ) - m 5 (G k F k 
+F k G k - F k G k ) + m 6 (g k 8 c F k + 8<T\.g k - g k 8 c F k ) - m s (2GG k ~G't) + n me* - m l3 F c 
= m U )(F k 8 x F t + 8 r F t F k — F k 8 x F t — 8^F k 8 x fi — 8 (> F k 8 x f • + 6^F k 6 x f t ) 

+m 7 {G k S y F k + 6 y F k G k — G k Syf k — 8 ( 7 k 8 y g k — 8(F k S y g k + 8^F k 8 y g k ) (78. b) 


^-momentum equation 


g k - g fc _! - ^ 1 [G k + Cu-i) = 0 (79. a) 

A((C k A ( 6*) + ni]{f k 8 c G k + 6<G k fk ~ AAG'*) ~ m X2GG k - G k ) — m 4 (G k F k 
+ F k G k - F k G k ) + m 6 {g k 8 c G k + 8 c G k g k - g k 6<G k ) - m 9 (2FF k - F k ) + m 12 c k - m 13 G< 


= mio(F k 8 x Gi + 8 x G t F k - F k 8 x Gi — 8 Cj G k 8 x f l — 8^G k 8 x j ■ + 8 c G k 8 x f t ) 

+m 7 (G k 8yG k + 6„G k G k - G k 8 y G k — 8^G k 8 y g k — 8^G k 8 y g k + 8^G k 8 y g k ) (79. b) 
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energy equation 


&<E k )+n2,k S^Ek + V*3,* ~ m t3£< = m io{EkS x Ei-6 < Ek6 x fi)+m 7 (Gi 1 SyEj-6 < Ek8ygj) 

( 80 ) 

where the overlined quantities are evaluated irom the previous iteration. The energy equation 
does not require linearization since it is solved after the momentum equations. 

The finite- difference momentum equations, Eqs. (78) and (79), are rearranged into a 2x2 
block tridiagonal form as 

hk — hk - i 4 ^ {Hk + Hk- 1) (81. a) 

— Ak-IIk - 1 + Bkl'Ik — Cklh+i + a t;h k — Dk (81. b) 


where 


h k = 


h 

9k 


Ih 


F b 

Ch 


/4fc, Bk, Ck , (ik are 2x2 matrices, and is a. vector. 


These equations are solved by the Davis Modified Tridiagonal Algorithm (See Appendix A). 
The finite- difference energy equation, Eq. (SO), is arranged into a linear tridiagonal matrix 
equation form as 


BkEk- 1 + D k E k + A k Ek+i = Ck (82) 

where At, Bk , Ck, and Dk are scalars. 

This equation is solved using the Thomas Algorithm. The momentum equations and the 
energy equation are solved iteratively in an uncoupled manner until the converged solution is 
obtained. The converged solution is usually obtained within ten iterations. If the converged 
solution is not obtained within twenty iterations, the flow may have been separated. 
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5. RESULTS AND DISCUSSION 


5.1 Boundary-layer Parameter Definition 


The skin friction coefficients are defined and calculated from the following equations. 

(a) for 2-D, axisymmetric flow and flow along a plane of symmetry: 

_ (ndu/d:) w _ 2p w u e ( OF/d ( ) w {p/ p e ) w {p e u e / p, -s) l/2 
L 'J X ~ l n T/2 - n t/2 1 

2 re 1 e re v e 

(b) for the flow along a leading-edge attachment line: 

^ {flOu/dz) w 2n w 'U e (dF'/ d( > ) w (p/ Pe)w{pt^e ! Pe) ^ / 

= l n 1/2 “ n T/2 1 

2 re * c re ’ e 


(c) for swept wing flow: 


^ {pdu/dz 

C-'/r \ n 

ifidv/dz 

fy ~ I/? t » a 

2 ' c 1 f 


2/« w | «e I {dF/dQ w (p/P')w(pt I tf e | 

Pe U l 

2fl u ,V ( (dG/dC)w{p/pe)w{pe | | 


(85. a) 
(85. b) 


The derivative terms; (dF/0 () w . and (■ 0G/dQ w are evaluated by second order one-sided 
differences at the wall, i.e., 

, OF (ACi + AC 2 ) 2 C 2 - (ACi) 2 F 3 /fiR , 


OF (ACi+AC 2 ) 2 F 2 -(AG) 2 F 3 

( <?C A (A Ci + ACi ) 2 - (AC, + AC 2 )(AC ,) 2 

OG _ (ACi + AC2) 2 C2 — (A(,i) 2 G’ 3 

( ^c )u ’ “ ACi(ACi + ACi ) 2 - (A Ci + ACi)(ACi ) 2 

Displacement thickness is defined and calculated from the following equations, 
(a) all the flows except the swept wing flow: 


(86. a) 

(86. b) 


/"< i 

Jo 


(b) for swept wing flow: 



fd - 

pit 


Jo 

Pe 

II 

* i 

< -o 

Hi- 

p r 


Jo 

Pe V e 


J Frr)i: = fd - ~F),h 

Pe V'e JO Pe 


)dz= (1 -J-F)d: 

Jo p e 



Momentum thickness is defined and calculated in a similar manner as the displacement 
thickness. 


Heat transfer is calculated from: 

(a) for 2-D, axisynnnetric flow and flow along plane of symmetry: 


t w O v / w \ 

ay rr p e p e s 

(b) for the flow along leading-edge attachment line:: 




ay rr p e p e oc. 


(c) for swept wing flow: 


((w = A'i 


dT _ ( P /'»■■ f P V , Pe I «e \ . U2 dT 
■a )n> — n { )w\ ) l .-)> )b 

ay rr p f p € x 


where (dT/d() w is obtained from: 


, > _ [(Ac,,) 2 ~ (ACi + AC 2 ) 2 ]7\ + (ACi + AC 2) 2 ^2 — (ACi ) 2 T’ 3 

1 dC )w ~ AC, (AC, + AC *) 2 - (AC, + ACCHACi ) 2 


(89) 


(90) 


(91) 


(92) 


5.2 Test. Cases 

5.2.1 Two- Dimensional Flow 

The subsonic flow past a NACA 0012 airfoil at 0 degrees angle of attack was selected as 
a two-dimensional flow test case with the following flow conditions (This flow condition is 
the same as that, used in Ref. [4]): 

A /00 = 0.5 
a = 0° 

Pot, = 2116 lb/ ft 2 

7^ = 520 °R 
T = T 
c = 0.28266 ft 

The free stream Reynolds number based on the chord length is 10'\ 
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The inviscid solution was obtained using the Euler code developed by Drela [13]. Figure 
7 shows the skin friction coefficients. In this figure, skin friction based on the free stream 
velocity (C/oo = ) is plotted to compare the results with Ref. [4]) These results were 

obtained using uniform grid spacing in the normal direction with A£i = 0.1 and £ e = 7-0. 
The inviscid grid was used in the ^'-direction. The results are in good agreement with the 
result of Ref. [4] (not shown). 

5.2.2 Axisymmetric Flow 

The hypersonic flow over the sharp cone (half cone angle of 5 degrees) with mass transfer 
at the wall was selected as an axisymmetric flow test case with the following flow conditions 
(This flow condition is the same as that used in Ref. [4] and [14]): 

Moo = 7.4 

half cone angle =5° 

cv = 0° 

Poo = 701.4 N/m 2 
^ = 69.7 °K 
T w = 316.65° 

From X = 0 to X = 0.096 m, there is no mass transfer at the wall. From X = 
0.096 m, two types of mass transler exist: {pw) w — 0; ( pw) w = —0.0901 1 7 N scc/m 3 or 
C q = -0.0020754 (wall suction). 

The inviscid pressure was obtained using the Euler code developed by M. D. Salas (un- 
published work). The boundary-layer edge velocity (u e /V x .) from the Euler code is about 
0.99176 and the initial velocity {u e /V c 0 at i=l) calculated using Eq. (Bl) is 0.99184. There- 
fore, for this flow, the initial edge velocity calculated assuming isentropic condition is very 
close (within 0.01 percent difference) to the inviscid velocity given by the Euler code. This 
seems to be because the shock angle is small (for this flow, about 9.2 degrees). 

Figure 8 shows the skin friction coefficients on the cone with and without wall suction 
condition. These results were obtained using non-uniform (stretched) grid spacing in the 



normal direction with A(i = 0.005 and A((k + l)/£((k) = 1.05. The step-size in the 
x-direction (A.r) varies from 0.001 (near the nose) to 0.01. The grid distribution used to 
obtain Fig. 8 can be seen in PART II, section 6.3. The result is in good agreement with other 
result [4] (not shown). 


5.2.3 Flow along Plane of Symmetry 

An ellipsoid of revolution having a four to one ratio of major to minor axis (a — 1 m, 
b = 1/4 777 ) was selected as a test case, and the boundary-layer solutions were obtained for 
incompressible flo w(A/ OC) = 0.05) at a = 6°. 

For this body, the angle between the x and y coordinates ( 0 ) is w/2. The metric coef- 
ficients can be obtained exactly, and the velocity components can be obtained analytically 
for the incompressible flow [15]: 


b 2 
u. e 


/ 1 + (A~/» — 1) 2 (/ 2 — 1) ) 1/2 
1 1 - (X/a - l) 2 j 

byft - (X/a - l) 2 

( Vo ( f ) cos o' cos ft — 14)o {t) sin o sin ft cos d>) 


(93. a) 

(93. b) 
(93. c) 


i’ f = Kx.(l4 J0 (/) sin o sin <f>) 


(93. d) 


where t = b/a. Note that r is measured along the axes X although .T-direction is along the 
body surface; therefore .v = A in quantity. Here ft is the angle between the line tangent to 
the ellipse and the positive A axis; it is given by 


cos ft 


ft - (A 7« - 1)2 


ft + (. Y /« - 1) 2 (< 2 - 1 ) 

ft < 0 if X/a > 1, and ft > 0 if X/a < 1 


(94. a) 
(94. b) 


The parameters Vo ( t , ) and are functions of t and are defined by 

(1 - t 2 ) 3 / 2 


Vo(t) 






2Vo(0 - 1 
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(95. a) 


(95. b) 
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The skin friction coefficients (C fooVR e oo = yv'' where a — 1) as a function 

of X at an angle of attack of 6 degrees are shown in Fig. 9. The present numerical results 
were obtained using the following grid distributions: Ar=0.002 near the nose followed by 
Ax=0.02 downstream, A y = 5°, and 61 grid points normal to surface uniformly spaced such 
that A£=0T. The result is in good agreement with result given by Iyer [4] (note that Iyer’s 
result must be multiplied by l/\/2 for comparison due to the different definition of Re 0 o). 

5.2.4 Flow Along A Leading-edge Attachment Line 

The subsonic flow past an infinite swept cylinder with a diameter of 0.75 feet and a sweep 
angle of 60° a.t 0 degrees angle of attack was selected as a test case with the following flow 
conditions: 

A/oo = 0.126675 
A = 60° 
a = 0° 

Poo = 2116.8 lb/ ft 2 
= 520 °R 
T w = T aw 

The inviscid solution was obtained from the analytical potential solution along with 
the sweep theory. For this flow, the ^-coordinate is along the leading-edge, and the y- 
coordinate is normal to the leading-edge and measured from the center of cylinder. The 
metric coefficients are defined h\ = 1.0 and = 0.375 where h 2 is the radius of the cylinder. 
The skin friction coefficient is the same along the leading-edge attachment line. The velocity 
profiles at various suction levels are plotted instead of the skin-friction coefficient for this 
test case. Figure 10 shows the results of applying suction C q = 0 , —0.0005, —0.001, —0.002, 
—0.005, and —0.01. The Reynolds number based on the momentum thickness at the zero 
suction condition is about 203. 
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5.2.5 Swept-Wing Flow 


The subsonic flow past a clean-up glove ol the swept-back tapered h -11 wing ( [16]) was 
selected as a test case. The flow conditions are the same as that used as case 5 in Ref. [16]: 

Moc =0.80 

A, = 22°, A 2 = 4.68° 
a = 0.53° 

Poo = 785.3 lb/ ft 2 
Too = 430" R 
T w = T aw 
c = 5.4256 ft 

which makes the chord Reynolds number 15 Million. 

Figure 11 shows the pressure distribution (measured) on this glove. The skin friction 
coefficient (C/y) lor zero suction as calculated using the present method and the Kaups- 
Cebeci code an' compared in Fig. 12 and are in good agreement. The velocity profile (t>) 
and its first and second normal derivatives at X/c = 0.405 ( 40-th station ) calculated using 
two codes are shown in Fig. 13. l'hc agreement is very good with the exception of small 
differences in the second normal derivative near the wall. 

The logarithmic disturbance amplication factor N was calculated with the envelope 
method using the (.' 0 S A F [18] linear stability analysis code to assess the effect of the differ- 
ent boundarv-layer codes. The envelope method calculates the wave length and orientation 
angle combination that maximizes the disturbance growth for a. fixed frequency. The cur- 
vature effect and non-parellel effect was not considered in this code. The N factor for a 
disturbance frequency of 3000 Ilz was obtained using the result of the present method and 
of the Kaups-Cebeci code and arc' compared in Fig. 14. The difference is within 2 percent. 
For this flow, the most amplified disturbaces are cross-flow vortices (the wave angle and the 
inviscid streamline direction is about S5 degrees) up to the 60 percent chord location. 

The skin friction coefficient obtained using Kaups-Cebeci code for this test case does not 
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show wiggles when suction is not applied. The skin friction coefficient obtained using Kaups- 
Cebeci code shows wiggles for some flow cases. The wiggles are propound especially when the 
pressure distribution near the stagnation point is not smooth or when suction is applied. The 
skin friction coefficient distributions with suction C q = -0.0007 for X/c = 0.051 to 0.332 
obtained using the two boundary-layer codes are compared in the Fig. 15. Immediately 
after the suction region, i.e., from X/c = 0.33, the skin friction coefficient obtained using 
the Kaups-Cebeci code shows an oscillation. The velocity and its normal derivatives at 
X/c = 0.405 (40-th step) are compared on Fig. 16. The first and second normal derivatives 
obtained using Kaups-Cebeci code are not smooth compared with those obtained using the 
present method. 
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6. CONCLUDING REMARKS 


A computer program has been presented for solving the compressible, laminar boundary- 
layer equations for two-dimensional, axisymmetric, and quasi-three-dimensional flows which 
includes flow along the plane ol symmetry, flow along the leading-edge attachment line, 
and swept wing flows with appropriate conical flow approximations. The method has been 
applied to a number of two-dimensional, axisymmetric, and quasi-three-dimensional flow 
cases over a range of speeds and validated. The results indicate that this boundary-layer 
code gives velocity and temperature profiles which are accurate, smooth, and continuous 
through the first and second normal derivatives. This code can be coupled with a stability 
analysis code and used to predict the onset of the boundary-layer transition which enables 
the assessment of the laminar flow control techniques. The CPU time for a typical swept 
wing with 48 streamwise station and 50 to SS grid points across the boundary layer was 
about 4 seconds on the CRAY- 2 or 1.5 seconds on the CRAY Y-MP computer. 
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Appendix A. Block Tridiagonal Matrix Algorithm 


The two vector equations to be solved are 

h = h k . , + A|zl (// , + Hk-i) ( Al.a) 

-A k H k - X + Bklh - C k H k+1 + a k h k = D k (Al.b) 

where A k , Cjt, and a k are 2x2 matrices, Ilk, hk, and Dk are vectors. These equations 
are solved using the Davis Modified Tridiagonal Algorithm. The variables E k , e k and dk are 
introduced such that 

H k = E k H k - 1 +e k h k - l +d k (A2) 

where Ek and e k are 2x2 matrices, and dk is a vector. 

Using Eq. (A2), Eq. (Al.b) becomes 

- AkHk-i + Bklh - CkEk+xlh - C k e k ^hk - C k d k + 1 + a k hk = D k (A3) 

Define 

Rk = a k — CkCk+i (A4) 

Then, Eq. (A3) may be written as 

““ Akll k— 1 + (5fc — Cjb£jf+1 )//fc + Rkhk — Dk + C/jrf/t+i (A5) 

Substituting Eq. (Al.a.) into Eq. (A5) gives 

(~A k + ^±±R k )H k - X + {B k - C k E k+ i + ^t±R k )H k + /?,/>,_! - £», - CY-4+i = 0 (A6) 
Next, define 

Pk = B k - C k E k+ i + ^^-Rk (A7) 

Solving Eq. (A6) for 77^, 

H k = p.T^/U- - ^ i 7?,)77,_ 1 - Pk l Rkh k -i + P k \D k + C k d k+1 ) (AS) 
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Equating Eqs. (A2) and (AS) term by term yields 


c k = ~ Pk R k 


r* — i 4 i ^ Ck— l 

Ek — Vk Ak -\ e k 

<lk — p k l {E) k + C k d k + i ) 


The boundary condition at the edge of the boundary-layer (k 

H kmai' ~ 

This provides the conditions 


is 


1.0 

^■e/krc/ Ol Vyei^ref 


dfctniijr 


1.0 

Ve/VrcJ Ol* l 

V /V ref 


(A 9. a) 
(AD.b) 
( A9.c) 


Climax = Ekmnr = 0 ( A 1 0) 

The parameters of Eq. (A9) are first determined for decreasing values of k (kmax-1, kmax- 
2 ,... ,2) beginning at the edge of the boundary-layer. Then Eqs. (Al.b) and (A2) are solved 
for increasing values of k (k=2,3,..., kmax) using the boundary conditions at the wall, 


III = hi = 0 


(All) 


Appendix B. Boundary-Layer Edge Conditions 


B.l Two-dimensional and Axisvmmetric Flow 

For an isentropic flow (when there is no shock between the free stream and the body), 
the boundary-layer edge velocity is obtained from C v by the following equation. 



However, for a non- isentropic flow (when there is a shock between the free stream and the 
body), the edge velocity is obtained from the C p distribution by integration of the following 
surface Euler equation. 

(l{u e /V<x>) m pood Op /DQ\ 

7 = 7 — ( b 2 ) 

a x p e a x 

which is another form of Eq. (21). This equation is integrated using the fourth-order Runge- 
Kutta method. For a blunted nose body at a supersonic free stream condition, we assume 
that there is a stagnation point where u e is zero and the stagnation pressure was calculated 
assuming that the streamline to the stagnat ion point passed through the normal shock. For 
a sharp nose body at a supersonic free stream condition, the initial edge velocity (n e at i — 1) 
was obtained from the Eq. (Bl). This approximation is based on the assumption that the 
shock angle relative to the free stream is small. It would be better to replace this value if 
an accurate value (u t at i = 1) is available and especially when shock angle is not small. To 
replace this value, put the value in the subroutine EDGECON. 

In the code, if the free-stream Mach number is greater than RISENTD (0.7 is set in 
the code), the flow is assumed to be non-isentropic; and Eq. (B2) is integrated to obtain 
the boundary-layer edge velocity. However, if the velocity (w e ) is available along with the 
pressure distribution from the inviscid code, integration of Eq. (B2) is not needed. In that 
case, use option KUE=1. 

B.2 Swept-wing Flow 
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The method to calculate the coordinate y is the same as used in Ref. [5], The nondimen- 
sional wing thickness distribution zjc is assumed specified as a function of Xjc on the airfoil 
surface. The boundary-layer calculations are done along the arc formed by the intersection 
of a sphere of radius x = Xq and the conical wing surface. This chord intersects the sphere 
at the wing leading-edge as shown in Figs. 3, 4, and 5. A new variable cj> is defined by (see 
Fig. 6) 

Xjc = 1 — coscj) (B3) 

A negative value of Xjc implies negative value of <f>, and <f> — 0 designates the leading-edge. 
The variable y is obtained from the integration of the following equations. 


dy 


d<j> 


dx 


dy, 


dz . 




(16 ' 


d6 ‘ 


For a swept-back wing with the spamvise conical flow assumption (Fig. 3), 


xq ■ 

cos\\(tan\ 

dlv 

-i 

d6 

cosXi Z? 1 / 2 

dy 

Xq dll 

d6 

2/i 3 / 2 (16 

dz 

c 

dcj> 

“ cosXih 1 /* < 


csincj) + 


[.rosmAi — c(l — coscj > )] dh 

2 h dj> 


) 


— 1 + / . — )o + { — 

(X(J COS At Y C { X ( 


— COS6 ) 


.To cos Aj J 

For a swept-forward wing with the spamvise conical flow assumption (Fig. 4), 


XQ - 

CO, 

dx _ 

sX-zitanX 

-1 

dcj) cc 

>sX 2 h'/ 2 

dy _ 

Xq (lh 

d 6 

2/) 3 /' 2 (16 

dz _ 

C 

(16 ci 

AsA 2 //>/ 2 ‘ 

h — 1 -f 

c 2 

{ xq CO. sA 


c sincp + 


[.r 0 sm( — A 2 ) — c(l — coscj))} dh 
2 h ~M> 


d(z/c) 0 

d6 


_ {zfc ) o dh\ 
2 h (16) 


c( 1 — cos6) 
Xq COs\ 2 


(B4) 


(B5.a) 

(B5.b) 

(B5.c) 

(B5.d) 

(B5.e) 


(B6.a) 

(B6.b) 

(B6.c) 

(B6.d) 

(B6.e) 
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For a swept-back wing with the streamwise conical flow assumption (Fig. 5), 


Xq : 

cos A] 

dx 

-1 

d <fi 

s/nAj Z; 1 / 2 

dy 

x 0 dh 

d(j> 

2/j 3 / 2 d<t> 

dz 

c 

d<j) 

" sinXi h 1 / 2 


csin<p + 


[a’oco^Aj — c(l — cos(f > )] dh 
2 h d<f> 


d(z/c ) o _ (z/c ) o dh 
d(j> 2 h dcj) 


h = 1 + 


r(”)o + \ c °t^i “ 


c( 1 — COS<f)) 

xq sinX\ 


(B7.a) 

(B7.b) 

(B7.c) 

(B7.d) 

(B7.e) 


(.r 0 sinXi) 2 c‘ 

(z/c ) o denotes the ordinate of the defining airfoil and is obtained by cubic polynomial in- 
terpolation. Integration of the equations above starts from the stagnation line (leading-edge 
attchment line) where y — 0. 

We solve the following surface Euler equations which are derived from Eq.(33). 

diUe/V^) 


dy 

d(vJV^) 2 


v e /\ 


/ 

CO 


= -2 u e v e /V£ 


p. x d ( 'p 


(BS.a) 

(BS.b) 


dy ^ p e dy 

The equations above are integrated by the fourth order Runge-Kutta method with the initial 
conditions: 


UfjVi jo = ± \J 1 — C pa , v e = 0 for isentropic flow (B9.a) 

WeAoo = i sinXi , v e = 0 for non-isentropic flow (B9.b) 

where C ps is C p at the stagnation line. Because the coordinate systems for swept wing flow 

are defined as in Figs. 3, 4, and 5, u e is negative for swept-back wing flow with the spanwise 

conical flow assumption (Fig. 3) and u e is positive for the other swept wing flow cases(Figs. 4 

and 5). When there is a shock in front of the wing, it is better to replace the initial condition 

(ue/V x — ±sin\i) by a more accurate value if it is available because this approximation is 

based on the assumption that shock is parallel to the stagnation line. To replace the initial 

condition, put the value in the subroutine WING. 
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It is not necessary to integrate Equation (RS.b) for isentropic flow since v e can be obtained 
from known C p and u e . In the code, if free stream Mach number is less than RISEN WG (0.7 
is set in the code), the flow is assumed to be isentropic. 

B.3 Temperature 

The temperature at the edge of the houndary-layer can be obtained using the total 
in viscid velocity, i.e., 

T./T„ = 1 + - (=t-) ! ] (BIO) 

Equation (BIO) is derived from the inviscid energy equation and is valid for all speed 
regimes. Temperature at the boundary- layer edge is calculated using the above equation 
in the boundary-layer code. 
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PART II. USER’S MANUAL 


1. PROGRAM DESCRIPTION 

BLSTA was developed to solve the compressible, laminar boundary-layer equations for the 
two-dimensional flow (KASE=0), axisvmmetric flow (I\ASE=1), and quasi-three-dimensional 
flows such as the flow along the planes of symmetry (I\ASE=2), the flow along the leading- 
edge attachment line (KASE—3), swept-back and swept-forward wing flows with the span- 
wise conical flow assumption (I\ASE=4), and swept-back delta wing flow with the streamwise 
conical flow assumption (KASE=5). Using BLSTA, the flows for wide speed range (from sub- 
sonic to hypersonic) can be solved; the post-shock conditions are not needed as a free-stream 
condition when there is a shock in front of the subject for supersonic flow. The velocity 
profiles and their normal derivatives are calculated to be used by the stability analysis codes 
as follows: e Malik [17] code for KASE=0,1,2,3, and COSAL [18] or COSCUR [19] code for 
KASE=4,5. 

The governing boundary-layer equations are in dimensional form; consequently, all in- 
puts to BLSTA must be consistently dimensional. This code operates in either English 
units(/G lb, sec , °R) or SI(MKS) units(m, kg, sec , °I\ ). However, for swept wing flow ( 
KASE= 4 and 5), use English units to be consistent with the stability analysis codes, COSAL 
or COSCUR, as these stability analysis codes (as written) use English units. 

The code block, COMBLCK, which lists the common blocks, is designed for flexibility 
in changing the dimensions in the spatial <;oordinat,es(;r and c ) and to avoid listing the 
common blocks in each subroutine. This COMBLCK is included in the main program 
and most of the subroutines (except subroutine SY and interpolation subroutines) by an 
‘INCLUDE’ statement. The dimensions ol the common blocks are controlled by changing 
the parameters IMANF and KMAXF in the COMBLCK. If either IMAX or KM AX, which 
are the actual number ol grid points in the x and ^-direction, respectively, happens to 



exeed IMAXF or KMAXF, the boundary-layer calculation should stop and the corresponding 
parameters(IMAXF or KMAXF) should be increased so that they are greater than IMAX 
or KMAX. Also, it should be noted that KMAX, the number of grid points in the normal 
direction, may be increased automatically in the boundary-layer code initially (at i=l ) and 
as marching downstream because of increasing ( e . 

The flow chart of the BLSTA is shown on Fig. 17. The convergence of the boundary-layer 
solution is checked at each i-th step in the subroutine SOLVE. If the convergence fails, the 
flow is most probably separated. The skin friction coefficients (C/ x or C/ y ) are a measure of 
the boundary-layer separation and it decreases rapidly near the separation. 


2. SUBROUTINE DESCRIPTION 


BLPARA 

BLPAR3 

BLPAR4 

COEF 

COEF 1 

COEFLE 

COEFLE 1 

COWING 

COWING 1 

CUBINT 

DERIVY 

EDGECON 

INPUT 

INTEG 

LININT 

NTRID 

OUTCOS 


OUTEMAL 

OUTPUT 


calculates the boundary-layer parameters for I\ASE=0,1,2. 
calculates the boundary-layer parameters for I\ASE=3. 
calculates the boundary-layer parameters for I\ASE=4,5. 
calculates the coefficients t?7j through 77 ?. i 3 lor I\ASE=0,1,2, i>2 
calculates the coefficients m\ through ?// 13 lor I\ASE=0,1,2, i=l 
calculates the coefficients ni\ through 17/13 lor KASE=3, i>2 
calculates the coefficients m j through 77/13 for KASE=3, i=l 
calculates the coefficients mi through mi 3 for KASE=4,5, i>2 
calculates the coefficients ? 7 ?i through 77/13 for I\ASE=4,5, i=l 
cubic polynomial interpolation subroutine 

calculates first and second normal derivatives to be used in the stability 
analysis codes 

For a given C p distribution, locates the stagnation point and calculates the 
boundary-layer edge velocity^ » e ) for KASE=0,1 
reads in input file 
integration subroutine 

linear interpolation subroutine (used if INTEV ^ 0 ) 

solves block tridiagonal matrix equations using Davis Modified Tridiagonal 
Algorithm 

writes the velocity profiles and its normal derivatives on COSLNPUT which 
is to be read by the GOSAL or COSCUR stability analysis code when 
KASE=4,5. 

writes the velocity profiles and its normal derivatives on BFLOW which 
is to be read by the c MaUk stability analysis code when KASE=0,1,2,3. 
writes the input echos, velocity profiles, boundary-layer edge conditions, 
and boundary-layer parameters on FORT. 30 
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SOLVE 


SPLINE 

STAGPT 

SY 

WING 


calculates the coefficients for the tridiagonal systems of equations and 
solves the momentum equations and energy equation iteratively, 
derivative subroutine 

solves the stagnation point flow (key parameter is C*). 

solves the tridiagonal systems of equations using Thomas Algorithm 

calculates velocity components for a given C v distribution (for KASE=4,5) 
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3. PARAMETER AND VARIABLE DIRECTORY 


A(I) 

Xjc for defining airfoil, where c=ehord. (for KASE=4,5) 

AINF 

speed of sound 

BC(I,K) 

C{ = Pn/pePe) at (Xi, z k ) 

BLTH(I) 

8, boundary-layer thickness where V/V e =BLTI\ 

BLTK 

V/V e value for defining boundary-layer thickness, 0.995 is 
the code. 

CAVD(I,J) 

V e , boundary-layer edge total velocity 

CFX(I) 

Cf x 

CFY(I) 

Cf y (for KASE=4,5) 

COSTII(I,J) 

cos 0 

CP 

c p , specific heat 

CPD(I,J) 

Cp, pressure coefficient 

CQD(I) 

Cq , suction coefficient 

CSTAR 

C-(= B/A) 

DK1DY 

d(I\ \ cos ())/ dy 

DSPTX(I) 

<$*, displacement thickness in .r-direction 

DSPTY(I) 

8*, displacement thickness in //-direction (for EASE— 4,5) 

DSTZ 

length scale (/) = y/v t x/u r (KASE=0,1,2,3), = 8* y (EASE 

DX 

A.r ; , = .r, - 

DY(J) 

A//j 

DZETA(K) 

= Cfc+i — G 

DZETA1 

ACi 

GAMMA 

7 

H(1,I,K) 

F at the point (.*,■, 

H(2,I,I\) 

G at, the point. (.)•;, r^.) 

H(3,I,K) 

F at the point (,r M c*.) 
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HN(1,I,K) 

HN(2,I,K) 

HS(1,I,I<) 

HS(2,I,I<) 

H 1 (I, J) 

H2(I,J) 

I 

IL 

IMAX 

IMAXF 

J 

K 

KASE 


IvMAX 


KMAXF 


KROW 


the solution for F from subroutine NTRID 
the solution for G from subroutine NTRID 
f at the point (x;, ~ k ) 
g at the point (x,, z k ) 
hi at the point (x ( , yj ) 
h -2 at the point ( x , , yj) 

index for the boundary-layer grid in the x direction 

i- th step where the boundary-layer calculation stops (not necessarily 

the same as IMAX or ISTOP) 

number of boundary-layer grid points in the x-direction 
maximum possible number of grid points in the x-direction, given 
in COMBLCK 

index for the boundary-layer grid in the y direction 

index for the boundary-layer grid in the ( direction 

index for case of the flow 

=0 for two-dimensional flow 

= 1 for axisymmetric flow 

=2 for the flow along the plane of symmetry 

=3 for the flow along the leading-edge attachment line 

=4 for the swept wing flow with the spamvise conical flow assumption 

=5 for the delta wing flow with the streamwise conical flow assumption 

number of grid points in the c-direction (may be changed as i increases) 

(KM AX < KMAXF) 

maximum possible number of grid points in the s-direction, given in 
COMBLCK 

= 1 when wall mass injection exists 
=0 when there is no wall mass injection 
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KSHALE 

= 1 when the shape of nose is sharp (when there is no stagnation point) 
—0 when the shape of nose is blunted (when there is a stagnation point) 

KTW 

= 1 when the wall temperature is given as a boundary condition 
=0 for adiabatic wall condition 

KUE 

= 1 when the inviscid velocity is given as input (for KASE=0,1) 
=0 when the inviscid velocity is not given as input 

MKS 

=0 when using the English units (ft, lb, sec, °R) 
= 1 when using SI units (m, kg, sec, °K) 


?7? 1 , .., 777-13 

NI 

number of input stations for the streamwise airfoil (for I\ASE=4,5) 

PE(I) 

p, pressure 

PI 

7 T 

PINF 

Pco, free stream pressure 

PR 

Pr, Prandtl Number 

QW(i) 

Q w 

RISENTD 

free-streani Mach number below which isentropic condition is used to calculate 
the edge velocity for 2-D and Axisymmetric flow (0.7 is set in the code) 

RISENWG 

free-stream Mach number below which isentropic condition is used to calculate 
the edge velocity for swept-wing flow (0.7 is set in the code) 

RI< 

k{= AOb+i/ACfc) 

RKAPA(I) 

Kappa value 

RMINF 

A/co 

RMYUED(I) 

Pe 

RNUINF 

V co 

ROED(I) 

Pe 

ROERO(I,I<) 

Pel P 

ROINF 

P CO 
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ROVVVV(I) 

RR 

SWLE 

SWTE 

S1(I,J) 

TD(I,I<) 

TE(I) 

THMOX(I) 

THMOY(I) 

TINF 

TWALL(I) 

TVVTADB(I) 

UE(I,J) 

UKMAXl 


VE(I,J) 

VINF 

VMAX(I) 

XC(I) 

XD(I) 

XI\I(I) 

XPD(I,.J) 

XO 

Y(I) 

YC(I) 


(pU’)w 

gas constant 

leading-edge sweep angle in degrees 
trai ling-edge sweep angle in degrees 

temperature inside the boundary-layer at (a*,-, z^) 
boundary-layer edge temperature 

momentum thickness in ^-direction, defined as — — (1 — —)dz 

momentum thickness in y-direction, defined as ^-^-(1 — — )dz 

for KASE=-1,5 

T w 

-f U: l f a w 
n, or u e j K x , 

KM AX and ( e (=(,"(KMAX)) are to be increased as going downstream 
so that (u / u e )k=h'M ax is greater than this value, typically 0.99999 
v ( or r t / 1 x 
Ko 

maximum crosswise velocity 

X/c (KASE=0,1,4,5) or x (KASE=2,3) 
x for KASE=0, 1,2,3, y for KASE=4,5 

cross-flow Reynolds number, defined as p e c mux &o.oi/f l e 
x' 

x (for KASE=4,5) 

z/c for defining airfoil (for KASE=4,5 only) 

y'/c for 2-1) flow, where c— chord 

(?■ — rif_)/c for axisymmetrie flow, where c= chord 
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y 


YD(J) 

YPD(I,J) y' for 2-D flow, r for axisymmotric flow 

ZACT(I,I<) s 

ZETA(K) C 

ZETAE Ce 



4. INPUT DESCRIPTION 


Inputs to BLSTA are read through subroutine INPUT. There are three input formats 
which differ depending on type of the flow. The necessary input parameters are described 
below. For supersonic flow, and when a shock is present between the free stream and the 
body, we do not have to replace the free stream condition by the post-shock condition, i.e., 
the undisturbed free stream condition is still used for the free stream condition. 

4.1 Two-dimensional and Axisym metric Flow 

KASE =0 for two-dimensional flow 

= 1 for axisymmetric flow 


ISTOP 

MI\S 


KSHALE 


KTW 


KROW 


KUE 


i-th step to stop the computaion 
— 1 when using the SI (MKS) Units (m, kg, sec, 0 I\) 

= 0 when using the English Units ( ft,lb,sec,° R) 

= 1 when the shape of the nose or leading-edge is sharp (if there is no 
stagnation point) 

=0 when the shape of the nose or leading-edge is blunted (if there is a 
stagnation point) 

= 1 when the wall temperature is given as a boundary condition 

(TWTADB(I) must be specified.) 

=0 when adiabatic wall condition is used (TWTADB(I) is not 

needed. Although it might be specified in the input, it will be neglected.) 

= 1 when wall mass injection exists (CQD(I) must be specified) 

=0 when there is no wall mass injection (CQD(I) is not needed. 

Although CQD(I) might be specified below, it will be neglected.) 

= 1 when the inviscid velocity is given as input 

(recommended when there is a shock in front of the body) 
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=0 when the inviscid velocity is not given as input 


PR 

RMINF 

PINF 

TINF 

UK MAXI 

ZETAE 

DZETA1 

RK 

IMAX 

CHORD 

RLE 


INTEV 

ISTABW 


XC(I) 


Prandtl Number 

(free stream Mach number) 

Pco (free stream pressure) in N/m 2 (if MI\S=1) or in lb/ ft 2 (if MKS=0) 
Ttx, (free stream temperature) in °I\ (if MKS=1) or in °R if (MI\S=0) 

ZETAE and KMAX will be increased in BLSTA so that (u/u e ) kmax-i is 
greater than this value, typically 0.99999. 

Ce at i=l (guess value, note that ZETAE and KMAX is increased as 
required even at i=l, and as i increases) 

AC(1) 

AC(fc + l)/AC(*) 


number of input stations where C v is specified 
chord in m(if MKS=1) or in //(if MI\S=0) 

0 for 2-D flow 

i'if , radius of leading-edge, measured from axes of rotation for axisymmetric 
body, in m(if MKS=1) or in //(if MI\S=0) 

number of additional points for the boundary-layer grid between each 
input stations 

=0 output stations for the stability analysis will be the same as the inviscid 
grid points(input stations). 

= 1 output stations for the stability analysis will be the same as the 
boundary-layer grid points (used when INTEV / 0) 

X/c, where c= chord 
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YC(I) 


?//c for 2-13 flow, where c= chord 


(r — ri e )/c for axisymmetric. flow, where c= chord 

cpd ( u ) c p (=(p-p^)/ l i p^yi) 

CQD(I) C„ (= {fnr) w /p. x Vco), C q < 0 for suction 

TWTADB(I) T w /T au , ('This is a required input if KTW=1) 

UE(I,1) u ( / l ’ x , (This is required input if KUE=1) 

For a blunt leading-edge body (I\SIIALE=0), the first input station does not have to 
be the stagnation point; however, the input stations must enclose the stagnation point. 
The code locates the stagnation point and starts the boundary-layer calculations from the 
stagnation point ( The code resets i=l at the stagnation point). For an axisymmetric flow, 
the first input station should be slightly downstream ot the nose tip to avoid /?2 = 0. For an 
axisymmetric nacelle application, use nonzero r/ f which is the radius of the nacelle leading- 
edge measured from the' center ol rotation ol the body. 

In the code, TWALL(I). which will be used as the wall boundary condition, is calculated 
from the input TWTADB(I) (= T w /T aw ) value and the adibatic wall temperature [T aw ) 
which is approximated bv the following equation. 

+ (II.l) 

This equation gives a fairly accutare value for T aw . However, to apply a more precise wall 
temperature condition, specify values directly to TWALL(I). 
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4.2 Flow Along Plane of Symmetry and Along Leading-Edge Attachment Line 


KASE 


ISTOP 

MKS 


I\ SHALE 


KTW 


KROW 


PR 

RMINF 

PINF 

TINF 


UK MAXI 


ZETAE 


=2 lor the flow along t he plane of symmetry 
=3 for the flow along leading-edge attachment line 

i-t.h step to stop the computaion 
= 1 when using the SI (MKS) Units (m, kg, sec,° I\) 

=0 when using the English Units (ft, lb, sec, ° R) 

= 1 when the shape of nose or leading-edge is sharp (if there is no 
stagnation point) 

=0 when the shape of nose or leading-edge is blunted (if there is a 
stagnation point ) 

= 1 when the wall temperature is given as a boundary condition 
(TWTADB(I) must be specified.) 

=0 when adiabatic wall condition is used (TVVTADB(I) is not 

needed. Although it might be specified in the input, it will be neglected.) 
= 1 when wall mass injection exists (CQD(I) must be specified) 

=0 when there is no wall mass injection (CQD(I) is not needed. 

Although CQD(l) might be specified below, it will be neglected.) 

Prandtl Number 

il/cc (free stream Mach number) 

Pot, (free stream pressure) in N/m 2 (if MKS=1) or in lb/ ft 2 (if MKS=0) 

T. Jx, (free stream temperature) in °I\ (if MI\S=1) or in °R if (MI\S=0) 

ZETAE and KM AX will be increased in BLSTA so that (u/u e )/CM.4A'-i is 
greater than this value, typically 0.99999. 

Q at i=l (guess value, note that ZETAE and KMAX is increased as 
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DZETA1 


required even at i= l , and as i increases) 

AC(1) 

AC(*+1)/AC(*) 


RK 

IMAX number of input stations where C p is specified 

XD(I) 

YD(J) 

XPD(I,J) 

YPD(I,J) 

S1(I,J) 

UE(I,J) 

VE(I,J) 

H1(I,J) 

H2(I,J) 

COSTH(I,J) 

CPD(I,1) C P (=(p-p»)/J/>»V«) 

CQD(I) C q (— (pw)w! P ooVqq)) C q < 0 for suction 

TWTADB(I) T w / T aw (This is a required input if KTW— 1) 

The location of the first input station should be slightly downstream of the stagnation 
point or nose tip for the flow along the plane of symmetry (I\ASE=2). j = l designates the 
plane ot symmetry or the leading-edge attachment line. j=2 is off the plane of symmetry or 
off the leading-edge attachment line. 

To apply a more precise wall temperature condition, specify values directly to TWALL(I) 
because T aw was approximated by the Eq. (II. 1). 


x 

y 

x f (for reference) 
y f (for reference) 

Ue/Vn 

Ve/Vo* 

hi 
h 2 

cosO 
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4.3 Swept Wing Flow 


Input parameters written in the input file for swept wing flows are the same as Kaups- 
Cebeci code [5] input except the first two lines. 

KASE =4 for swept wing flow with the spanwise conical flow assumption 

=5 for delta wing flow with the streamwise conical flow assumption 


KTW 


KROW 


INTEV 


ISTABW 


= 1 when the wall temperature is given as a boundary condition 
(TWTADB(I) must be specified.) 

=0 when adiabatic wall condition is used (TWTADB(I) is not 

needed. Although it might be specified in the input, it will be neglected.) 
= 1 when wall mass injection exists (CQD(l) must be specified below) 

=0 when there is no wall mass injection (CQD(I) is not needed. 

Although CQD(I) might be specified below, it will be neglected.) 
number of additional points for the boundary-layer grid between each 
inviscid input station. If you want 2 more boundary-layer grid points 
between each inviscid input station, set INTEV=2. This is useful if the 
inviscid grid is sparce or highly nonuniform. If you want the boundary-layer 
grid the same as the inviscid grid, set INTEV=0. Setting INTEV > 0 may 
cause early convergence failure if the pressure distribution is not reasonable 
near the stagnation point. 

-0 output stations for the stability analysis will be the same as the inviscid 
grid points(input stations). 

= 1 output stations for the stability analysis will be the same as the 
boundary-layer grid points (used when INTEV / 0) 


NI 


number of input stat ions tor the airfoil 

number of input stations where pressure is specified 
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IMAX 



ZETAE 

DZETA1 

RK 

CHORD 

SVVLE 

SWTE 

RMINF 
U REE 
PINF 
TINF 
PR 

A(I) 

Y ( I ) 

XC( I) 

CPI)(I,1) 

CQD(I) 

TYVTADB(I) 

RKAPA(I) 


Est imat'd value of (,' e at the first station. 

AC(1) 

AC(k+l)/Aak) 


chord length in ft for the airfoil 

leading-edge sweep angle in degrees (negative for swept-forward wing) 

< railing-edge sweep angle in degrees (negative for swept-forward wing) 

(not used for KASE=5) 
d /, x , (free stream Mach number) 

l x (free* stream velocity) in J 't/sec, input only if = 0. 

7 J X (free' stream pressure) in lb/ ft 2 
/’ x , (free stream temperature) in °R 
Prandtl Number 

X/c for defining airfoil, where c= chord. A(I)=0.0 must be input if 
calculations contain the leading edge. Total of NI points 
c/e for defining airfoil, where c= chord. Y(I)=0.0 must be input if 
calculations contain the leading edge, dotal of NI points 

X/c of input stations, total of IMAX points, must include the stagnation 
line at 7 = 1. 

C v values at X/c of input stations, total of IMAX points 
( values (= {pw) w / P<x>V<x,)i C’ y < 0 for suction. Total of IMAX points. 

7,, / values, total of IMAX points. Required if I\TW=1 

Kapj^a value', total of IMAX points. Required for COSCIJR stability code. 
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The following concerning sweep angles Aj and \ 2 must be observed. 


Aj >0, Ai > X 2 for swept-back wing flow (KASE=4) (II. 2. a) 

A! < 0, A 2 < 0, Aj > A 2 lor swept-forward wing flow (I\ASE=4) (II. 2. b) 
Aj > 0 for delta wing flow (KASE=5) (II. 2. c) 


As described in PART I, Appendix B, this code does not assume isentropic. flow to calcu- 
late the velocity components when the free stream Mach number is greater than RISEN WG; 
therefore, supersonic flow can be solved without replacing the free stream condition by the 
condition downstream ol the shock. The sweep angle docs not need to be changed either. 
Airfoil coordinates and the pressure input, for KASE=4 are along the streamwisc direction 


as shown on Figs. 3 and 1. However, airfoil coordinates and the pressure input for KASE — 5 
should be along the direction as indicated on Fig. 5. A few extra points ahead of the 
stagnation point should be included in the definition of the airfoil for proper interpolation 
of the geometrical data near the stagnation point. 

The first, input station must be the stagnation point and the pressure coefficient at this 
point can be obtained analytically by the following equation for isentropic flow (subsonic 
flow). 

c p , = {(1 + Micas 2 - 1 )l\Ml (II 3) 

For the swept wing flows, a variable normal grid is necessary because Q increases rapidly 


as i increases. The present code employs a stretched normal grid which applies a. constant 
ratio between two adjacent grid spacings, i.e., A((k + ])/A((k) = k — const. I he total 
number of points across the boundary layer, KMAX, can be calculated by the following 
equation: 


KMAX = 


in\\ + (k. — 1 )(Ce/ACi )] 
In k 


( 1 1 - 4 ) 


A uniform normal grid, of course, can be obtained by setting k = 1. For typical swept 
wing flows, k = 1.05 is a reasonable value. For a quick estimate of input ZETAE (( e ) for a 
typical swept wing, the following guide may be useful although it may differ depending on 


51 



the velocity gradient at the stagnation point. 

A| — /\'2 ~ 0.01 (in degrees) — > ( e < 0.01 (II. 5. a) 

A| - A 2 ~ 1.0 (in degrees) — > Q < 0.1 (II. 5. b) 

A[ - A 2 ss 10. (in degrees) — > Q < 0.2 (II. 5. c) 

Ai — A 2 ~ 20. (in degrees) — > ( £ < 0.4 (II. 5. d) 

Aj — A 2 > 30. (in degrees) — ► C e < 1.0 (II. 5. e) 


Assuming k = 1.05, if we choose A;*! = de/500, we have 67 grid points at the stagnation 
point. 

In the code, TYVALL(I), which will be used as the wall boundary condition, is calculated 
from the input TWTADB(l) (= r I' u ,/T aw ) value and the adibatic wall temperature (T aw ) 
which is ajjproxi mated by the Kq. (II. 1). Equation (II. 1) gives a fairly accutare value for 
T aw . However, to apply a more precise wall temperature condition, specify values directly 
to TYVALL(I). 

BLSTA can use either SI unit or English unit. However, since the stability analysis codes, 
COSAL and COSCUR use English units, English units should be used for this case (swept 
wing flows, I\ASE=4 and 5). 
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5. OUTPUT DESCRIPTION 


5.1 Standard Output. 

The standard output iroin BEST A is written on FORT. 30 by subroutine OUTPUT.. 

(1) The flow condition and other input parameters are echoed. This includes 

KASE, ISTOP, MKS, K SHALE, KTW, KROW, PR, RMINF, PINF, TINF, 
UKMAX1, ZETAE, DZETAl, RK, IMAX, CHORD, RLE, INTEV, and ISTABW. 

NI, SWLE, SWTE is added for KASE=1 and 5. 

Also the calculated free-stream conditions, CP(c ; ,, specific heat), ROINF(/) cc ), 
RMYUINF (/loo), RNUINF(//oo), AINF(« 00 ), \TNF(I4,) are printed. 

(2) The velocity and temperature profiles tor all the i step (i=l,..,IL) are printed: 

c, F(= u/tte), G (= v/Vrcf or vj V ref ), T/T e (= pjp) for k=f ,2,..,KMAX, i=l,2,..,IL 

(3) The boundary-layer edge conditions are printed: 

it e , e e , , 1 e , Pc , Pc , »'c, fei and Afc lot i 1,2,..,IL 

(4) The boundary-layer parameters are printed: 

XC (A'/c), CFX(C;,.), BLTH(^), DSPTX(<5*), THMOX(x-momentnm thickness), 
CFY(C/j,), DSPTY(^*), TITMOY(y-momentum thickness), QVV ! {q w ), T\VALL(T’ W ) 
for i=l,2,..,lL 

The unit for q u , is I V/ni* (il MKS=1) or Btu/sec/ ft 2 (if MKS=0). 

5.2 Output to Stability Analysis Codes 

The profiles of the flows of KASE=0,1,2,3 are written by the subroutine OUTEMAL on 
BFLOW which can be read by the 2-D and Axisymmetric, spatial, linear stability analysis 
code, e Malik [17]. 
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BFLQW 


IMAXW 

IPS 

KMAX 

DSTZ 

REY 

RES 

EPSXR 


DRDX 

YPD(I,1) 

THMOX(I) 

BLTH(I) 

RET 1 1 ETA 

RETIIM 

PR 

MKS 


mmibcr of output stations 
station number 

number of points inside the boundary layer 
length scale / = \J jit, 

Reynolds number, = u e l/v c where / = \J u t x j a. 

Length Reynolds number, = ( u e l/v e ) 2 
nondimensional streamwise curvature, = Ik 
k >0 for convex, /« <0 for concave 

dvfda, where /■ is local radius of body and s measured along the body 
r; local radius of body 
m om e n t. u m t h i c k j less 
boundary- layer thickness 

momentum thickness Reynolds number (Reg) 

Rtfl/M, 

Prandtl Number 

= 1 when using the SI (MKS) Units (in, kg,sec,° I\) 

=0 when using the English Units ( ft,lb,sec, u R ) 


TE(I) 

RMYUED(I) 

UE(I,i) 

S1(I,1) 

XC(I) 


boundary-layer edge temperature (2’ e ) 
boundary-layer edge viscosity (/< c ) 

u, _ 

s 

X/c (added for reference) 


ZACT(I,I\)/DSTZ actual distance normal to the wall devided by length scale, zjl 
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I 


11(1, 1, K) 

U/U c 

UP(K) 

d(u/u t )/d(z/l) 

IJDP(K) 

cP(u/u t )/d(zllf 

ROERO(I,K) 

pjp{ = T/T e ) 

TP(K) 

d(T/'I\)/d{:/l) 

TDP(K) 

d 2 (T/T c )/d(=/l) 


The profiles of the flows for I\ASE=1,5 (swept wings) arc written by subroutine OUTCOS 
on COSLNPUT which can be read by the 3-D, temporal, linear stability analysis codes, 
COSAL [18] and COSCIMl [19]. 

CQSbNPUT 


IMAXVV 

number of out put, stat ions 

XO 

x 

CHORD 

chord length in j f ior t he airfoil 

SWLE 

lea.cling-<'dge swt'cp angle 

VINF 

Vco 

XC(I) 

X/c 

XD(I) 

y 

VE(I,1) 

V c 

RRW(I) 

s : arc length along x coordinate 

IPS 

station number 

KM AX 

number of points inside the boundary 

DSTZ 

lengt h scale / = S* 

REY 

Reynolds number, = v f l/u e , where / = 



XI\I(I) cross-flow Reynolds number, - p e c max 6 0 . oi/p-e, for COSCUR code 

EPS 6* R K A P A ( I ) / ( c s i n ( S W L E ) ) , for COSCUR code 

UTDEL S/c, for COSCUR code 

cavd(i,i) v; 

PE(1) Pe 

TE(I) T f 

ROED(I) Pt 

RMYUED(I) p c 

TWALL(I) T w 

R0ED(I)/R0ER0(I,1) p u . 

ZACT(I,K)/DSTZ zjl 
VTEMP(K) v/v £ 

VP(K) d(v/v e )/d(z/l) 

VDP(K) d 2 (v/v e )/d(z/l) 2 

UTEMP(K) -u/v e 

UP(K) d(-u/v,)fd{z/l) 

UDP(K) <P(-u/v e )/d{z/l ) 2 

ROERO(I,K) p e /p(=T/T e ) 

TP(K) d(T/T e )/d(z/l) 

TDP(K) d 2 (T/T e )/d(z/l) 2 
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6. SAMPLE CASE INPUT 


6.1 Flat plate 


Flat plate 
0 


istop mks kshale 

ktw krow kue 



21 0 

1 

0 0 0 



pr rminf 


pint tint 

ukmaxl zetae 

0.72 2.0 


25334.3 518.4 

0.999999 

8.0 

imax chord 

rle 

intev istabw 



21 1 . 

0 . 

0 0 



x/c 

y/c 

cp 

cq 

twtadb 

0.000000 

0.0 

0.000000 

0.000000 

1.0 

0.010000 

0.0 

0.000000 

0.000000 

1.0 

0.020000 

0.0 

0.000000 

0.000000 

1.0 

0.030000 

0.0 

0.000000 

0.000000 

1.0 

0.040000 

0 . 0 

0.000000 

0.000000 

1.0 

0.050000 

0.0 

0.000000 

0.000000 

1.0 

0.060000 

0.0 

0.000000 

0.000000 

1.0 

0.070000 

0.0 

0.000000 

0.000000 

1.0 

0.080000 

0.0 

0.000000 

0.000000 

1.0 

0.090000 

0.0 

0.000000 

0.000000 

1.0 

0.100000 

0 . 0 

0.000000 

0.000000 

1.0 

0.110000 

0.0 

0.000000 

0.000000 

1.0 

0.120000 

0.0 

0.000000 

0.000000 

1.0 

0 . 130000 

0.0 

0 . 000000 

0.000000 

1.0 

0.140000 

0.0 

0.000000 

0.000000 

1.0 

0.150000 

0.0 

0.000000 

0.000000 

1.0 

0.160000 

0.0 

0.000000 

0.000000 

1.0 

0.170000 

0.0 

0.000000 

0.000000 

1.0 

0.180000 

0.0 

0.000000 

0.000000 

1.0 

0.190000 

0 . 0 

0.000000 

0.000000 

1.0 

0.200000 

0 . 0 

0.000000 

0.000000 

1.0 


dzetal rk 
0.1 1.0 
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6.2 NACA 0012 airfoil 


This input file was used to generate the result presented in PART I, section 5.2.1. 


NACA 0012 airfoil 
0 


istop mks kshale ktw krow kue 
85 0 0 0 0 0 

pr rminf pinf tinf ukmaxl zetae 

dzetal rk 

0.72 0.5 

2116.0 

520.0 0, 

.99999 7.0 

0.100 1 . 

imax chord 
85 0.28266 

x/c 

rle intev 
0 . 0 
y/c 

istabw 

0 

cp 

cq 

twtadb 

0.000000 

0.000000 

1.064070 

0.000000 

1.000000 

0.000775 

0.004841 

0.892081 

0.000000 

1.000000 

0.002955 

0.009464 

0 . 610767 

0.000000 

1.000000 

0.006215 

0.013564 

0.324677 

0.000000 

1.000000 

0.010168 

0.017176 

0 . 104429 

0.000000 

1.000000 

0.014600 

0.020384 

-0 .047884 

0.000000 

1.000000 

0.019392 

0.023265 

-0.152775 

0.000000 

1.000000 

0.024453 

0.025893 

-0.232184 

0.000000 

1.000000 

0.029746 

0.028308 

-0.292468 

0.000000 

1.000000 

0.035235 

0.030530 

-0.333274 

0.000000 

1.000000 

0.040894 

0.032593 

-0.363119 

0.000000 

1 .000000 

0.046701 

0.034525 

-0.389434 

0.000000 

1.000000 

0.052644 

0.036337 

-0 .411557 

0.000000 

1.000000 

0.058707 

0.038036 

-0 .428407 

0.000000 

1.000000 

0.064883 

0.039631 

-0.440971 

0.000000 

1.000000 

0.071161 

0.041134 

-0 .450779 

0.000000 

1.000000 

0.077527 

0.042547 

-0 .457964 

0.000000 

1.000000 

0.083979 

0.043883 

-0.464400 

0.000000 

1.000000 

0.090506 

0.045144 

-0.470200 

0.000000 

1.000000 

0.097100 

0.046333 

-0.475096 

0.000000 

1.000000 

0 .103757 

0.047452 

-0.477730 

0.000000 

1.000000 

0.117234 

0.049494 

-0 .479044 

0.000000 

1.000000 

0.130892 

0.051303 

-0 .479344 

0.000000 

1.000000 

0 .144698 

0.052898 

-0.477948 

0.000000 

1.000000 

0.158623 

0.054293 

-0 .474022 

0.000000 

1.000000 

0.172644 

0.055505 

-0 .469354 

0.000000 

1.000000 

0 .186741 

0.056548 

-0.463392 

0.000000 

1 .000000 

0.200899 

0.057432 

-0.456391 

0.000000 

1.000000 

0.215106 

0.058169 

-0 .449168 

0.000000 

1.000000 

0.229351 

0.058767 

-0 .439979 

0.000000 

1.000000 

0.243625 

0.059238 

-0.430423 

0.000000 

1.000000 

0.257924 

0.059598 

-0 .422939 

0.000000 

1.000000 

0.272242 

0.059843 

-0 .413984 

0.000000 

1.000000 

0.286573 

0.059980 

-0 .403774 

0.000000 

1.000000 

0.300918 

0.060020 

-0.394188 

0.000000 

1.000000 

0.315270 

0.059967 

-0.384397 

0.000000 

1.000000 

0.329629 

0.059825 

-0.374046 

0.000000 

1.000000 

0.343992 

0.059600 

-0.363322 

0.000000 

1.000000 

0.358359 

0.059299 

-0.353406 

0.000000 

1.000000 

0.372730 

0.058926 

-0.343330 

0.000000 

1.000000 

0.387100 

0.058483 

-0.333330 

0.000000 

1.000000 

0.401472 

0.057975 

-0.323081 

0.000000 

1.000000 
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0.415846 
0.430217 
0.444588 
0.458957 
0.473326 
0.487695 
0.502059 
0.516423 
0 . 530785 
0.545145 
0.559501 
0.573856 
0.588207 
0 . 602554 
0.616900 
0 . 631242 
0.645577 
0.659910 
0.674235 
0.688554 
0.702865 
0 .717166 
0.731454 
0.745729 
0.759983 
0.774215 
0.788415 
0.802579 
0.816697 
0.830753 
0.844735 
0.858622 
0.872391 
0 .886011 
0.899451 
0.912667 
0.925617 
0.938253 
0 . 950523 
0 . 962379 
0 . 973779 
0.984674 
0 .995024 


0.057406 
0.056774 
0.056083 
0.055336 
0 .054544 
0.053701 
0.052810 
0.051879 
0.050902 
0.049879 
0.048818 
0.047721 
0.046587 
0 .045420 
0 . 044224 
0.042991 
0 .041724 
0 .040428 
0 . 039102 
0 . 037745 
0.036361 
0 .034953 
0.033514 
0.032049 
0 . 030562 
0.029045 
0.027504 
0.025945 
0 . 024364 
0 . 022760 
0.021138 
0.019503 
0.017860 
0.016207 
0.014549 
0.012891 
0.011241 
0 .009608 
0.007999 
0.006421 
0 . 004881 
0.003386 
0.001952 


- 0.313969 
- 0.303939 
- 0.293325 
- 0.282234 
- 0.273298 
- 0.263542 
- 0.252860 
- 0.244331 
- 0.234945 
- 0.224334 
- 0.214411 
- 0.205009 
- 0 .195108 
- 0.185450 
- 0.177266 
- 0.167714 
- 0.157628 
- 0.148539 
- 0.139062 
- 0 . 129007 
- 0 .119126 
- 0.110205 
- 0.099750 
- 0.089578 
- 0.080102 
- 0.068917 
- 0.057411 
- 0.046723 
- 0.035715 
- 0.023096 
- 0.009978 
0.003225 
0.016574 
0.031380 
0.047433 
0.064860 
0.084260 
0 .105485 
0 . 129255 
0 .157685 
0 .192978 
0.242692 
0.343002 


0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 

0.000000 


1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1 .000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 

1.000000 
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6.3 Hypersonic cone 


This input file was used to generate the result presented in PART I, section 5.2.2. 


Hypersonic cone 
1 


istop mks kshale 

ktw krow 

kue 



58 1 1 

1 1 


0 



pr rminf 

pinf tinf 

ukmaxl 

zetae zetal rk 

0.72 7.4 

701.4 69, 

, 7 

0 . 99999 

5.0 0. 

.005 1.05 

imax chord rle 

intev istabw 




58 1. 0. 

0 0 





x/c 

(r-rle) /c 


cp 

cq 

t wtadb 

0.005000 

0.000437 


0.020200 

0.000000 

0 .440700 

0.006000 

0.000525 


0.020200 

0.000000 

0 .440700 

0.007000 

0.000612 


0.020200 

0.000000 

0 . 440700 

0.008000 

0.000700 


0.020200 

0.000000 

0 . 440700 

0.009000 

0.000787 


0.020200 

0.000000 

0 .440700 

0.010000 

0.000875 


0.020200 

0.000000 

0 .440700 

0.012000 

0.001050 


0.020200 

0.000000 

0 .440700 

0.014000 

0.001225 


0.020200 

0.000000 

0 .440700 

0.016000 

0.001400 


0.020200 

0.000000 

0 . 440700 

0.018000 

0.001575 


0.020200 

0.000000 

0 .440700 

0.020000 

0.001750 


0.020200 

0.000000 

0 .440700 

0.022000 

0.001925 


0.020200 

0.000000 

0 . 440700 

0.024000 

0.002100 


0.020200 

0.000000 

0 .440700 

0.026000 

0.002275 


0.020200 

0.000000 

0 .440700 

0.028000 

0.002450 


0.020200 

0.000000 

0 . 440700 

0.030000 

0.002625 


0.020200 

0.000000 

0 .440700 

0.035000 

0.003062 


0.020200 

0.000000 

0 . 440700 

0.040000 

0.003500 


0.020200 

0.000000 

0 .440700 

0.045000 

0.003937 


0.020200 

0.000000 

0 .440700 

0.050000 

0.004374 


0.020200 

0.000000 

0 .440700 

0.055000 

0.004812 


0.020200 

0.000000 

0 .440700 

0.060000 

0.005249 


0.020200 

0.000000 

0 .440700 

0.065000 

0.005687 


0.020200 

0.000000 

0 .440700 

0.070000 

0.006124 


0.020200 

0.000000 

0 . 440700 

0.075000 

0.006562 


0.020200 

0.000000 

0 .440700 

0.080000 

0.006999 


0.020200 

0.000000 

0 .440700 

0.085000 

0.007437 


0.020200 

0.000000 

0 .440700 

0.090000 

0.007874 


0.020200 

0.000000 

0 .440700 

0.095000 

0.008311 


0.020200 

0.000000 

0 .440700 

0.097000 

0.008486 


0.020200 

-0.0020754 

0 .440700 

0.099000 

0.008661 


0.020200 

-0.0020754 

0 .440700 

0.101000 

0.008836 


0.020200 

-0.0020754 

0 .440700 

0.103000 

0.009011 


0.020200 

-0.0020754 

0.440700 

0 . 105000 

0.009186 


0.020200 

-0.0020754 

0 . 440700 

0.107000 

0.009361 


0.020200 

-0.0020754 

0.440700 

0.109000 

0.009536 


0.020200 

-0.0020754 

0.440700 

0 . 111000 

0.009711 


0.020200 

-0.0020754 

0.440700 

0.113000 

0.009886 


0.020200 

-0.0020754 

0.440700 

0.115000 

0.010061 


0.020200 

-0.0020754 

0.440700 

0.125000 

0.010936 


0.020200 

-0.0020754 

0.440700 

0.135000 

0.011811 


0.020200 

-0.0020754 

0.440700 

0.145000 

0.012686 


0.020200 

-0.0020754 

0 .440700 
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i 


0.155000 

0.165000 

0.175000 

0.185000 

0.195000 

0.205000 

0.215000 

0.225000 

0.235000 

0.245000 

0.255000 

0.265000 

0.275000 

0.285000 

0.295000 

0.305000 


0.013561 
0.014436 
0.015311 
0.016185 
0.017060 
0.017935 
0.018810 
0.019685 
0.020560 
0.021435 
0.022310 
0 . 023184 
0 .024059 
0.024934 
0.025809 
0.026684 


0 . 020200 
0 . 020200 
0 . 020200 
0.020200 
0.020200 
0.020200 
0.020200 
0.020200 
0 . 020200 
0 . 020200 
0 . 020200 
0.020200 
0.020200 
0.020200 
0.020200 
0 . 020200 


0.0020754 

0.0020754 

0.0020754 

0.0020754 

0.0020754 

0.0020754 

0.0020754 

0.0020754 

0.0020754 

0.0020754 

0.0020754 

0.0020754 

0.0020754 

0.0020754 

- 0.0020754 

- 0.0020754 


0.440700 
0 .440700 
0.440700 
0 . 440700 
0 . 440700 
0 . 440700 
0 .440700 
0 .440700 
0.440700 
0. 440700 
0 .440700 
0 . 440700 
0 .440700 
0 .440700 
0.440700 
0 .440700 
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6.4 Axisymmetric engine nacelle 


Engine nacelle external flow 
1 

istop mks kshale ktw krow kue 
56 0 0 0 0 0 

pr rminf pinf tinf ukmaxl zetae dzetal rk 

0.72 0.82 498.145 394.013 0.99999 5.0 0.05 1.0 

imax chord rle intev istabw 
56 14.3 4.25 0 0 

x/c (r-rle) / c cp cq twtadb 

0.0238683256 - 0.0204869641 0.6151417248 0.0000000000 1.0 

0.0164744072 - 0.0168793666 0.7887464320 0.0000000000 1.0 

0.0103140282 - 0.0128975952 0.9718194839 0.0000000000 1 0 

0.0055795044 - 0.0086236824 1.1049125354 0.0000000000 l!o 

0.0021887510 - 0.0045658737 1.1124471984 0.0000000000 1.0 

0.0000000000 0.0000000000 0.7048278453 0.0000000000 l!o 

0.0014599675 0.0051387024 - 0.0842995853 0.0000000000 1 .' 0 

0.0056109715 0.0085813903 - 0.5854369694 0.0000000000 1.0 

0.0111811832 0.0113131167 - 0.6227023027 0.0000000000 10 

0.0176661613 0.0141906722 - 0.5492159541 0.0000000000 1 0 

0.0251549880 0.0170536214 - 0.5007490775 0.0000000000 1 0 

0.0335803651 0.0199254323 - 0.4527134120 0.0000000000 1 .' 0 

0.0428933413 0.0227864298 - 0.4200126030 0.0000000000 1 0 

0.0530431706 0.0256447006 - 0.4025171908 0.0000000000 1 .’ 0 

0.0639960788 0.0284675950 - 0.3966051721 0.0000000000 1 0 

0.0757221834 0.0312211987 - 0.3897372348 0.0000000000 1 .0 

0.0881880014 0.0339075743 - 0.3754882614 0.0000000000 1 0 

0.1013622347 0.0365418465 - 0.3572414199 0.0000000000 10 

0.1152177741 0.0391388779 - 0.3416229385 0.0000000000 1.0 

0.1297320899 0.0417051878 - 0.3326261400 0.0000000000 1.0 

0.1448862411 0.0442393140 - 0.3309526130 0.0000000000 1.0 

0.1606642699 0.0467312006 - 0.3353721378 0.0000000000 1 0 

0.1770521761 0.0491657984 - 0.3430597205 0.0000000000 10 

0.1940369963 0.0515277730 - 0.3516655119 0.0000000000 l'o 

0.2116064374 0.0538032291 - 0.3601120215 0.0000000000 1.0 

0.2297489695 0.0559787605 - 0.3682139709 0.0000000000 1.0 

0.2484536975 0.0580415091 - 0.3760610538 0.0000000000 10 

0.2677103674 0.0599788028 - 0.3837500080 0.0000000000 1.0 

0.2875092786 0.0617779173 - 0.3913308636 0.0000000000 1.0 

0.3078412841 0.0634258019 - 0.3988228827 0.0000000000 1 0 

0.3286977320 0.0649087944 - 0.4062423817 0.0000000000 1.0 

0.3500704366 0.0662123468 - 0.4135845297 0.0000000000 1.0 

0.3719516374 0.0673206586 - 0.4207987245 0.0000000000 1.0 

0.3943339817 0.0682164201 - 0.4277360508 0.0000000000 1.0 

0.4172104607 0.0688805972 - 0.4340833245 0.0000000000 l!o 

0.4405743102 0.0692928918 - 0.4394972038 0.0000000000 1.0 

0.4644189467 0.0694321954 - 0.4438936907 0.0000000000 1.0 

0.4887377690 0.0692782623 - 0.4479707927 0.0000000000 1.0 

0.5135242811 0.0688102521 - 0.4536920492 0.0000000000 1.0 

0.5387721440 0.0680053963 - 0.4641817715 0.0000000000 1.0 

0.5644758112 0.0668267161 - 0.4817351890 0.0000000000 1 0 

0.5906321308 0.0651894933 - 0.5025687672 0.0000000000 1 .' 0 

0.6172386384 0.0629741575 - 0.5127644663 0.0000000000 1 0 

0.6442919546 0.0600429897 - 0.4899800994 0.0000000000 1 .' 0 
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0.6717855835 
0.6997076802 
0 . 7280428406 
0.7567745415 
0.7858881302 
0.8153729797 
0.8452233688 
0.8754380855 
0.9060191458 
0.9369702317 
0 . 9682957072 
1.0000000000 


0.0562771302 
0.0516424902 
0 .0461947096 
0.0400791629 
0.0335065360 
0.0267229901 
0.0199787301 
0 .0135004870 
0 .0074787590 
0.0020876729 
- 0 .0024966698 
- 0.0060581764 


- 0 .4192777548 
- 0.3117699386 
- 0 . 1975039140 
- 0 .0960932589 
- 0 . 0073313301 
0 . 0716443622 
0.1392549527 
0 .1947019608 
0.2388075692 
0.2754120815 
0.2995486827 
0.3146850653 


0.0000000000 

0.0000000000 

0.0000000000 

0.0000000000 

0.0000000000 

0.0000000000 

0.0000000000 

0.0000000000 

0.0000000000 

0.0000000000 

0.0000000000 

0.0000000000 


1 . 0 
1 . 0 
1.0 
1 . 0 
1 . 0 
1 . 0 
1 . 0 
1 . 0 
1 . 0 
1.0 
1.0 
1.0 
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6.5 Ellipsoid of revolution along windward plane of symmetry 

This is reduced input file from which was used to generate the result presented in PART I, 
section 5.2.3. 


Ellipsoid of revolution, windward plane of symmetry 
2 

istop mks kshale ktw krow 
20 0 0 00 

pr rminf pinf tinf ukmaxl zetae dzetal rk 

0.72 0.05 2116. 520.0 0.99999 6.0 0.1 1.0 

imax 

20 

0 . 400000E-02 0 . 450000E-02 0.500000E-02 0.550000E-02 0.600000E-02 

0 . 800000E-02 0 . 100000E-01 0.120000E-01 0.140000E-01 0.160000E-01 

0 . 180000E-0 1 0.200000E-01 0.220000E-01 0.240000E-01 0.260000E-01 

0 . 280000E-01 0.300000E-01 0.320000E-01 0.340000E-01 0.360000E-01 

0 . OOOOOOE+OO 0 . 872665E-01 

1 1 0 . 4 00000E-02 0. 000000E+00 0.226936E-01 0.180332E + 00 

0. 000000E+00 0 . 2 96069E+0 1 0.223383E-01 0. 000000E+00 

1 2 0 . 400000E-02 0.194691E-02 0.226936E-01 0.181028E+00 

0.169429E-01 0.296069E+01 0.223383E-01 0. 000000E+00 

2 1 0 . 450000E-02 0. 000000E+00 0.241363E-01 0.201077E+00 

0. 000000E+00 0 . 281027E+01 0.236904E-01 0. 000000E+00 

2 2 0 . 450000E-02 0.206475E-02 0.241363E-01 0.201768E+00 

0. 16942 9E-01 0.281027E+01 0.236904E-01 0. OOOOOOE+OO 

3 1 0 . 500000E-02 0. 000000E+00 0.255099E-01 0.220376E+00 

0. OOOOOOE+OO 0 . 2 6838 7E+0 1 0.249687E-01 0. OOOOOOE + OO 

3 2 0 . 500000E-02 0.217617E-02 0.255099E-01 0.221063E+00 

0 . 169429E-01 0.268387E+01 0.249687E-01 0. OOOOOOE+OO 

4 1 0 . 550000E-02 0. OOOOOOE+OO 0.268248E-01 0.238433E+00 

0. 000000E+00 0.257 58 5E+0 1 0.261841E-01 0. 000000E+00 

4 2 0.550000E-02 0.228210E-02 0.268248E-01 0.239115E+00 

0.169429E-01 0.257585E+01 0.261841E-01 0. 000000E+00 

5 1 0 . 600000E-02 0. 000000E+00 0.280893E-01 0.255407E+00 

0. OOOOOOE+OO 0.248224E+01 0.273450E-01 0. OOOOOOE+OO 

5 2 0 . 600000E-02 0.238328E-02 0.280893E-01 0.256084E+00 

0 . 1 6942 9E-0 1 0.248224E+01 0.273450E-01 0. 000000E+00 

6 1 0 . 800000E-02 0. 000000E+00 0.327760E-01 0.314700E + 00 

0. 000000E+00 0 . 220441E+01 0.315595E-01 0. 000000E + 00 

6 2 0 . 800000E-02 0.275059E-02 0.327760E-01 0.315359E+00 

0 . 169429E-01 0.220441E+01 0.315595E-01 0. OOOOOOE+OO 

7 1 0 . 100000E-01 0. OOOOOOE+OO 0.369999E-01 0.363744E + 00 

0. OOOOOOE+OO 0.20194 6E + 0 1 0.352668E-01 0. OOOOOOE+OO 

7 2 0 . 100000E-01 0.307371E-02 0.369999E-01 0.364387E+00 

0 . 169429E-01 0.201946E+01 0.352668E-01 0. OOOOOOE+OO 

8 1 0 . 120000E-01 0. OOOOOOE+OO 0.409054E-01 0.405468E + 00 

0. OOOOOOE+OO 0 . 188610E+01 0.386135E-01 0. OOOOOOE+OO 

8 2 0.120000E-01 0.336539E-02 0.409054E-01 0.406095E+00 

0.169429E-01 0.188610E+01 0.386135E-01 0. OOOOOOE+OO 

9 1 0.140000E-01 0. OOOOOOE+OO 0.445763E-01 0.441657E+00 

0. OOOOOOE + OO 0.17847 6E+0 1 0.416863E-01 0. OOOOOOE+OO 

9 2 0.140000E-01 0.363320E-02 0.445763E-01 0.442270E+00 

0 . 169429E-01 0.178476E+01 0.416863E-01 0. OOOOOOE+OO 
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10 


1 0.160000E-01 0 . 000000E+00 0.480659E-01 0.473498E+00 

0 . 000000E+00 0 . 170481E + 01 0.445421E-01 0.000000E + 00 

10 2 0 . 160000E-01 0 . 388210E-02 0.480659E-01 0.474097E+00 

0.169429E-01 0.170481E+01 0.445421E-01 0.000000E+00 

11 1 0 . 180000E-01 0 . 000000E+00 0.514106E-01 0.501827E+00 

0 . 000000E+00 0 . 163993E+01 0.472202E-01 0.000000E+00 

11 2 0 . 180000E-01 0 . 411551E-02 0.514106E-01 0.502413E+00 

0.169429E-01 0.163993E+01 0.472202E-01 0.000000E+00 

12 1 0 . 200000E-01 0 . OOOOOOE+OO 0.546367E-01 0.527258E+00 

0 . 000000E+00 0 . 158612E+01 0.497494E-01 0. 000000E+00 

12 2 0 . 200000E-01 0.433594E-02 0.546367E-01 0.527832E+00 

0 . 169429E-01 0.158612E+01 0.497494E-01 0. OOOOOOE+OO 

13 1 0 . 220000E-01 0. OOOOOOE+OO 0.577635E-01 0.550259E+00 

0. OOOOOOE+OO 0 . 154070E+01 0.521512E-01 0. OOOOOOE+OO 

13 2 0 . 220000E-01 0.454528E-02 0.577635E-01 0.550822E+00 

0 . 1 6942 9E-01 0.154070E + 01 0.521512E-01 0. OOOOOOE + OO 

14 1 0 . 240000E-01 0. OOOOOOE+OO 0.608060E-01 0.571195E+00 

0. OOOOOOE+OO 0 . 150180E + 01 0.544426E-01 0. OOOOOOE + OO 

14 2 0 . 240000E-01 0.474499E-02 0.608060E-01 0.571747E+00 

0 . 169429E-01 0.150180E+01 0.544426E-01 0. OOOOOOE+OO 

15 1 0.260000E-01 0. OOOOOOE+OO 0.637758E-01 0.590356E+00 

0. OOOOOOE+OO 0 . 146808E+01 0.566370E-01 0. OOOOOOE+OO 

15 2 0.260000E-01 0.493624E-02 0.637758E-01 0.590897E+00 

0 . 169429E-01 0.146808E+01 0.566370E-01 0. OOOOOOE+OO 

16 1 0.280000E-01 0. OOOOOOE+OO 0.666825E-01 0.607975E+00 

0. OOOOOOE+OO 0 . 143855E + 01 0.587452E-01 0. OOOOOOE + OO 

16 2 0 . 280000E-01 0.511998E-02 0.666825E-01 0.608506E+00 

0. 16942 9E-01 0.143855E+01 0.587452E-01 0. OOOOOOE+OO 

17 1 0 . 300000E-01 0. OOOOOOE+OO 0.695335E-01 0.624245E+00 

0. OOOOOOE+OO 0.14124 6E+01 0.607762E-01 0. OOOOOOE+OO 

17 2 0 . 300000E-01 0.529700E-02 0.695335E-01 0.624767E+00 

0 . 169429E-01 0.141246E+01 0.607762E-01 0. OOOOOOE+OO 

18 1 0 . 320000E-01 0. OOOOOOE + OO 0.723352E-01 0.639327E+00 

0. OOOOOOE+OO 0 . 138922E+01 0.627375E-01 0. OOOOOOE+OO 

18 2 0 . 320000E-01 0.546794E-02 0.723352E-01 0.639840E+00 

0 . 169429E-01 0.138922E+01 0.627375E-01 0. OOOOOOE+OO 

19 1 0 . 340000E-01 0. OOOOOOE+OO 0.750928E-01 0.653354E + 00 

0. OOOOOOE+OO 0 . 136840E+01 0.646355E-01 0. OOOOOOE + OO 

19 2 0 . 340000E-01 0.563336E-02 0.750928E-01 0.653859E+00 

0 . 169429E-01 0.136840E+01 0.646355E-01 0. OOOOOOE+OO 

20 1 0 . 360000E-01 0. OOOOOOE+OO 0.778108E-01 0.666441E+00 

0. OOOOOOE+OO 0.134 962E + 01 0.664756E-01 0. OOOOOOE + OO 

20 2 0 . 360000E-01 0.579373E-02 0.778108E-01 0.666938E+00 

0 . 169429E-01 0 . 134962E+01 0.664756E-01 0. OOOOOOE+OO 

0 . 967481E+00 0.959568E+00 0.951434E+00 0.943149E+00 0.934767E+00 

0 . 900964E+0 0 0.867690E + 00 0.835596E+00 0.804939E+00 0.775800E + 00 

0 . 748170E+00 0.721999E+00 0.697215E+00 0.673736E+00 0.651480E+00 

0.630367E+00 0.610318E+00 0.591261E+00 0.573128E+00 0.555857E+00 

0. OOOOOOE + OO 0. OOOOOOE + OO 0. OOOOOOE+OO 0. OOOOOOE+OO 0. OOOOOOE+OO 

0. OOOOOOE+OO 0. OOOOOOE+OO 0. OOOOOOE+OO 0. OOOOOOE+OO 0. OOOOOOE+OO 

0. OOOOOOE + OO 0. OOOOOOE + OO 0. OOOOOOE + OO 0. OOOOOOE+OO 0. OOOOOOE + OO 

0. OOOOOOE + OO 0. OOOOOOE+OO 0. OOOOOOE + OO 0. OOOOOOE+OO 0. OOOOOOE + OO 

0 . 100000E+01 0 . 100000E+01 0.100000E+01 0.100000E+01 0.100000E+01 

0 . 100000E+01 0 . 100000E+0 1 0.100000E+01 0.100000E+01 0.100000E+01 

0 . 100000E+0 1 0 . 100000E+01 0.100000E+01 0.100000E+01 0.100000E+01 

0 . 100000E+0 1 0 . 100000E+01 0.100000E+01 0.100000E+01 0.100000E+01 
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6.6 Ellipsoid ol revolution along leeward plane of symmetry 

This is reduced input, file lrom which was used to generate the result presented in PART I, 
section 5.2.3. 


Ellipsoid of revolution, leeward plane of symmetry 

istop mks kshale ktw krow 
20 0 0 00 

pr rminf pinf tinf ukmaxl zetae 

0.72 0.05 2116. 520.0 0.99999 6 0 

imax 
20 


dzetal 

0.1 


rk 

1.0 


0 . 400000E-02 
0 . 800000E-02 
0 . 180000E-01 
0 .280000E-01 
0.314159E+01 


0 . 4 5000 OE-02 
0 . 100000E-01 
0 .200000E-01 
0 . 300000E-01 
0 . 305433E+01 


0 . 500000E-02 
0 . 120000E-01 
0.220000E-01 
0 . 320000E-01 


0.550000E-02 
0.140000E-01 
0 . 240000E-01 
0 .340000E-01 


0 . 600000E-02 
0 . 160000E-01 
0 .260000E-01 
0 .360000E-01 


1 37 0 . 400000E-02 -0 . 7 66260E-1 6 

-0 . 666833E-15 0.296069E+01 

1 36 0.400000E-02 0.194691E-02 

0.169429E-01 0.296069E+01 

2 37 0.450000E-02 -0 . 8 12 63 9E-1 6 

-0 . 666833E-15 0.281027E+01 

2 36 0 . 450000E-02 0.206475E-02 

0.169429E-01 0.281027E+01 

3 37 0.500000E-02 -0 . 856490E-16 

-0 . 666833E-15 0.268387E+01 

3 36 0 . 500000E-02 0.217617E-02 

0.169429E-01 0.268387E+01 

4 37 0 . 550000E-02 -0 . 898181E-16 

-0 . 666833E-15 0.257585E+01 

4 36 0 . 550 000E-02 0.228210E-02 

0.169429E-01 0.257585E+01 

5 37 0.600000E-02 -0 . 938 002E-1 6 

-0 . 666833E-15 0.248224E+01 

5 36 0.600000E-02 0.238328E-02 

0 . 169429E-01 0.248224E+01 

6 37 0.800000E-02 -0 . 108257E-15 

-0 . 666833E-15 0.220441E+01 

6 36 0.800000E-02 0.275059E-02 

0 . 169429E-01 0.220441E+01 

7 37 0 . 100000E-01 -0 . 120974E-15 

-0 . 666833E-15 0.201946E+01 

7 36 0.100000E-01 0.307371E-02 

0.169429E-01 0.201946E+01 

8 37 0 . 120000E-01 -0 . 132454E-15 

-0 . 666833E-15 0.188610E+01 

8 36 0 . 120000E-0 1 0.336539E-02 

0 . 169429E-01 0.188610E+01 

9 37 0.140000E-01 -0 . 142994E-15 

-0 . 666833E-15 0.178476E+01 

9 36 0 . 140000E-01 0.363320E-02 

0 . 169429E-01 0.178476E+01 


0.226936E-01 0.546279E+00 

0 . 223383E-01 O.OOOOOOE+OO 
0 . 226936E-01 0.545583E+00 

0.223383E-01 0.000000E+00 

0.241363E-01 0.564425E+00 

0 . 236904E-01 0.000000E+00 

0.241363E-01 0.563734E+00 

0 . 236904E-01 0.000000E+00 

0.255099E-01 0.581176E+00 

0 . 249687E-01 O.OOOOOOE+OO 
0.255099E-01 0.580490E+00 

0 . 249687E-01 O.OOOOOOE+OO 
0.268248E-01 0.596734E+00 

0.261841E-01 0.000000E+00 

0 . 268248E-01 0.596053E+00 

0 . 261841E-01 0.000000E + 00 

0.280893E-01 0.611256E+00 

0 . 273450E-01 O.OOOOOOE+OO 
0 . 280893E-01 0.610579E+00 

0 . 273450E-01 0.000000E+00 

0.327760E-01 0.661190E+00 

0 . 315595E-01 0.000000E+00 

0.327760E-01 0.660531E+00 

0 . 315595E-01 0.000000E+00 

0 . 369999E-01 0.701525E+00 

0.352668E-01 0.000000E+00 

0.369999E-01 0.700883E+00 

0 . 352668E-01 0.000000E+00 

0 . 409054E-01 0.735119E+00 

0 . 386135E-01 0.000000E+00 

0 . 4 09054E-0 1 0.734492E+00 

0 . 386135E-01 0.000000E+00 

0 . 4 4 5 763E-0 1 0.763693E + 00 

0 . 416863E-01 O.OOOOOOE+OO 
0.445763E-01 0.763081E+00 

0.416863E-01 0.000000E+00 
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10 37 0.160000E-01 -0 . 152791E-15 0.480659E-01 0.788382E+00 

-0 . 666833E-15 0.170481E+01 0.445421E-01 0.000000E+00 

10 36 0 . 160000E-01 0.388210E-02 0.480659E-01 0.787783E+00 

0 . 169429E-01 0.170481E+01 0.445421E-01 0.000000E+00 

11 37 0.180000E-01 -0 . 1 61 977E-1 5 0.514106E-01 0.809974E+00 

-0 . 666833E-15 0.163993E+01 0.472202E-01 0.000000E+00 

11 36 0 . 180000E-01 0.411551E-02 0.514106E-01 0.809388E+00 

0 . 1 6942 9E-01 0.163993E+01 0.472202E-01 0.000000E+00 

12 37 0 . 200000E-01 -0 . 1 70 653E-1 5 0.546367E-01 0.829047E+00 

-0 . 666833E-15 0.158612E+01 0.497494E-01 0.000000E+00 

12 36 0.200000E-01 0.433594E-02 0.546367E-01 0.828473E+00 

0 . 1 6942 9E-01 0.158612E+01 0.497494E-01 O.OOOOOOE+OO 

13 37 0 . 220000E-01 -0 . 178892E-15 0.577635E-01 0.846032E+00 

-0 . 666833E-15 0.154070E+01 0.521512E-01 0.000000E+00 

13 36 0 . 220000E-01 0.454528E-02 0.577635E-01 0.845469E+00 

0 . 169429E-01 0.154070E+01 0.521512E-01 O.OOOOOOE+OO 

14 37 0 . 240000E-01 -0 . 186752E-15 0.608060E-01 0.861264E+00 

-0 . 666833E-15 0.150180E+01 0.544426E-01 O.OOOOOOE+OO 

14 36 0 . 240000E-01 0.474499E-02 0.608060E-01 0.860712E+00 

0 . 169429E-01 0.150180E+01 0.544426E-01 O.OOOOOOE+OO 

15 37 0 . 2 60000E-01 -0 . 194279E-15 0.637758E-01 0.875005E+00 

-0 .666833E-15 0.146808E+01 0.566370E-01 O.OOOOOOE+OO 

15 36 0 . 260000E-01 0.493624E-02 0.637758E-01 0.874464E+00 

0 . 169429E-01 0.146808E+01 0.566370E-01 O.OOOOOOE+OO 

16 37 0 . 280000E-01 -0 . 201511E-15 0.666825E-01 0.887468E+00 

-0 . 666833E-15 0.143855E+01 0.587452E-01 O.OOOOOOE+OO 

16 36 0 . 280000E-01 0.511998E-02 0.666825E-01 0.886936E+00 

0 . 169429E-01 0.143855E+01 0.587452E-01 O.OOOOOOE+OO 

17 37 0 . 300000E-01 -0 . 208 4 7 8E-15 0.695335E-01 0.898823E+00 

-0 . 666833E-15 0.141246E+01 0.607762E-01 O.OOOOOOE+OO 

17 36 0 . 300000E-01 0.529700E-02 0.695335E-01 0.898300E+00 

0.169429E-01 0 . 141246E+01 0.607762E-01 O.OOOOOOE+OO 

18 37 0 . 320000E-01 -0 . 215205E-15 0.723352E-01 0.909211E+00 

-0 . 666833E-15 0.138922E+01 0.627375E-01 O.OOOOOOE+OO 

18 36 0 . 320000E-01 0.546794E-02 0.723352E-01 0.908697E+00 

0 . 169429E-01 0.138922E+01 0.627375E-01 O.OOOOOOE+OO 

19 37 0 . 340000E-01 -0 . 22171 6E-15 0.750928E-01 0.918751E+00 

-0 . 666833E-15 0.136840E+01 0.646355E-01 O.OOOOOOE+OO 

19 36 0.340000E-01 0.563336E-02 0.750928E-01 0.918246E+00 

0.169429E-01 0.136840E+01 0.646355E-01 O.OOOOOOE+OO 

20 37 0 . 360000E-01 -0 .228028E-15 0.778108E-01 0.927540E+00 

-0 . 666833E-15 0.134962E+01 0.664756E-01 O.OOOOOOE+OO 

20 36 0 . 360000E-01 0.579373E-02 0.778108E-01 0.927044E+00 

0 . 169429E-01 0.134962E+01 0.664756E-01 O.OOOOOOE+OO 

0 . 701579E+00 0.681424E+00 0.662234E+00 0.643908E+00 0.626366E+00 

0 . 562828E+00 0.507862E+00 0.459600E+00 0.416773E+00 0.378454E+00 

0 . 343942E+00 0.312681E+00 0.284230E+00 0.258225E+00 0.234365E+00 

0.212401E+00 0 . 192118E+00 0.173336E+00 0.155897E+00 0.139669E+00 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 

0.100000E+01 0 . 100000E+01 0.100000E+01 0.100000E+01 0.100000E+01 

0 . 100000E+01 0 . 100000E+01 0.100000E+01 0.100000E+01 0.100000E+01 

0 . 100000E+01 0 . 100000E+01 0.100000E+01 0 . 100000E+01 0.100000E+01 

0 . 100000E+01 0 . 100000E+01 0.100000E+01 0.100000E+01 0 . 100000E+01 
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6.7 Infinite swept cylider along the leading-edge attachment line 


This input file was used to generate the result presented in PART I, section 5.2.4. 

Infinite swept cylinder leading-edge attachment line 
3 

istop mks kshale ktw krow 
4 0 0 0 0 

pr rminf pinf tinf ukmaxl zetae dzetal rk 

0.72 0.126675 2116.8 520.0 0.99999 3.5 0.001 1.05 

imax 

10 

0 . 100000E-01 0 . 200000E-01 0.300000E-01 0.400000E-01 0.500000E-01 
0 . 600000E-0 1 0 . 700000E-01 0.800000E-01 0.900000E-01 0.100000E+00 
0 . 000000E+00 0 . 174533E-01 

1 1 0 . 000000E+00 0 . 000000E+00 0.100000E-01 0.866025E+00 

0 . 359535E-11 0.100000E+01 0.375000E+00 0.000000E+00 

1 2 0.000000E+00 0 . 000000E+00 0.100000E-01 0.866025E+00 

0 . 174547E-01 0.100000E+01 0.375000E+00 O.OOOOOOE+OO 

2 1 0.000000E+00 0.000000E+00 0.200000E-01 0.866025E+00 

0 . 359535E-1 1 0.100000E+01 0.375000E+00 0.000000E+00 

2 2 O.OOOOOOE+OO 0.000000E+00 0.200000E-01 0.866025E+00 

0 . 17454 7E-01 0.100000E+01 0.375000E+00 0.000000E+00 

3 1 0.000000E+00 0.000000E+00 0.300000E-01 0.866025E+00 

0 . 359535E-11 0.100000E+01 0.375000E+00 0.000000E+00 

3 2 0.000000E+00 0.000000E+00 0.300000E-01 0.866025E+00 

0.17 4 54 7E-01 0.100000E + 01 0.375000E+00 0.000000E + 00 

4 1 O.OOOOOOE+OO 0.000000E+00 0.400000E-01 0.866025E+00 

0 . 359535E-11 0.100000E+01 0.375000E+00 0.000000E+00 

4 2 0.000000E+00 O.OOOOOOE+OO 0.400000E-01 0.866025E+00 

0 . 17454 7E-01 0.100000E+01 0.375000E+00 0.000000E+00 

5 1 0.000000E+00 O.OOOOOOE+OO 0.500000E-01 0.866025E+00 

0 . 359535E-11 0.100000E+01 0.375000E+00 0.000000E+00 

5 2 0.000000E+00 O.OOOOOOE+OO 0.500000E-01 0.866025E+00 

0 . 17454 7E-01 0.100000E+01 0.375000E+00 O.OOOOOOE+OO 

6 1 O.OOOOOOE+OO O.OOOOOOE+OO 0.600000E-01 0.866025E+00 

0 . 359535E-11 0.100000E+01 0.375000E+00 O.OOOOOOE+OO 

6 2 O.OOOOOOE+OO O.OOOOOOE+OO 0.600000E-01 0.866025E+00 

0 . 174547E-01 0.100000E+01 0.375000E+00 O.OOOOOOE+OO 

7 1 O.OOOOOOE+OO O.OOOOOOE+OO 0.700000E-01 0.866025E+00 

0 . 359535E-11 0.100000E+01 0.375000E+00 O.OOOOOOE+OO 

7 2 O.OOOOOOE+OO O.OOOOOOE+OO 0.700000E-01 0.866025E+00 

0 . 174547E-01 0.100000E+01 0.375000E+00 O.OOOOOOE+OO 

8 1 O.OOOOOOE+OO O.OOOOOOE+OO 0.800000E-01 0.866025E+00 

0 . 359535E-11 0.100000E+01 0.375000E+00 O.OOOOOOE+OO 

8 2 O.OOOOOOE+OO O.OOOOOOE+OO 0.800000E-01 0.866025E+00 

0 . 174547E-01 0.100000E+01 0.375000E+00 O.OOOOOOE+OO 

9 1 O.OOOOOOE+OO O.OOOOOOE+OO 0.900000E-01 0.866025E+00 

0 . 359535E-11 0.100000E+01 0.375000E+00 O.OOOOOOE+OO 

9 2 O.OOOOOOE+OO O.OOOOOOE+OO 0.900000E-01 0.866025E+00 

0.174 547E-01 0.100000E + 01 0.375000E + 00 O.OOOOOOE+OO 

10 1 O.OOOOOOE+OO O.OOOOOOE+OO 0.100000E+00 0.866025E+00 

0 . 359535E-1 1 0.100000E+01 0.375000E+00 O.OOOOOOE+OO 

10 2 O.OOOOOOE+OO O.OOOOOOE+OO 0.100000E+00 0.866025E+00 

0.17 4 54 7E-0 1 0.100000E+01 0.375000E+00 O.OOOOOOE + OO 
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0 .250251E+00 
0 .250251E+00 
0 .OOOOOOE+OO 
0 . OOOOOOE+OO 
0 .100000E+01 
0 .100000E+01 


0 .250251E+00 
0 .250251E+00 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 . lOOOOOE+Ol 
0 . lOOOOOE+Ol 


0 .250251E+00 
0 .250251E+00 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 .lOOOOOE+Ol 
0 . lOOOOOE+Ol 


0 .250251E+00 
0 .250251E+00 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 .lOOOOOE+Ol 
0 .lOOOOOE+Ol 


0 .250251E+00 
0 .250251E+00 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 . lOOOOOE+Ol 
0 .lOOOOOE+Ol 
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6.8 Svvept-back wing 


This input file was used to generate the result presented in PART I, section 5.2.5. 


F-14 VSTFE CLEAN-UP GLOVE 
4 

0 10 0 

51 48 0.20000 0.0010000 1.05 

5.425594 22.18 4.68 .80 .0 785.30000 430.00000 


.1616606 

-.0262935 

.1414880 

-.0257465 

.1213054 

-.0250774 

.1011129 

-.0242459 

.0809103 

-.0231646 

.0606977 

-.0216108 

.0404044 

-.0194493 

.0202426 

-.0156186 

.0100719 

-.0121270 

.0050211 

-.0091629 

.0019336 

-.0058652 

.0000000 

.0000000 

.0019336 

. 0061704 

.0050211 

.0097646 

.0100719 

. 0136796 

.0202426 

. 0192767 

.0404044 

. 0271152 

.0606977 

.0329864 

.0809103 

. 0379677 

. 1011129 

. 0422232 

.1213054 

. 0459822 

.1414880 

.0493646 

.1616606 

.0524201 

.1818233 

.0551935 

.2019760 

.0576876 

.2221187 

.0599615 

.2422515 

.0619935 

.2623744 

.0638032 

.2824873 

. 0653983 

.3025903 

.0667787 

.3226834 

.0679643 

.3427665 

.0688854 

.3628398 

.0697315 

.3829032 

.0703230 

.4029567 

. 0707373 

.4230003 

. 0709645 

.4430340 

.0710145 

.4630579 

.0708873 

.4830719 

.0705982 

.5030761 

.0701466 

.5230704 

.0695354 

.5630296 

.0678540 

. 6029494 

.0656117 

.6428300 

. 0628557 

.6826715 

. 0596285 

.7025776 

. 0577982 
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. 72 


i 


7224739 

.0559902 

7423604 

.0540278 

7622372 

.0519857 

7821042 

.0498664 

8019615 

. 0476850 

0001219 

. 9828369 

0008643 

. 9291969 

0016067 

. 8558530 

0023491 

.7860392 

0030916 

.7235947 

0038340 

. 6685212 

0060612 

.5306466 

0082885 

. 414 6636 

0105158 

. 3174372 

0127431 

.2479026 

0149703 

.1841853 

0171976 

. 1250204 

0194249 

. 0695431 

0246218 

- . 0133144 

0298188 

-.0815164 

0350158 

-.1439954 

0402127 

-.2019712 

0454097 

-.2352954 

0506067 

-.2656347 

0558036 

-.2944504 

0669400 

-.3357120 

0780763 

- . 3613739 

0892127 

-.3837189 

1003490 

-.4043051 

1114854 

-.4173073 

1226217 

-.4298466 

1337581 

-.4457382 

1448944 

-.4610751 

1656823 

-.4949900 

1864701 

-.5317120 

2072580 

-.5662568 

,2280459 

- . 5959700 

,2488337 

-.6247629 

,2696216 

-.6515539 

,2904094 

-.6774718 

.3111973 

-.7032281 

.3319851 

-.7290463 

.3527730 

-.7545299 

.3787578 

-.7868917 

.4047426 

-.8187169 

.4307275 

-.8478254 

.4567123 

-.8746795 

.4826971 

-.8885065 

.5086819 

-.8964149 

. 5230704 

-.8887564 

.5391213 

-.8802513 

.5695606 

-.8642401 

.6000000 

-.8483981 


.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1 .000000 

-.0007000 

1.000000 

-.0007000 

1.000000 

-.0007000 

1 .000000 

- . 0007000 

1.000000 

-.0007000 

1.000000 

- .0007000 

1.000000 

-.0007000 

1.000000 

- .0007000 

1.000000 

-.0007000 

1.000000 

-.0007000 

1.000000 

-.0007000 

1.000000 

-.0007000 

1.000000 

- .0007000 

1.000000 

-.0007000 

1.000000 

-.0007000 

1.000000 

-.0007000 

1.000000 

-.0007000 

1.000000 

- .0007000 

1.000000 

-.0007000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 

.0000000 

1.000000 
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7. SAMPLE CASE OUTPUT 


To reduce the size of output data, only the boundary-layer parameters written on file 
FORT. 30 will be shown. 


7.1 Flat plate 


i x/c 

cf y 


cf x 
dspty 


blth 

thmoy 


dsptx 

twall 


thmox 

qw 


1 0.000000E+00 
0.000000E+00 

2 0.100000E-01 
0.000000E+00 

3 0 . 200000E-01 
0 . 000000E+00 

4 0.300000E-01 
O.OOOOOOE+OO 

5 0 . 400000E-01 
O.OOOOOOE+OO 

6 0 . 500000E-0 1 
0.000000E+00 

7 0 . 600000E-01 
O.OOOOOOE+OO 

8 0 . 700000E-01 
O.OOOOOOE+OO 

9 0 . 800000E-01 
0.000000E+00 

10 0 . 900000E-01 
0.000000E+00 

11 0 . 100000E+00 
0 . 000000E+00 

12 0 . 110000E+00 
0.000000E+00 

13 0 . 120000E+00 
0 . OOOOOOE+OO 

14 0 . 130000E+00 
O.OOOOOOE+OO 

15 0 . 140000E+00 
O.OOOOOOE+OO 

16 0 . 150000E+00 
O.OOOOOOE+OO 

17 0 . 160000E+00 
O.OOOOOOE+OO 

18 0 . 170000E+00 
0 . 000000E+00 

19 0 . 180000E+00 
O.OOOOOOE+OO 

20 0.190000E+00 
O.OOOOOOE+OO 

21 0 . 200000E+00 
O.OOOOOOE+OO 


0 . 000000E+00 
O.OOOOOOE+OO 
0 . 484239E-03 
O.OOOOOOE+OO 
0 .342408E-03 
O.OOOOOOE+OO 
0 .279575E-03 
O.OOOOOOE+OO 
0 . 242119E-03 
O.OOOOOOE+OO 
0 .216558E-03 
O.OOOOOOE+OO 
0 .197690E-03 
O.OOOOOOE+OO 
0 .183025E-03 
O.OOOOOOE+OO 
0 .171204E-03 
0 . OOOOOOE+OO 
0 . 161413E-03 
O.OOOOOOE+OO 
0 .153130E-03 
O.OOOOOOE+OO 
0 .146003E-03 
O.OOOOOOE+OO 
0 .139788E-03 
O.OOOOOOE+OO 
0 . 134304E-03 
O.OOOOOOE+OO 
0 . 129418E-03 
O.OOOOOOE+OO 
0 . 125030E-03 
O.OOOOOOE+OO 
0 . 121060E-03 
O.OOOOOOE+OO 
0 . 117445E-03 
O.OOOOOOE+OO 
0 . 114136E-03 
O.OOOOOOE+OO 
0 . 111092E-03 
O.OOOOOOE+OO 
0 . 108279E-03 
0 .OOOOOOE+OO 


O.OOOOOOE+OO 
O.OOOOOOE+OO 
0 . 517312E-04 
O.OOOOOOE+OO 
0 . 731 58 9E-04 
O.OOOOOOE+OO 
0 . 896010E-04 
0 . OOOOOOE+OO 
0 . 1034 62E-03 
O.OOOOOOE+OO 
0 . 115674E-03 
O.OOOOOOE+OO 
0 . 126715E-03 
O.OOOOOOE+OO 
0 .136868E-03 
O.OOOOOOE+OO 
0 . 146318E-03 
O.OOOOOOE+OO 
0 .155194E-03 
O.OOOOOOE+OO 
0 . 163588E-03 
O.OOOOOOE+OO 
0 . 171573E-03 
O.OOOOOOE+OO 
0 . 179202E-03 
O.OOOOOOE+OO 
0.186519E-03 
O.OOOOOOE+OO 
0 . 193560E-03 
O.OOOOOOE+OO 
0 .200354E-03 
0 .OOOOOOE+OO 
0 .206925E-03 
O.OOOOOOE+OO 
0 .213293E-03 
O.OOOOOOE+OO 
0 .219477E-03 
O.OOOOOOE+OO 
0 .225491E-03 
O.OOOOOOE+OO 
0 .231349E-03 
O.OOOOOOE+OO 


O.OOOOOOE+OO 
O.OOOOOOE+OO 
0 .251080E-04 
0 . 868923E+03 
0 . 355080E-0 4 
0 . 868923E+03 
0 . 434 883E-0 4 
0 . 868923E+03 
0 . 502160E-04 
0 .868923E+03 
0 .561431E-04 
0 .868923E+03 
0 .615017E-04 
0 .868923E+03 
0 . 664295E-0 4 
0 .868923E+03 
0 .710161E-04 
0 . 868923E+03 
0 . 753239E-04 
0.868923E+03 
0 .793984E-04 
0 . 868923E+03 
0 .832737E-04 
0 .868923E+03 
0 .869766E-04 
0 .868923E+03 
0 . 905281E-04 
0 .868923E+03 
0 .939455E-04 
0 .868923E+03 
0 . 972428E-04 
0 .868923E+03 
0 .100432E-03 
0 .868923E+03 
0.103523E-03 
0 . 868923E+03 
0 . 106524E-03 
0 . 868923E+03 
0.109443E-03 
0 . 868923E+03 
0 . 112286E-03 
0.868923E+03 


O.OOOOOOE+OO 
O.OOOOOOE+OO 
0 .483610E-05 
-0 . 150288E+01 
0 . 683927E-05 
-0 . 106270E+01 
0 .837636E-05 
-0 .867689E+00 
0 . 967219E-05 
-0 .75144 1E+00 
0 . 108138E-04 
-0 . 672109E+00 
0 . 118460E-04 
-0 .613549E+00 
0 . 127951E-04 
-0 . 568036E+00 
0 . 136785E-04 
-0 .531349E+00 
0 .145083E-04 
-0 . 500960E+00 
0 . 152931E-04 
-0 .475253E+00 
0 . 160395E-04 
-0 . 453136E+00 
0 . 167527E-04 
-0 .433844E+00 
0 .174368E-04 
-0 . 416824E+00 
0 . 180950E-04 
-0 . 401662E+00 
0 . 187301E-04 
-0 . 388042E+00 
0 . 193444E-04 
-0.375720E+00 
0 . 199397E-04 
-0 .364502E+00 
0 .205178E-04 
-0 .354233E+00 
0 . 210801E-04 
-0 . 344785E+00 
0 . 216277E-04 
-0 . 336054E+00 
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I 


7.2 NACA 0012 airfoil 


i x/c 

cf y 


cfx blth 

dspty thmoy 


dsptx 

twall 


1 0 . 000000E+00 
O.OOOOOOE+OO 

2 0 . 775000E-03 
O.OOOOOOE+OO 

3 0 . 295500E-02 
0.000000E+00 

4 0 . 621500E-02 
O.OOOOOOE+OO 

5 0 . 101680E-01 
0.000000E+00 

6 0 . 146000E-01 
O.OOOOOOE+OO 

7 0 . 193920E-01 
O.OOOOOOE+OO 

8 0 . 244 530E-01 
0.000000E+00 

9 0 . 297460E-01 
O.OOOOOOE+OO 

10 0 . 352350E-01 
0 . 000000E+00 

11 0 . 408940E-01 
O.OOOOOOE+OO 

12 0 . 467010E-01 
O.OOOOOOE+OO 

13 0 . 526440E-01 
O.OOOOOOE+OO 

14 0 . 587070E-01 
O.OOOOOOE+OO 

15 0 . 64 8830E-0 1 
O.OOOOOOE+OO 

16 0.711610E-01 
O.OOOOOOE+OO 

17 0.775270E-01 
O.OOOOOOE+OO 

18 0 . 839790E-01 
O.OOOOOOE+OO 

19 0 . 905060E-01 
O.OOOOOOE+OO 

20 0 . 971000E-01 
O.OOOOOOE+OO 

21 0 . 103757E+00 
O.OOOOOOE+OO 

22 0 . 117234E+00 
O.OOOOOOE+OO 

23 0 . 130892E+00 
O.OOOOOOE+OO 

24 0 . 144698E+00 
O.OOOOOOE+OO 

25 0 . 158623E+00 
O.OOOOOOE+OO 


O.OOOOOOE+OO 
0 . 000000E+00 
0 . 495157E-01 
O.OOOOOOE+OO 
0 . 250882E-01 
O.OOOOOOE+OO 
0 . 166742E-01 
O.OOOOOOE+OO 
0 . 121680E-01 
O.OOOOOOE+OO 
0 . 941305E-02 
O.OOOOOOE+OO 
0 .766982E-02 
O.OOOOOOE+OO 
0 . 650778E-02 
O.OOOOOOE+OO 
0 .560452E-02 
O.OOOOOOE+OO 
0 . 488171E-02 
O.OOOOOOE+OO 
0 . 437027E-02 
0 . OOOOOOE+OO 
0 .399905E-02 
O.OOOOOOE+OO 
0 .367386E-02 
O.OOOOOOE+OO 
0 . 337700E-02 
0 .OOOOOOE+OO 
0 . 312142E-02 
O.OOOOOOE+OO 
0 .290185E-02 
0 .OOOOOOE+OO 
0 . 272025E-02 
O.OOOOOOE+OO 
0 . 257701E-02 
O.OOOOOOE+OO 
0 .245206E-02 
O.OOOOOOE+OO 
0 .232289E-02 
O.OOOOOOE+OO 
0 .219169E-02 
O.OOOOOOE+OO 
0 . 198242E-02 
O.OOOOOOE+OO 
0 . 182297E-02 
O.OOOOOOE+OO 
0 . 167648E-02 
O.OOOOOOE+OO 
0 . 154640E-02 
O.OOOOOOE+OO 


O.OOOOOOE+OO 
0 . OOOOOOE+OO 
0 . 852824E-04 
0 . OOOOOOE+OO 
0 . 102382E-03 
0 . OOOOOOE+OO 
0 . 119106E-03 
O.OOOOOOE+OO 
0 . 138414E-03 
O.OOOOOOE+OO 
0 . 159176E-03 
O.OOOOOOE+OO 
0 . 179854E-03 
0 .OOOOOOE+OO 
0 . 200040E-03 
0 .OOOOOOE+OO 
0 .220015E-03 
0 .OOOOOOE+OO 
0 .240048E-03 
0 .OOOOOOE+OO 
0 . 259342E-03 
0 .OOOOOOE+OO 
0 .277574E-03 
0 .OOOOOOE+OO 
0 .295311E-03 
O.OOOOOOE+OO 
0 .313138E-03 
O.OOOOOOE+OO 
0 .330562E-03 
0 .OOOOOOE+OO 
0 . 347783E-03 
O.OOOOOOE+OO 
0.364579E-03 
0 .OOOOOOE+OO 
0 .380758E-03 
O.OOOOOOE+OO 
0 . 3 96351E-03 
O.OOOOOOE+OO 
0 . 412124E-03 
O.OOOOOOE+OO 
0 . 427771E-03 
O.OOOOOOE+OO 
0 . 458722E-03 
0 .OOOOOOE+OO 
0 .488208E-03 
0 .OOOOOOE+OO 
0 . 517222E-03 
0 .OOOOOOE+OO 
0 . 545901E-03 
O.OOOOOOE+OO 


O.OOOOOOE+OO 
O.OOOOOOE+OO 
0 . 215699E-04 
0 . 545375E+03 - 
0 . 266100E-04 
0 . 544426E+03 - 
0.318478E-04 
0 . 54340 6E+03 - 
0 .381869E-04 
0 . 542555E+03 - 
0 .452794E-04 
0 . 541932E+03 - 
0 .525144E-04 
0 . 541490E+03 - 
0 .595256E-04 
0 . 541160E+03 - 
0 . 66634 7E-0 4 
0 . 540912E+03 - 
0 .739434E-04 
0 . 540734E+03 - 
0 .809198E-04 
0 . 540 601E+0 3 - 
0 .873951E-04 
0 . 540491E+03 - 
0 . 937532E-04 
0 . 540399E+03 - 
0.100198E-03 
0 . 540327E+03 - 
0 . 106644E-03 
0 . 540271E+03 - 
0 . 113038E-03 
0 . 540228E+03 - 
0 . 119258E-03 
0 . 540195E+03 - 
0 . 125154E-03 
0 . 540167E+03 - 
0 . 130836E-03 
0 . 540143E+03 - 
0 .136626E-03 
0 . 540123E+03 - 
0 . 142646E-03 
0 . 540110E+03 - 
0 .154454E-03 
0 . 540099E+03 - 
0 .165712E-03 
0 .540097E+03 
0 .176968E-03 
0 .540104E+03 
0 .188243E-03 
0 .540118E+03 


thmox 

qw 

O.OOOOOOE+OO 
0 .OOOOOOE+OO 
0 . 952800E-05 
0 . 729836E+01 
0 . 114849E-04 
0 . 130195E+02 
0 . 134165E-04 
0 . 148705E+02 
0 . 156829E-04 
0 . 121798E+02 
0 . 181766E-04 
0 . 877449E+01 
0 .207107E-04 
0 . 642159E+01 
0 .231665E-04 
0 . 497353E+01 
0 . 256091E-04 
0 .363998E+01 
0 .280756E-04 
0 .253464E+01 
0 .304633E-04 
0 .201267E+01 
0 .327215E-04 
0 . 175704E+01 
0 . 349129E-04 
0 .143119E+01 
0 .370848E-04 
0 . 108580E+01 
0 .392359E-04 
0 . 826581E+00 
0 .413584E-04 
0 . 629495E+00 
0 .434362E-04 
0 . 505497E+00 
0 . 454434E-04 
0 .456180E+00 
0 . 473884E-04 
0 .403273E+00 
0 . 493140E-04 
0 .289249E+00 
0 .512520E-04 
0 .156901E+00 
0 . 550697E-04 
0 .379250E-01 
0 .587479E-04 
0 .133113E-01 
0 .623528E-04 
0 . 910654E-01 
0 . 659050E-04 
0 .151775E+00 
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26 0 . 172644E+00 
0 .OOOOOOE+OO 

27 0.186741E+00 
0 .OOOOOOE+OO 

28 0 . 200899E+00 
0 .OOOOOOE+OO 

29 0 . 215106E+00 
0 .OOOOOOE+OO 

30 0 . 229351E+00 
0. OOOOOOE+OO 

31 0 . 243625E+00 
0. OOOOOOE+OO 

32 0 . 257924E+00 
0 .OOOOOOE+OO 

33 0 . 272242E+00 
0. OOOOOOE+OO 

34 0.286573E+00 
0. OOOOOOE+OO 

35 0 . 300 91 8E+00 
0. OOOOOOE+OO 

36 0.315270E+00 
0 .OOOOOOE+OO 

37 0 . 329 62 9E + 00 
0. OOOOOOE+OO 

38 0 . 343992E+00 
0 .OOOOOOE+OO 

39 0 . 358359E+00 
0 .OOOOOOE+OO 

40 0 . 372730E+00 
0 .OOOOOOE+OO 

41 0 . 387 10 OE+O 0 
0 .OOOOOOE+OO 

42 0 . 401472E+00 
0 .OOOOOOE+OO 

43 0 . 4 1584 6E+0 0 
0 .OOOOOOE+OO 

44 0 . 430217E+00 
0 .OOOOOOE+OO 

45 0 . 444588E+00 
0 .OOOOOOE+OO 

46 0 . 458957E+00 
0 .OOOOOOE+OO 

47 0 .47332 6E+00 
0 .OOOOOOE+OO 

48 0 . 487695E+00 
0 .OOOOOOE+OO 


0 . 143449E-02 
0 .OOOOOOE+OO 
0 .132994E-02 
0 .OOOOOOE+OO 
0 .123895E-02 
0 . OOOOOOE+OO 
0 . 114714E-02 
0 .OOOOOOE+OO 
0 .105472E-02 
0 . OOOOOOE+OO 
0 . 998273E-03 
0 . OOOOOOE+OO 
0 . 948 108E-03 
0. OOOOOOE+OO 
0 .872308E-03 
0 .OOOOOOE+OO 
0 .806571E-03 
0. OOOOOOE+OO 
0.754748E-03 
0 .OOOOOOE+OO 
0 . 696915E-03 
0 .OOOOOOE+OO 
0 . 634864E-03 
0. OOOOOOE+OO 
0 . 586320E-03 
0 .OOOOOOE+OO 
0 . 546420E-03 
0. OOOOOOE+OO 
0 .502311E-03 
0 .OOOOOOE+OO 
0 . 455670E-03 
0 .OOOOOOE+OO 
0 .422837E-03 
0 .OOOOOOE+OO 
0 .386999E-03 
0 .OOOOOOE+OO 
0 .321864E-03 
0 .OOOOOOE+OO 
0 .242323E-03 
0 .OOOOOOE+OO 
0 .206234E-03 
0 .OOOOOOE+OO 
0 . 184615E-03 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 


0 .573779E-03 
0 .OOOOOOE+OO 
0 .601694E-03 
0 .OOOOOOE+OO 
0 . 628882E-03 
0. OOOOOOE+OO 
0 . 656482E-03 
0 . OOOOOOE+OO 
0 . 684224E-03 
0 . OOOOOOE+OO 
0 . 710749E-03 
0 . OOOOOOE+OO 
0 . 736188E-03 
0 .OOOOOOE+OO 
0 . 763204E-03 
0 .OOOOOOE+OO 
0.79044 9E-03 
0. OOOOOOE+OO 
0 . 816888E-03 
0. OOOOOOE+OO 
0 . 844 151E-03 
0 .OOOOOOE+OO 
0 .872225E-03 
0 .OOOOOOE+OO 
0 .899933E-03 
0 .OOOOOOE+OO 
0 . 926862E-03 
0 .OOOOOOE+OO 
0 . 954566E-03 
0 .OOOOOOE+OO 
0 . 983070E-03 
0 .OOOOOOE+OO 
0 . 101051E-02 
0 .OOOOOOE+OO 
0 .103841E-02 
0 .OOOOOOE+OO 
0 .107186E-02 
0 .OOOOOOE+OO 
0 . 11122 8E-02 
0 .OOOOOOE+OO 
0 . 114607E-02 
0 .OOOOOOE+OO 
0 . 117 689E-02 
0 .OOOOOOE+OO 
0 .OOOOOOE + OO 
0 .OOOOOOE+OO 


0 . 199399E-03 
0 .540138E+03 
0 .210683E-03 
0 .540164E+03 
0 .221905E-03 
0 .540192E+03 
0 .233532E-03 
0 .540225E+03 
0 .245795E-03 
0 .540263E+03 
0 .256948E-03 
0 .540300E+03 
0 .267649E-03 
0 . 540336E+03 
0 .279913E-03 
0 . 540374E+03 
0 . 292504E-03 
0 .540416E+03 
0.304733E-03 
0 . 540457E+03 
0 . 317695E-03 
0 . 540499E+03 
0 .331716E-03 
0.540543E+03 
0 .345480E-03 
0 . 540587E+03 
0 .358838E-03 
0 .540630E+03 
0 .372916E-03 
0 .540673E+03 
0 .387973E-03 
0 .540716E+03 
0 .402210E-03 
0 . 540758E+03 
0 . 417052E-03 
0 .540799E+03 
0 . 437166E-03 
0 .540842E+03 
0 . 463352E-03 
0 . 540889E+03 
0 .484547E-03 
0 .540933E+03 
0 .502658E-03 
0 .540973E+03 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
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0 . 693890E-04 
0 . 189371E+00 
0 . 728367E-04 
0 .230475E+00 
0 .762361E-04 
0 .252147E+00 
0 .796277E-04 
0 .286715E+00 
0 .830536E-04 
0 .319619E+00 
0 .863605E-04 
0 .290747E+00 
0 .895620E-04 
0 . 279839E+00 
0 . 928435E-04 
0 .315149E+00 
0 .961393E-04 
0 .316246E+00 
0 .993817E-04 
0 .303006E+00 
0 . 102641E-03 
0 .303635E+00 
0 . 105947E-03 
0 .302272E+00 
0 . 109221E-03 
0 . 284722E+00 
0 . 112436E-03 
0 .266901E+00 
0 . 115655E-03 
0 .255835E+00 
0 . 118891E-03 
0 .243281E+00 
0 .122065E-03 
0 .223411E+00 
0 .125216E-03 
0 .209451E+00 
0 .128571E-03 
0 .196891E+00 
0 .132141E-03 
0 . 165820E+00 
0 .135414E-03 
0 . 135557E+00 
0 .138538E-03 
0 .117190E+00 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 


I 


7.3 Hypersonic cone 


i x/c 

cf y 


cf X 

blth 

dsptx 

thmox 

dspty 

thmoy 

twall 

qw 


1 0 . 500000E-02 
0 . OOOOOOE+OO 

2 0 . 600000E-02 
0. OOOOOOE+OO 

3 0 . 700000E-02 
0 .OOOOOOE+OO 

4 0 . 800000E-02 
0 .000000E+00 

5 0 . 900000E-02 
0 .000000E+00 

6 0 . 100000E-01 
0 .OOOOOOE+OO 

7 0 . 120000E-01 
0 .OOOOOOE+OO 

8 0.140000E-01 
0 .OOOOOOE+OO 

9 0 . 160000E-01 
0 .OOOOOOE+OO 

10 0 . 180000E-01 
0 .OOOOOOE+OO 

11 0 . 200000E-01 
0. OOOOOOE+OO 

12 0 . 220000E-01 
0 .OOOOOOE+OO 

13 0.240000E-01 
0 .OOOOOOE+OO 

14 0 . 260000E-01 
0 .OOOOOOE+OO 

15 0 . 280000E-01 
0. OOOOOOE+OO 

16 0 . 300000E-01 
0. OOOOOOE+OO 

17 0.350000E-01 
0 .OOOOOOE+OO 

18 0.400000E-01 
0. OOOOOOE+OO 

19 0 . 450000E-01 
0 .OOOOOOE+OO 

20 0 . 500000E-01 
0 .OOOOOOE+OO 

21 0.550000E-01 
0 .OOOOOOE+OO 

22 0 . 600000E-0 1 
0 .OOOOOOE+OO 

23 0 . 650000E-0 1 
0 .OOOOOOE+OO 

24 0 . 700000E-01 
0. OOOOOOE+OO 

25 0 . 750000E-01 
0 .OOOOOOE+OO 


0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 .628591E-02 
0 .OOOOOOE+OO 
0 . 474071E-02 
0 .OOOOOOE+OO 
0 . 407778E-02 
0 . OOOOOOE+OO 
0 . 367952E-02 
0 .OOOOOOE+OO 
0 .340323E-02 
0 . OOOOOOE+OO 
0 .302557E-02 
0 .OOOOOOE+OO 
0 .276607E-02 
0 .OOOOOOE+OO 
0 .256892E-02 
0 .OOOOOOE+OO 
0 .241072E-02 
0 .OOOOOOE+OO 
0 .227945E-02 
0 .OOOOOOE+OO 
0 .216802E-02 
0 . OOOOOOE+OO 
0 .207179E-02 
0 .OOOOOOE+OO 
0 . 198757E-02 
0 .OOOOOOE+OO 
0 . 191302E-02 
0 .OOOOOOE+OO 
0 . 184 64 3E-02 
0 .OOOOOOE+OO 
0 . 170687E-02 
0 .OOOOOOE+OO 
0 .159549E-02 
0 .OOOOOOE+OO 
0 .150338E-02 
0 .OOOOOOE+OO 
0 .142559E-02 
0 .OOOOOOE+OO 
0 . 135890E-02 
0. OOOOOOE+OO 
0 .130070E-02 
0 .OOOOOOE+OO 
0 .124949E-02 
0 .OOOOOOE+OO 
0 . 120381E-02 
0 .OOOOOOE+OO 
0.116288E-02 
0 .OOOOOOE+OO 


0. OOOOOOE+OO 
0. OOOOOOE+OO 
0 . 130625E-03 
0. OOOOOOE+OO 
0 . 172452E-03 
0. OOOOOOE+OO 
0 . 199896E-03 
0 .OOOOOOE+OO 
0 .221067E-03 
0 .OOOOOOE+OO 
0 .239033E-03 
0 .OOOOOOE+OO 
0 .268638E-03 
0 .OOOOOOE+OO 
0 .293948E-03 
0 .OOOOOOE+OO 
0 .316509E-03 
0 .OOOOOOE+OO 
0 .337936E-03 
0 .OOOOOOE+OO 
0 . 357772E-03 
0 .OOOOOOE+OO 
0 .376294E-03 
0 .OOOOOOE+OO 
0 .393850E-03 
0 .OOOOOOE+OO 
0 . 410864E-03 
0 . OOOOOOE+OO 
0 .427084E-03 
0. OOOOOOE+OO 
0 . 442612E-03 
0. OOOOOOE+OO 
0 .478839E-03 
0. OOOOOOE+OO 
0 . 512077E-03 
0 .OOOOOOE+OO 
0 . 543287E-03 
0 .OOOOOOE+OO 
0 . 573147E-03 
0 .OOOOOOE+OO 
0 . 601427E-03 
0 .OOOOOOE+OO 
0 . 628 4 4 9E-03 
0 .OOOOOOE+OO 
0 . 654282E-03 
0 .OOOOOOE+OO 
0 . 679158E-03 
0. OOOOOOE+OO 
0 .703091E-03 
0 .OOOOOOE+OO 


0 .OOOOOOE+OO 
0 .316999E+03 
0 .100913E-03 
0 .316999E+03 
0 . 133553E-03 
0 . 316999E+03 
0 . 154 997E-03 
0 .316999E+03 
0 .171620E-03 
0 .316999E+03 
0 .185525E-03 
0 . 316999E+03 
0 .208605E-03 
0 .316999E+03 
0 .228179E-03 
0 . 316999E+03 
0 .245801E-03 
0 .316999E+03 
0 .262080E-03 
0 . 316999E+03 
0 . 277311E-03 
0 .316999E+03 
0 .291676E-03 
0 .316999E+03 
0 .305306E-03 
0 . 316999E+03 
0 .318301E-03 
0 .316999E+03 
0 . 33074 4E-03 
0 . 316999E+03 
0 .342701E-03 
0 . 316999E+03 
0 .370748E-03 
0 .316999E+03 
0 .396637E-03 
0 .316999E+03 
0 . 420941E-03 
0 .316999E+03 
0 .443917E-03 
0 . 316999E+03 
0 .465707E-03 
0 .316999E+03 
0 .486554E-03 
0 . 316999E+03 
0 .506501E-03 
0 .316999E+03 
0 .525724E-03 
0 .316999E+03 
0 .544229E-03 
0 . 316999E+03 


0 . OOOOOOE+OO 
0. OOOOOOE+OO 
0 . 538282E-05 
0 . 101339E+06 
0 . 711912E-05 
0 . 764025E+05 
0 . 825610E-05 
0 . 657031E+05 
0 . 913739E-05 
0 . 592819E+05 
0 . 987604E-05 
0 . 548330E+05 
0 . 111012E-04 
0 . 487516E+05 
0 . 121416E-04 
0 . 445744E+05 
0 . 130806E-04 
0 .414023E+05 
0 . 139495E-04 
0 .388572E+05 
0 . 147 630E-04 
0 .367445E+05 
0 . 155302E-04 
0 .349501E+05 
0 . 162579E-04 
0 .333998E+05 
0 .169515E-04 
0 .320423E+05 
0 .176154E-04 
0 .308406E+05 
0 .182531E-04 
0 .297668E+05 
0 . 197485E-04 
0 . 275162E+05 
0 . 211283E-04 
0 .257198E+05 
0 . 224236E-04 
0 .24234 6E+05 
0 .236481E-04 
0 .229803E+05 
0 .248093E-04 
0 .219052E + 05 
0 .259202E-04 
0 .209668E+05 
0 .269832E-04 
0 .201412E+05 
0 .280075E-04 
0 . 194048E+05 
0 .289936E-04 
0.187450E+05 
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26 0 . 800000E-01 
0 . OOOOOOE+OO 

27 0.850000E-01 
0 .OOOOOOE+OO 

28 0 . 900000E-01 
0. OOOOOOE+OO 

29 0 . 950000E-01 
0. OOOOOOE+OO 

30 0 . 97000 OE-0 1 
0. OOOOOOE+OO 

31 0 . 99000 OE-0 1 
0. OOOOOOE+OO 

32 0 . 101000E+00 
0 .OOOOOOE+OO 

33 0 . 103000E+00 
0 .OOOOOOE+OO 

34 0 . 105000E+00 
0 .OOOOOOE+OO 

35 0 . 107000E+00 
0 .OOOOOOE+OO 

36 0 . 109000E+00 
0 .OOOOOOE+OO 

37 0 . 111000E+00 
0. OOOOOOE+OO 

38 0 . 113000E+00 
0 .OOOOOOE+OO 

39 0.115000E+00 
0. OOOOOOE+OO 

40 0 . 125000E+00 
0 .OOOOOOE+OO 

41 0 . 135000E+00 
0 .OOOOOOE+OO 

42 0 . 145000E+00 
0. OOOOOOE+OO 

43 0 . 155000E+00 
0. OOOOOOE+OO 

44 0 . 165000E+00 
0. OOOOOOE+OO 

45 0 . 175000E+00 
0. OOOOOOE+OO 

46 0 . 185000E+00 
0. OOOOOOE+OO 

47 0 . 195000E+00 
0 .OOOOOOE+OO 

48 0 . 205000E+00 
0 .OOOOOOE+OO 

49 0 . 215000E+00 
0 .OOOOOOE+OO 

50 0 . 225000E+00 
0 .OOOOOOE+OO 

51 0 . 235000E+00 
0 .OOOOOOE + OO 

52 0 . 245000E+00 
0. OOOOOOE+OO 

53 0 . 255000E+00 
0 .OOOOOOE+OO 


0 . 112581E-02 
0 .OOOOOOE+OO 
0 . 109213E-02 
0 .OOOOOOE+OO 
0 . 106125E-02 
0 .OOOOOOE+OO 
0 . 103284E-02 
0. OOOOOOE+OO 
0 . 173760E-02 
0. OOOOOOE+OO 
0 .206746E-02 
0 .OOOOOOE+OO 
0 .227347E-02 
0 .OOOOOOE+OO 
0 . 241863E-02 
0 .OOOOOOE+OO 
0 . 252802E-02 
0 .OOOOOOE+OO 
0 .261392E-02 
0 .OOOOOOE+OO 
0 .268318E-02 
0 .OOOOOOE+OO 
0 .274005E-02 
0 .OOOOOOE+OO 
0 . 278733E-02 
0 .OOOOOOE+OO 
0 .282702E-02 
0 . OOOOOOE+OO 
0 . 2 95952E-02 
0. OOOOOOE+OO 
0 . 302481E-02 
0 .OOOOOOE+OO 
0 .305333E-02 
0 .OOOOOOE+OO 
0 . 306254E-02 
0 .OOOOOOE+OO 
0 .306189E-02 
0 .OOOOOOE+OO 
0 .305632E-02 
0 .OOOOOOE+OO 
0 .304839E-02 
0 .OOOOOOE+OO 
0 .303956E-02 
0 .OOOOOOE+OO 
0 .303049E-02 
0 .OOOOOOE+OO 
0 .302154E-02 
0 .OOOOOOE+OO 
0 .301291E-02 
0 .OOOOOOE+OO 
0 .300468E-02 
0. OOOOOOE+OO 
0 .299689E-02 
0. OOOOOOE+OO 
0 . 298954E-02 
0 .OOOOOOE+OO 


0 .726265E-03 
0 .OOOOOOE+OO 
0 . 7 48665E-03 
0 .OOOOOOE+OO 
0 .770447E-03 
0 .OOOOOOE+OO 
0 .791632E-03 
0 .OOOOOOE+OO 
0 .753936E-03 
0 .OOOOOOE+OO 
0 .718314E-03 
0 .OOOOOOE+OO 
0 . 690004E-03 
0 .OOOOOOE+OO 
0 .668187E-03 
0 .OOOOOOE+OO 
0 . 650370E-03 
0 . OOOOOOE+OO 
0 . 636766E-03 
0 . OOOOOOE+OO 
0 . 625154E-03 
0 .OOOOOOE+OO 
0 . 615407E-03 
0 .OOOOOOE+OO 
0 . 607552E-03 
0 .OOOOOOE+OO 
0 . 600424E-03 
0 .OOOOOOE+OO 
0 . 574461E-03 
0 .OOOOOOE+OO 
0 . 559251E-03 
0 .OOOOOOE+OO 
0 .552099E-03 
0 .OOOOOOE+OO 
0 .550301E-03 
0 .OOOOOOE+OO 
0 . 551563E-03 
0 .OOOOOOE+OO 
0 . 555210E-03 
0 .OOOOOOE+OO 
0 .559412E-03 
0 . OOOOOOE+OO 
0 . 564093E-03 
0 . OOOOOOE+OO 
0 . 569255E-03 
0 .OOOOOOE+OO 
0 . 574 115E-03 
0. OOOOOOE+OO 
0 . 578621E-03 
0 .OOOOOOE+OO 
0 . 583684E-03 
0 .OOOOOOE+OO 
0 . 588395E-03 
0 .OOOOOOE+OO 
0 .592743E-03 
0 .OOOOOOE + OO 


0 .562157E-03 
0 .316999E+03 
0 .579492E-03 
0 .316999E+03 
0 .596355E-03 
0 .316999E+03 
0 .612760E-03 
0 .316999E+03 
0 .573713E-03 
0 .316999E+03 
0 .539133E-03 
0 .316999E+03 
0 . 512221E-03 
0 .316999E+03 
0 . 4 91556E-03 
0 . 316999E+03 
0 . 475395E-03 
0.316999E+03 
0 .462452E-03 
0 . 316999E+03 
0 .451868E-03 
0 .316999E+03 
0 . 443074E-03 
0 . 316999E+03 
0 . 435683E-03 
0 . 316999E+03 
0 .429418E-03 
0 .316999E+03 
0 .407479E-03 
0 . 316999E+03 
0 .396181E-03 
0 .316999E+03 
0 .391620E-03 
0 .316999E+03 
0 .390896E-03 
0 .316999E+03 
0 . 392237E-0 3 
0 . 316999E+03 
0.394 64 7E-0 3 
0 . 316999E+03 
0 . 397593E-03 
0 . 316999E+03 
0 .400756E-03 
0 . 316999E+03 
0 . 403985E-03 
0 .316999E+03 
0 .407193E-03 
0 .316999E+03 
0 . 410333E-03 
0 .316999E+03 
0 .413379E-03 
0 .316999E+03 
0 .416318E-03 
0 .316999E+03 
0 . 419147E-03 
0 .316999E+03 
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0 . 299488E-04 
0 . 181473E+05 
0 .308726E-04 
0 . 17 604 4E + 0 5 
0 .317710E-04 
0 .171067E+05 
0 .326452E-04 
0 .166487E+05 
0 .319390E-04 
0 .268610E+05 
0 .306585E-04 
0 .317821E+05 
0 .294100E-04 
0 . 349898E+05 
0 .283353E-04 
0 . 373565E+05 
0 . 274318E-04 
0 .392156E+05 
0 .266692E-04 
0 .407277E+05 
0 .260194E-04 
0 . 419839E+05 
0 .254606E-04 
0 . 430422E+05 
0.249768E-04 
0 .439427E+05 
0 .245556E-04 
0 . 447144E+05 
0 .229890E-04 
0 .473867E+05 
0 .220893E-04 
0 .488197E+05 
0 . 216500E-04 
0 . 495480E+05 
0 .214881E-04 
0 . 498811E+05 
0 . 214830E-04 
0 . 499955E+05 
0 .215631E-04 
0 . 499912E+05 
0 .216884E-04 
0 . 499233E+05 
0 .218349E-04 
0 .498252E+05 
0 .219903E-04 
0 . 497136E+05 
0 .221477E-04 
0 . 495980E+05 
0 .223033E-04 
0 . 494836E+05 
0 .224548E-04 
0 .493732E+05 
0 .226011E-04 
0 .492682E+05 
0 .227417E-04 
0 .491690E+05 


i 


54 0.265000E+00 
0 . OOOOOOE+OO 

55 0 . 275000E+00 
0. OOOOOOE+OO 

56 0 . 285000E+00 
0. OOOOOOE+OO 

57 0 . 295000E+00 
0. OOOOOOE+OO 

58 0 . 305000E+00 
0. OOOOOOE+OO 


0 .298258E-02 
0 .OOOOOOE+OO 
0 .297607E-02 
0 .OOOOOOE+OO 
0 .296995E-02 
0 .OOOOOOE+OO 
0 .296419E-02 
0 .OOOOOOE+OO 
0 .295875E-02 
0 .OOOOOOE+OO 


0 . 596735E-03 
0. OOOOOOE+OO 
0 . 600944E-03 
0. OOOOOOE+OO 
0 . 605068E-03 
0. OOOOOOE+OO 
0 . 608888E-03 
0. OOOOOOE+OO 
0 . 612409E-03 
0. OOOOOOE+OO 


0 .421872E-03 
0 .316999E+03 
0 .424479E-03 
0 .316999E+03 
0 .426978E-03 
0 .316999E+03 
0 .429374E-03 
0 .316999E+03 
0 .431673E-03 
0 .316999E+03 


0 .228769E-04 
0 . 490751E+05 
0 .230059E-04 
0 . 489877E+05 
0 .231291E-04 
0 . 489059E+05 
0 .232468E-04 
0 . 488291E+05 
0 .233594E-04 
0 . 487570E+05 



7.4 Axisymmetric engine nacelle 


i x/c 

cf y 

1 0 . 373416E-02 
0 .OOOOOOE+OO 

2 0 . 218 87 5E-02 
0 . OOOOOOE+OO 

3 0. OOOOOOE+OO 
0 .000000E+00 

4 0 . 145997E-02 
0. OOOOOOE+OO 

5 0 . 561097E-02 
0 .OOOOOOE+OO 

6 0 . 111812E-01 
0. OOOOOOE+OO 

7 0 . 176662E-01 
0 .OOOOOOE+OO 

8 0 . 251550E-0 1 
0 .OOOOOOE+OO 

9 0 . 335804E-01 
0 .OOOOOOE+OO 

10 0 . 428933E-01 
0 .OOOOOOE+OO 

11 0 . 530432E-01 
0 .OOOOOOE+OO 

12 0 . 639961E-01 
0 .OOOOOOE+OO 

13 0.757222E-01 
0 .OOOOOOE+OO 

14 0 . 881880E-01 
0 . OOOOOOE+OO 

15 0.101362E+00 
0 .OOOOOOE+OO 

16 0 . 115218E+00 
0 . OOOOOOE+OO 

17 0 . 129732E+00 
0. OOOOOOE+OO 

18 0 . 144886E+00 
0 .OOOOOOE+OO 

19 0 . 160664E+00 
0 . OOOOOOE+OO 

20 0 . 177052E+00 
0 .OOOOOOE+OO 

21 0 . 194037E+00 
0 .OOOOOOE+OO 

22 0 . 211 606E+00 
0 .OOOOOOE+OO 

23 0 . 22974 9E + 00 
0 .OOOOOOE+OO 

24 0.248454E+00 
0 .OOOOOOE+OO 

25 0.267710E+00 
0 .OOOOOOE+OO 


c f x 
dspty 

0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 .162810E-01 
0. OOOOOOE+OO 
0 . 604084E-02 
0 .OOOOOOE+OO 
0 .367783E-02 
0 .OOOOOOE+OO 
0 .235163E-02 
0 .OOOOOOE+OO 
0 . 1254 68E-02 
0 .OOOOOOE+OO 
0 . 701414E-03 
0 . OOOOOOE+OO 
0 . 503271E-03 
0 .OOOOOOE+OO 
0 . 410599E-03 
0 . OOOOOOE+OO 
0 .389267E-03 
0 .OOOOOOE+OO 
0 .395935E-03 
0 . OOOOOOE+OO 
0 .392345E-03 
0 . OOOOOOE+OO 
0 .359495E-03 
0 .OOOOOOE+OO 
0 .312134E-03 
0 .OOOOOOE+OO 
0 . 276467E-03 
0 .OOOOOOE+OO 
0 .265231E-03 
0 .OOOOOOE+OO 
0 .273282E-03 
0 .OOOOOOE+OO 
0 .289138E-03 
0 .OOOOOOE+OO 
0 .302779E-03 
0 .OOOOOOE + OO 
0 .308940E-03 
0 .OOOOOOE+OO 
0 .307959E-03 
0 .OOOOOOE+OO 
0 .302694E-03 
0 .OOOOOOE+OO 
0 .295679E-03 
0 .OOOOOOE+OO 
0 .288261E-03 
0 .OOOOOOE+OO 
0 .280948E-03 
0 .OOOOOOE+OO 


blth 

thmoy 

0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 . 649704E-03 
0 .OOOOOOE+OO 
0 .778570E-03 
0 .OOOOOOE+OO 
0 . 951187E-03 
0 .OOOOOOE+OO 
0 . 130 662E-02 
0 .OOOOOOE+OO 
0 . 174349E-02 
0 .OOOOOOE+OO 
0 .217831E-02 
0 . OOOOOOE+OO 
0 .257521E-02 
0 . OOOOOOE+OO 
0 . 294645E-02 
0 .OOOOOOE+OO 
0 . 327812E-02 
0 .OOOOOOE+OO 
0.357713E-02 
0 .OOOOOOE+OO 
0 . 387022E-02 
0 .OOOOOOE+OO 
0 . 418669E-02 
0 .OOOOOOE+OO 
0 . 452897E-02 
0 .OOOOOOE+OO 
0 . 486913E-02 
0 .OOOOOOE+OO 
0 . 517285E-02 
0 .OOOOOOE+OO 
0 . 543096E-02 
0 .OOOOOOE+OO 
0 . 565398E-02 
0 .OOOOOOE+OO 
0 . 585929E-02 
0 .OOOOOOE+OO 
0 . 606272E-02 
0 .OOOOOOE+OO 
0 . 627200E-02 
0 . OOOOOOE+OO 
0 . 648871E-02 
0 . OOOOOOE+OO 
0 . 670 995E-02 
0. OOOOOOE+OO 
0 . 693435E-02 
0. OOOOOOE+OO 
0 . 715934E-02 
0. OOOOOOE+OO 
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dsptx 

twall 

0 .OOOOOOE+OO 
0. OOOOOOE+OO 
0 . 164296E-03 
0 . 446464E+03 - 
0 .201441E-03 
0 . 443914E+03 - 
0 . 263736E-03 
0 . 438654E+03 - 
0 . 398167E-03 
0 . 434196E+03 - 
0 . 598226E-03 
0 . 433090E+03 - 
0 .823670E-03 
0 . 433571E+03 
0 . 102513E-02 
0 . 434200E+03 
0 . 120198E-02 
0 .434770E+03 
0 . 133903E-02 
0 .435183E+03 
0 . 144198E-02 
0 .435428E+03 
0 .153890E-02 
0 .435557E+03 
0 .165871E-02 
0 . 435659E+03 
0 .180597E-02 
0 .435793E+03 
0 .195564E-02 
0 .435946E+03 
0 .207565E-02 
0 .436074E+03 
0 .215548E-02 
0.436149E+03 
0 .220556E-02 
0 . 436161E+03 - 
0 . 224417E-02 
0 . 436118E+03 - 
0 . 228 630E-02 
0 . 436040E+03 - 
0 . 233826E-02 
0.435946E+03 - 
0 . 239965E-02 
0 . 435851E+03 - 
0 .246756E-02 
0 . 435760E + 03 - 
0 .253932E-02 
0 . 435672E+03 - 
0 .261331E-02 
0.435588E+03 - 


thmox 

qw 

0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 .727296E-04 
0 . 193601E+00 
0 .856589E-04 
0 . 109754E+01 
0 .100272E-03 
0 . 212606E+01 
0 . 134 613E-03 
0 . 112155E+01 
0 . 185629E-03 
0 . 349651E-01 
0 .239134E-03 
0 . 867581E-01 
0 .287813E-03 
0 . 584770E-01 
0 . 334 112E-03 
0 .426648E-01 
0 .377128E-03 
0 .279242E-01 
0 .416248E-03 
0 . 151245E-01 
0 .452737E-03 
0 . 904173E-02 
0 .489669E-03 
0 . 130545E-01 
0 .528775E-03 
0 . 182800E-01 
0 . 568240E-03 
0 . 181670E-01 
0 .605288E-03 
0 .136696E-01 
0 . 638392E-03 
0 . 669523E-02 
0 . 667044E-03 
0 .144470E-02 
0 . 691956E-03 
0 . 864108E-02 
0 .714780E-03 
0 .129254E-01 
0 .737101E-03 
0 . 144695E-01 
0 . 759753E-03 
0 . 146568E-01 
0 . 7 82 92 IE-03 
0 . 145427E-01 
0 .806510E-03 
0 .144964E-01 
0 .830392E-03 
0 . 145445E-01 


I 


26 0 . 287509E+00 
0 .OOOOOOE+OO 

27 0 . 307 84 1E+00 
0 .000000E+00 

28 0 . 328698E+00 
0 .OOOOOOE+OO 

29 0.350070E+00 
0 .OOOOOOE+OO 

30 0 . 371952E+00 
0 .OOOOOOE+OO 

31 0 . 394334E+0 0 
0 .OOOOOOE+OO 

32 0 . 417210E+00 
0 .OOOOOOE+OO 

33 0 . 440574E+00 
0 .OOOOOOE + OO 

34 0 . 464419E + 00 
0 .OOOOOOE+OO 

35 0 . 488738E+00 
0 .OOOOOOE + OO 

36 0 . 513524E+00 
0 .OOOOOOE+OO 

37 0 . 538772E+00 
0 .OOOOOOE+OO 

38 0.564476E+00 
0. OOOOOOE + OO 

39 0 . 590632E+00 
0 .OOOOOOE+OO 

40 0 . 617239E+00 
0. OOOOOOE + OO 

41 0 . 644292E+00 
0 .OOOOOOE + OO 


0 . 273879E-03 
0 . OOOOOOE+OO 
0 .267085E-03 
0 .OOOOOOE+OO 
0 .260551E-03 
0 .OOOOOOE+OO 
0 . 254176E-03 
0 .OOOOOOE+OO 
0 . 2 4 771 9E-03 
0 .OOOOOOE+OO 
0 .2407 4 9E-0 3 
0 .OOOOOOE+OO 
0 .232856E-03 
0 .OOOOOOE+OO 
0 . 224 17 8E-03 
0 .OOOOOOE+OO 
0 . 216113E-03 
0 .OOOOOOE+OO 
0 . 211817E-03 
0 .OOOOOOE+OO 
0 . 215325E-03 
0 .OOOOOOE+OO 
0 . 227689E-03 
0 .OOOOOOE+OO 
0 .240595E-03 
0 .OOOOOOE+OO 
0 . 2334 91E-03 
0 . OOOOOOE+OO 
0 . 173675E-03 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 


0 .738459E-02 
0 .OOOOOOE+OO 
0 .761104E-02 
0 .OOOOOOE+OO 
0 .783843E-02 
0 .OOOOOOE+OO 
0 . 806708E-02 
0 .OOOOOOE+OO 
0 . 829838E-02 
0. OOOOOOE+OO 
0 .853478E-02 
0 .OOOOOOE+OO 
0 . 877912E-02 
0 .OOOOOOE+OO 
0 . 903265E-02 
0. OOOOOOE+OO 
0 . 929317E-02 
0 .OOOOOOE+OO 
0 . 954 539E-02 
0 . OOOOOOE+OO 
0 . 976756E-02 
0 .OOOOOOE+OO 
0 . 994758E-02 
0 .OOOOOOE+OO 
0 .101185E-01 
0 . OOOOOOE+OO 
0 . 103707E-01 
0 .OOOOOOE+OO 
0 . 109444E-01 
0 .OOOOOOE+OO 
0 . OOOOOOE+OO 
0 .OOOOOOE+OO 


0 .268874E-02 
0 . 435507E+03 
0 .276523E-02 
0 . 435427E+03 
0 .284269E-02 
0 .435349E+03 
0 . 292135E-02 
0 . 435272E+03 
0 .300203E-02 
0 .435197E+03 
0 .308644E-02 
0 . 435126E+03 
0 .317678E-02 
0 .435062E+03 
0 .327404E-02 
0 .435006E+03 
0 .337450E-02 
0 .434959E+03 
0 .346570E-02 
0 .434909E+03 
0 . 352770E-02 
0 .434839E+03 
0 .354934E-02 
0 . 434 723E+03 
0 . 355615E-02 
0 .434554E+03 
0 .362707E-02 
0 .434386E+03 
0 .394885E-02 
0 .434385E+03 
0 .OOOOOOE+OO 
0 .OOOOOOE+OO 


0 .854485E-03 
-0 . 146367E-01 
0 .878752E-03 
-0 . 147443E-01 
0 . 903194E-03 
-0 . 148422E-01 
0 . 927858E-03 
-0 . 148569E-01 
0 . 952871E-03 
-0 . 146227E-01 
0 . 978486E-03 
-0 .138521E-01 
0 . 100505E-02 
-0 . 123398E-01 
0 . 103282E-02 
-0 . 103677E-01 
0 .106154E-02 
-0 . 909932E-02 
0 . 108994E-02 
-0 . 107338E-01 
0 . 111567E-02 
-0 . 181953E-01 
0 .113653E-02 
-0 .331862E-01 
0 .115374E-02 
-0 . 481426E-01 
0 . 117656E-02 
-0 . 374389E-01 
0 . 123058E-02 
0 . 149628E-01 
0 .OOOOOOE+OO 
0 . OOOOOOE+OO 
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7.5 Ellipsoid of revolution along windward plane of symmetry 


i x/c cfx blth dsptx thmox 

cfy dspty thmoy twall qw 

1 0.400000E-02 0 . OOOOOOE+OO 0.000000E+00 0.000000E+00 0 . OOOOOOE+OO 

O.OOOOOOE+OO 0. OOOOOOE+OO 0. 000000E+00 0. 000000E+00 0. 000000E+00 

2 0.450000E-02 0.762597E-01 0.115106E-02 0.294283E-03 0.130723E-03 

O.OOOOOOE+OO 0.000000E+00 0.000000E+00 0.520259E+03 -0 . 274508E-02 

3 0 . 50000 0E-02 0.728950E-01 0.106555E-02 0.273060E-03 0.120332E-03 

O.OOOOOOE+OO 0.000000E+00 0.000000E+00 0.520259E+03 -0 . 335102E-02 

4 0 . 550000E-02 0.680647E-01 0.101057E-02 0.264054E-03 0.115054E-03 

O.OOOOOOE+OO 0.000000E+00 O.OOOOOOE+OO 0.520258E+03 -0 . 377599E-02 

5 0 . 600000E-02 0.634669E-01 0.984121E-03 0.260747E-03 0.112725E-03 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 0 . 520258E+03 -0 . 414369E-02 

6 0 . 800000E-02 0.502922E-01 0.951446E-03 0.260638E-03 0.111134E-03 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 0.520256E+03 -0 . 535532E-02 

7 0 . lOOOOOE-Ol 0 . 420724E-01 0.965137E-03 0.267585E-03 0.113716E-03 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 0.520255E+03 -0 . 624270E-02 

8 0 . 120000E-01 0 . 364407E-01 0.995366E-03 0.276117E-03 0.117206E-03 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 0.520254E+03 -0 . 687253E-02 

9 0 . 140000E-01 0 . 323392E-01 0.102208E-02 0.284682E-03 0.120716E-03 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 0.520253E+03 -0 . 731235E-02 

10 0 . 160000E-01 0 . 292091E-01 0.105039E-02 0.292961E-03 0.124077E-03 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 0.520252E+03 -0 . 760958E-02 

11 0 . 180000E-01 0 . 267302E-01 0.107384E-02 0.300960E-03 0.127294E-03 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 0.520251E+03 -0 . 779689E-02 

12 0 . 200000E-01 0 . 247096E-01 0.110036E-02 0.308734E-03 0.130398E-03 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 0.520250E+03 -0 . 789926E-02 

13 0 . 220000E-01 0.230253E-01 0.112447E-02 0.316328E-03 0.133413E-03 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 0.520249E+03 -0 . 793640E-02 

14 0.240000E-01 0.215956E-01 0.114627E-02 0.323768E-03 0.136357E-03 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 0.520248E+03 -0 . 792337E-02 

15 0.260000E-01 0.203640E-01 0.116935E-02 0.331073E-03 0.139239E-03 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 0.520247E+03 -0 . 787136E-02 

16 0.280000E-01 0.192900E-01 0.119281E-02 0.338258E-03 0.142065E-03 

O.OOOOOOE + OO O.OOOOOOE + OO O.OOOOOOE+OO 0 . 52024 6E+03 -0 . 7 78 97 6E-02 

17 0 . 30000 OE-0 1 0 . 183441E-01 0.121460E-02 0.345326E-03 0.144841E-03 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 0.520245E+03 -0 . 7 68674E-02 

18 0.320000E-01 0.175036E-01 0.123509E-02 0.352284E-03 0.147568E-03 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 0.520244E+03 -0 . 756745E-02 

19 0 . 34000 OE-0 1 0.167511E-01 0.125490E-02 0.359137E-03 0.150251E-03 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 0.520244E+03 -0 . 743644E-02 

20 0 . 36000 OE-0 1 0.161442E-01 0.127630E-02 0.365285E-03 0.152765E-03 

O.OOOOOOE+OO O.OOOOOOE+OO O.OOOOOOE+OO 0 . 520243E+03 -0 . 750395E-02 


SO 


7.6 Ellipsoid of revolution along leeward plane of symmetry 


i 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 


x/c 
cf y 

0 .400000E-02 
0 .000000E+00 
0 .450000E-02 
0 . 000000E+00 
0 .500000E-02 
0 .OOOOOOE+OO 
0 .550000E-02 
0 .OOOOOOE+OO 
0 .600000E-02 
0 .OOOOOOE+OO 
0 . 800000E-02 
0 .000000E+00 
0 .100000E-01 
0 .000000E+00 
0 . 120000E-01 
0. 000000E+00 
0 . 140000E-01 
0. OOOOOOE+OO 
0 . 160000E-01 
0 .OOOOOOE+OO 
0 .180000E-01 
0 .OOOOOOE+OO 
0 .200000E-01 
0 .OOOOOOE+OO 
0 .220000E-01 
0. OOOOOOE+OO 
0 . 240000E-01 
0 .OOOOOOE+OO 
0 .260000E-01 
0 .OOOOOOE+OO 
0 . 280000E-01 
0. OOOOOOE+OO 
0.300000E-01 
0. OOOOOOE+OO 
0 . 32000 OE-0 1 
0. OOOOOOE+OO 
0 . 340000E-01 
0. OOOOOOE+OO 
0 . 360000E-01 
0 .OOOOOOE+OO 


cf x 
dspty 

0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 . 321218E-01 
0 .OOOOOOE+OO 
0 . 292291E-01 
0 .OOOOOOE+OO 
0 .274982E-01 
0 .OOOOOOE+OO 
0 .261933E-01 
0 .OOOOOOE+OO 
0 .225559E-01 
0 .OOOOOOE+OO 
0 .201938E-01 
0 .OOOOOOE+OO 
0 .184437E-01 
0 .OOOOOOE+OO 
0 .170506E-01 
0 .OOOOOOE+OO 
0 .158981E-01 
0 .OOOOOOE+OO 
0 . 149219E-01 
0 .OOOOOOE+OO 
0 . 140812E-01 
0 .OOOOOOE+OO 
0 .133477E-01 
0 .OOOOOOE+OO 
0 . 127009E-01 
0 .OOOOOOE+OO 
0 . 121253E-01 
0. OOOOOOE+OO 
0 . 116091E-01 
0 .OOOOOOE+OO 
0 . 11142 9E-01 
0 .OOOOOOE+OO 
0 .107194E-01 
0 .OOOOOOE+OO 
0 . 10332 6E-01 
0 .OOOOOOE+OO 
0 .100190E-01 
0 .OOOOOOE+OO 


blth 

thmoy 

0. OOOOOOE+OO 
0. OOOOOOE+OO 
0 . 784167E-03 
0 .OOOOOOE+OO 
0 . 821668E-03 
0. OOOOOOE+OO 
0 . 850371E-03 
0. OOOOOOE+OO 
0 .875234E-03 
0. OOOOOOE+OO 
0 . 945531E-03 
0. OOOOOOE+OO 
0 . 996408E-03 
0. OOOOOOE+OO 
0 . 103881E-02 
0 .OOOOOOE+OO 
0 .107739E-02 
0 .OOOOOOE+OO 
0 . 111213E-02 
0 .OOOOOOE+OO 
0 . 114456E-02 
0 .OOOOOOE+OO 
0 . 117782E-02 
0 .OOOOOOE+OO 
0 . 120929E-02 
0 .OOOOOOE+OO 
0 . 123931E-02 
0 .OOOOOOE+OO 
0 . 126807E-02 
0 .OOOOOOE+OO 
0 . 129757E-02 
0 .OOOOOOE+OO 
0 . 132689E-02 
0 .OOOOOOE+OO 
0 . 135515E-02 
0 .OOOOOOE+OO 
0 .138247E-02 
0 .OOOOOOE+OO 
0 . 140831E-02 
0 .OOOOOOE+OO 


dsptx 

twall 

0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 .217507E-03 
0 .520253E+03 
0 .231428E-03 
0 .520250E+03 
0 .240952E-03 
0 .520248E+03 
0 .248208E-03 
0 .520247E+03 
0 .268665E-03 
0 . 520243E+03 
0 .283017E-03 
0 .520241E+03 
0 .294983E-03 
0 .520239E+03 
0 .306002E-03 
0 . 520238E+03 
0 .31657 6E-03 
0 . 520236E+03 
0 .326881E-03 
0 .520235E+03 
0 .336977E-03 
0 .520233E+03 
0 .346896E-03 
0 .520232E+03 
0 .356655E-03 
0 .520231E+03 
0 .366265E-03 
0 .520230E+03 
0 .375739E-03 
0 .520229E+03 
0 . 385087E-03 
0 . 520229E+03 
0 .394315E-03 
0 . 520228E+03 
0 . 403430E-03 
0 . 520227E+03 
0 . 411904E-03 
0 . 520227E+03 


thmox 

qw 

0 .OOOOOOE+OO 
0 .OOOOOOE+OO 
0 . 918 1 93E-0 4 
-0 .876983E-02 
0 .970987E-04 
-0 .839545E-02 
0 . 101107E-03 
-0 .850976E-02 
0 . 104313E-03 
-0 .864227E-02 
0 . 113260E-03 
-0 . 905064E-02 
0 . 119414E-03 
-0 .926889E-02 
0 . 124 4 14E-03 
-0 . 929852E-02 
0 .128910E-03 
-0 . 918767E-02 
0 . 133153E-03 
-0 . 898171E-02 
0 .137243E-03 
-0 . 871521E-02 
0 . 141224E-03 
-0 . 841304E-02 
0 . 145116E-03 
-0 . 809191E-02 
0 . 148931E-03 
-0 . 776332E-02 
0 . 152677E-03 
-0 . 743518E-02 
0 .156361E-03 
-0 . 711269E-02 
0 .159989E-03 
-0 . 67 9867E-02 
0 . 1 63563E-03 
-0 . 64 9590E-02 
0 . 167088E-03 
-0 . 620 523E-02 
0 . 170474E-03 
-0 . 612996E-02 
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7.7 Infinite swept, cyl icier along the leading-edge attachment line 


i x/c 

cf y 

1 0 . 100000E-0 1 
0 .000000E+00 

2 0 . 2 00 00 OE-0 1 
0 .OOOOOOE+OO 

3 0.300000E-01 
0. 000000E+00 

4 0 . 400000E-01 
0. OOOOOOE+OO 


cfx blth dsptx thmox 

dspty thmoy twall qw 


0 .227267E-02 
0 .000000E+00 
0 .227267E-02 
0. 000000E+00 
0 .227267E-02 
0. 000000E+00 
0 .227267E-02 
0 .000000E+00 


0 .212016E-02 
0 .OOOOOOE+OO 
0 . 212016E-02 
0. 000000E+00 
0.212016E-02 
0 .000000E+00 
0 .212016E-02 
0 .OOOOOOE+OO 


0 . 663990E-03 
0 . 521493E+03 
0 . 663990E-03 
0 . 521493E+03 
0 . 663990E-03 
0 . 521493E+03 
0.663990E-03 
0 . 521493E+03 


0 .260441E-03 
-0 .888370E-08 
0 .26044 IE-03 
0 .222093E-07 
0 . 260441E-03 
0 .888370E-08 
0 .26044 1E-0 3 
0 . 177674E-07 


7.8 Swept-back wing 


i x/c 

cf y 

1 0 . 121900E-03 
0 .OOOOOOE+OO 

2 0 . 864300E-03 
0 . 23694 5E-0 1 

3 0 . 160670E-02 
0 . 150014E-01 

4 0 . 234 910E-02 
0.114673E-01 

5 0 . 309160E-02 
0 . 940938E-02 

6 0.383400E-02 
0 . 809523E-02 

7 0 . 606120E-02 
0 . 601189E-02 

8 0.828850E-02 
0 . 490637E-02 

9 0 . 105158E-01 
0 . 411003E-02 

10 0 . 127431E-01 
0 .355803E-02 

11 0.149703E-01 
0 .321903E-02 

12 0.17197 6E-01 
0 .296348E-02 

13 0 . 194249E-01 
0 .271542E-02 

14 0 . 246218E-01 
0 .224944E-02 

15 0 . 298188E-01 
0 .198989E-02 

16 0.350158E-01 
0 . 181834E-02 

17 0 . 402127E-01 
0 .163004E-02 

18 0 . 454097E-01 
0 . 145231E-02 

19 0 . 506067E-01 
0 .183446E-02 

20 0.558036E-01 
0 .191825E-02 

21 0 . 669400E-01 
0 .189381E-02 

22 0 . 780763E-01 
0 . 182969E-02 

23 0 . 892127E-01 
0 . 17791 4E-02 

24 0 . 10034 9E + 0 0 
0 . 173047E-02 

25 0 . 111485E+00 
0 .168803E-02 


cf x 
dspty 

0 . 648001E-02 
0 . OOOOOOE+OO 
0 . 619005E-02 
0 . 103222E-03 
0 .594413E-02 
0 . 110595E-03 
0 .573117E-02 
0 . 117804E-03 
0 .551083E-02 
0 .125769E-03 
0 . 531356E-02 
0 .133637E-03 
0 . 488254E-02 
0 .153936E-03 
0 . 456890E-02 
0 . 172218E-03 
0 . 428030E-02 
0 . 191890E-03 
0 . 402506E-02 
0 .211491E-03 
0 .384535E-02 
0 .228496E-03 
0 .370054E-02 
0 .243966E-03 
0 .355980E-02 
0 .260184E-03 
0 .324529E-02 
0 .299421E-03 
0 .303573E-02 
0 .333533E-03 
0 .288465E-02 
0 .363505E-03 
0 . 273238E-02 
0 .396591E-03 
0 .257720E-02 
0 . 431381E-03 
0 .353363E-02 
0 .418741E-03 
0 .393445E-02 
0 .413101E-03 
0 .432097E-02 
0 .419399E-03 
0 .450060E-02 
0 .433497E-03 
0 . 462176E-02 
0 . 448325E-03 
0 .471450E-02 
0 . 463284E-03 
0 . 478874E-02 
0 . 476993E-03 


blth 

thmoy 

0 . 479884E-03 
0 .000000E+00 
0 .494110E-03 
0 .445461E-04 
0 .508753E-03 
0 .471944E-04 
0 .524839E-03 
0 .497426E-04 
0 .542727E-03 
0 .525154E-04 
0 .563707E-03 
0 . 552665E-04 
0 . 617817E-03 
0 .623367E-04 
0 . 670 657E-03 
0 . 685211E-04 
0.726227E-03 
0 . 749336E-04 
0 .780888E-03 
0 . 813126E-04 
0 .829711E-03 
0 .869714E-04 
0 . 874538E-03 
0 . 920822E-04 
0 . 918225E-03 
0 . 971654E-04 
0 . 102787E-02 
0 . 109347E-03 
0 .112438E-02 
0 .120144E-03 
0 . 121123E-02 
0 . 129691E-03 
0 . 130095E-02 
0 . 139528E-03 
0 . 139031E-02 
0 . 149611E-03 
0 . 140956E-02 
0 . 1 51 633E-03 
0 . 143290E-02 
0 . 151966E-03 
0 . 1 50220E-02 
0 . 154702E-03 
0 . 157116E-02 
0 . 159095E-03 
0 .163121E-02 
0.16378 5E-03 
0 . 1 69362E-02 
0 . 168500E-03 
0 . 174453E-02 
0 . 172948E-03 


dsptx 

twall 

0 . 152013E-03 
0 . 483938E+03 - 
0 .160975E-03 
0 . 483524E+03 - 
0 . 170342E-03 
0 . 483074E+03 - 
0 .179489E-03 
0 . 482617E+03 - 
0 . 189243E-03 
0 . 482191E+03 - 
0 . 198824E-03 
0 . 481805E+03 - 
0 .224182E-03 
0 . 480814E+03 - 
0 .247221E-03 
0 . 479963E+03 - 
0 .270905E-03 
0 . 479238E+03 - 
0 .293860E-03 
0 . 478655E+03 - 
0 .314394E-03 
0.478138E+03 - 
0 .333489E-03 
0 . 477658E+03 - 
0 .352891E-03 
0 . 477223E+03 - 
0 .397694E-03 
0 . 476467E+03 - 
0.437282E-03 
0 . 475862E+03 - 
0 . 473053E-03 
0 . 475315E+03 - 
0 .510916E-03 
0 . 474844E+03 - 
0 . 548686E-03 
0 . 474499E+03 - 
0 . 534 318E-03 
0 . 474589E+03 - 
0 .525487E-03 
0 . 474764E+03 - 
0 .523764E-03 
0 . 4 75206E+03 - 
0 .530200E-03 
0 . 475667E + 03 - 
0 .538715E-03 
0 . 476083E+03 - 
0 . 548166E-03 
0 . 476451E+03 - 
0 .556812E-03 
0.476799E+03 - 


thmox 

qw 

0 . 585255E-04 
0 .351638E-04 
0 . 617156E-04 
0 . 729616E-01 
0 . 649066E-04 
0 . 132228E+00 
0 . 679731E-04 
0 .171745E+00 
0 . 712194E-04 
0 . 184495E+00 
0 .744109E-04 
0 .187440E+00 
0 .827832E-04 
0 . 188232E+00 
0 . 902216E-04 
0 . 178061E+00 
0 . 97727 6E-04 
0 . 147379E+00 
0 . 105010E-03 
0 . 119317E+00 
0 . 111547E-03 
0 . 111213E+00 
0 . 117558E-03 
0 . 104321E+00 
0 . 123488E-03 
0 .895337E-01 
0 . 137203E-03 
0 . 578648E-01 
0 .149344E-03 
0 .489781E-01 
0 . 1 60250E-03 
0 .443263E-01 
0 . 171375E-03 
0 .330507E-01 
0 . 182438E-03 
0 .222588E-01 
0 . 187133E-03 
0 . 111386E+00 
0 . 189579E-03 
0 .123242E+00 
0 . 194247E-03 
0 . 121451E+00 
0 . 198671E-03 
0 . 116441E+00 
0 .202695E-03 
0 .114545E+00 
0 . 206462E-03 
0 . 111563E+00 
0 .209783E-03 
0 . 108508E+00 
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26 0.122622E+00 
0 . 1 66666E-02 

27 0.133758E+00 
0.16548 8E-02 

28 0 . 144894E+00 
0 . 164594E-02 

29 0 . 165682E+00 
0 . 163740E-02 

30 0 . 186470E + 00 
0 . 162572E-02 

31 0 . 207258E+00 
0 .160746E-02 

32 0 . 228 04 6E+00 
0 . 158923E-02 

33 0 . 248834E+00 
0 .157469E-02 

34 0 . 269622E+00 
0 .156244E-02 

35 0 . 290409E+00 
0 . 155357E-02 

36 0 . 311197E+00 
0 . 154753E-02 

37 0 . 331985E+00 
0 .154312E-02 

38 0 . 352773E+00 
0 . 860034E-03 

39 0 . 378758E+00 
0 . 676153E-03 

40 0.404743E+00 
0 .600065E-03 

41 0 . 430727E+00 
0 .549513E-03 

42 0 . 456712E+00 
0 . 498313E-03 

43 0 . 482697E+00 
0 . 442710E-03 

44 0 . 508682E+00 
0 .365094E-03 

45 0.523070E+00 
0 . 318873E-03 

46 0 . 539121E+00 
0 .285474E-03 

47 0 . 569561E+00 
0 .244581E-03 

48 0 . 600000E+00 
0 .215248E-03 


0 .487009E-02 
0 .488101E-03 
0 .495858E-02 
0 . 497413E-03 
0 . 504 627E-02 
0 . 505215E-03 
0 . 522451E-02 
0 . 517367E-03 
0 . 541044E-02 
0 . 528528E-03 
0 .559782E-02 
0 . 539779E-03 
0 . 578921E-02 
0 . 550490E-03 
0 .599498E-02 
0 . 560698E-03 
0 . 62134 6E-02 
0 . 570148E-03 
0 . 644 903E-02 
0 . 578784E-03 
0 . 670440E-02 
0 . 586686E-03 
0 . 698122E-02 
0.594063E-03 
0 . 394345E-02 
0 . 767111E-03 
0 .313802E-02 
0 . 911435E-03 
0 .285078E-02 
0 . 102506E-02 
0 .270359E-02 
0 .112561E-02 
0 .260168E-02 
0 .122766E-02 
0 .251322E-02 
0 . 133404E-02 
0 . 2 39502E-02 
0 . 147390E-02 
0 .232022E-02 
0 .155986E-02 
0 .227911E-02 
0 .164967E-02 
0 .226365E-02 
0 .180120E-02 
0 .229420E-02 
0 . 194051E-02 


0 .179429E-02 
0 . 1 767 65E-03 
0 . 183724E-02 
0 .179963E-03 
0 . 187373E-02 
0 . 182584E-03 
0 . 193499E-02 
0 . 186338E-03 
0 . 1 99060E-02 
0 . 189273E-03 
0 . 204009E-02 
0 . 1 91 920E-03 
0 . 208417E-02 
0 . 194349E-03 
0 .212460E-02 
0 . 196592E-03 
0 . 216111E-02 
0 . 198581E-03 
0 .219391E-02 
0 .200293E-03 
0 .222294E-02 
0 .201721E-03 
0 .224805E-02 
0 .202905E-03 
0 . 254270E-02 
0 .240279E-03 
0 .279693E-02 
0 .273455E-03 
0 . 302247E-02 
0 .302347E-03 
0 .323873E-02 
0 .328922E-03 
0 . 345275E-02 
0 .354951E-03 
0 .367205E-02 
0 .380867E-03 
0 . 3 9364 6E-02 
0 . 4 10803E-03 
0 . 407961E-02 
0 . 427664E-03 
0 . 423750E-02 
0 . 446153E-03 
0 . 451179E-02 
0 . 479014E-03 
0 . 475976E-02 
0 .510110E-03 


0 .564001E-03 
0 . 477103E+03 
0 . 570528E-03 
0 .477356E + 03 
0 .576336E-03 
0 .477570E+03 
0.586482E-03 
0 .477881E + 03 
0.596669E-03 
0.478120E+03 
0 . 606964E-03 
0 . 478334E+03 
0 . 616575E-03 
0 . 478548E+03 
0 . 62577 6E-03 
0 .478753E + 03 
0 . 63431 9E-0 3 
0 .478949E+03 
0 . 642279E-03 
0 . 479133E+03 
0 . 649792E-03 
0 . 479302E+03 
0 . 657027E-03 
0 . 479456E+03 
0 . 837342E-03 
0 . 478408E+03 
0 . 994587E-03 
0.476932E+03 
0 . 112313E-02 
0.475778E+03 
0 .123934E-02 
0 .474868E+03 
0 . 135522E-02 
0 .474151E+03 
0 . 147118E-02 
0 .473632E+03 
0 . 161111E-02 
0 . 473351E+03 
0 . 169164E-02 
0 . 4 73302E + 03 
0 . 177448E-02 
0 . 473282E+03 
0 . 191354E-02 
0 . 473267E+03 
0 .204084E-02 
0 .473277E+03 


0 .212525E-03 
-0 . 109966E+00 
0 .214842E-03 
-0 . 113198E+00 
0 . 216784E-03 
-0 . 116654E + 00 
0 . 2 19823E-03 
-0 . 125483E+00 
0 .222523E-03 
-0 . 133606E+00 
0 .225085E-03 
-0 .139685E+00 
0 .227413E-03 
-0 . 145350E+00 
0 .229560E-03 
-0 . 152037E+00 
0 .231463E-03 
-0 .159236E+00 
0 .233141E-03 
-0 . 167645E+00 
0 .234626E-03 
-0 . 177352E+00 
0 .235965E-03 
-0 . 188119E+00 
0 .272363E-03 
-0 . 670153E-02 
0 .306064E-03 
-0 . 559573E-02 
0 .336483E-03 
-0 . 507192E-02 
0 .365323E-03 
-0 .452597E-02 
0 .394160E-03 
-0 .323506E-02 
0 . 423057E-03 
-0 .164004E-02 
0 .456416E-03 
0 . 642064E-03 
0 .475463E-03 
0 .150493E-02 
0 .496192E-03 
0 .151861E-02 
0 .532796E-03 
0 .137205E-02 
0 .567603E-03 
0 . 12832 8E-02 
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Axisymmetric Flow 


Figure 1. Coordinate System lor 2-D and Axisymmetric Flow 
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Figure 2. Coordinate System lor the Flow Along Plane of Symmetry 
and Along Leading-Edge Attachment Line 
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Figure 3. Coordinate System for the Swept-Back Wing Flow 
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Figure 4. Coordinate System for the Swept- Forward Wing Flow 
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Figure 5. Coordinate System for the Delta Wing Flow 
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Figure 6. Airfoil Definition for Swept Wing Flow 
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Figure 12. Skin Friction Coefficient (Swept Wing, = 0.80, without suction) 
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Figure 14. Comparison ol (. alciilatecl N 1 actor (Trcci — -3000 Hz) 
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Figure 16. Velocity Profile and Its Normal Derivatives 

(Swept Wing, Moo = 0.80, with suction) 

103 



DO LOOP for I=2,IMAX 

I ’ 





















REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No 0704-0188 

O, 0 ... ~ -r^t.or ■'* rr-a*:~r •<. M i.e’ ’ " .j’ DO’ ’OSDC”se r- Udi~c: tr** re.'?*'" - ‘ s 

lV ^ f ; 3 ’ , nr ,pc -Vu nppd*»o a-o >:r-rrDieii° z .--c - ** -d .--.lecTiO^ O? = cements r^aromg this o-jrdoi ost r- O' in. ;.y*pr ( l s 

UntiiJ ■ V— o.ro.r t- .v„n.n,to» L -tkmn ■ ' or o'-’'™’ 

n^vis a ■ Su ’Jrj ; i 22202-S3C2 ana t: t*? o*fi 'ira rtuiaje*. p 5 dp ' c ' Roau'tion P'Cje-.t (G7C**-0 ■ 36). Aasn. -..r - . 

1. AGENCY USE ONLY (Leave blank) 2. REPORT DATE 3. REPORT TYPE AN( 

December 1992 Contractor R( 

3 OATES COVERED 

iDort 

4. TITLE AND SUBTITLE 

BLSTA - A Boundary Layer Code for Stability Analysis 

5. FUNDING NUMBERS 

C NAS1-19299 
WU 505-59-10 

6. AUTHOR(S) 

Yong-Sun Wie 

7. PERFORMING ORGANIZATION NAME(S) AND ADORE SS(ES) 

High Technology Corporation 
28 Research Drive 
Hampton, VA 23666 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

HTC-9204 

9. SPONSORING /MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Langley Research Center 
Hampton, VA 23681 

10. SPONSORING MONITORING 
AGENCY REPORT NUMBER 

NASA CR-4481 

11. SUPPLEMENTARY NOTES 

Langley Technical Monitor: Fayette S. Collier, Jr. 

Final Report 

12a. DISTRIBUTION / AVAILABILITY STATEMENT 1 1ib DISTRIBUTION CODE 


Unclassified - Unlimited 
Subject Category 43" S*/ 


13. AB5TRACT (Maximum 200 words) 

A computer program is developed to solve the compressible, laminar boundary- 
layer equations for two- dimensional flow, axisymmetric flow, and quasi-three- 
dimensional flows including the flow along the plane of symmetry, flow along the 
leading-edge attachament line, and swept wing flows with a conical flow approximatiorj 
The finite-difference numerical procedure used to solve the governing equations is 
second-order accurate. The flow over a wide range of speed, from subsonic to 
hypersonic speed with perfect gas assumption, can be calculated. Various wall 
boundary conditions, such as wall suction or blowing and hot or cold walls, can be 
applied. 

The results indicate that this boundary-layer code gives velocity and tempera 
ture profiles which are accurate, smooth, and continuous through the first and 
second normal derivatives. The code presented herein can be coupled with a 
stability analysis code and used to predict the onset of the boundary-layer trans- 
ition which enables the assessment of the laminar flow control techniques. A 
User's Manual is also included. 


14. SUBJECT TERMS 


laminar boundary layer, linear stability analysis, 
conical flow approximation 


15. NUMBER OF PAGES 
116 


16. PRICE CODE 

A06 


17. SECURITY CLASSIFICATION 
OF REPORT 

Unclassified 


18. SECURITY CLASSIFICATION 
OF THIS PAGE 

Unclassified 


19 


SECURITY CLASSIFICATION 
OF ABSTRACT 


20. LIMITATION OF ABSTRACT 


NSN 7540-0 1 -280-5500 


Standard Form 298 (Rev 2-89) 

by ANSI Std 239-’b 

298-107 


NASA-Langli’y. 1992 














National Aeronautics and 
Space Administration 
Code JTT 
Washington, D.C. 
20546-0001 

Official Business 

Penalty for Private Use, $300 



NASA 


POSTMASTER ^ Undeliverable (Section 158 
Postal Manual) Do Not Return 



